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CHAPTER  1 


Introduction  and  Review 

1 . 1  Purpose  and  Scope 

Suppose  one  wishes  to  determine  the  exact  probability  density 
function  of  the  product  of  two  random  variables,  X  and  Y,  with  known 
probability  density  functions,  fx(x)  and  fY(y)»  such  that  fx(x)  -  0 
for  x  <  0  and  f^(y)  ■  0  for  y  <  0.  If  X  and  Y  are  independent,  the 
desired  answer  is  the  inverse  Mellin  transform  of  the  product  of  the 
Mellin  transforms  of  f^(x)  and  fY(y).  This  process  is  well 
established  and  has  been  used  extensively  in  the  algebra  of 
independent  variables. 

Now,  suppose  that  X  and  Y  are  not  independent.  If  the 
bivariate  probability  density  function  fx>Y(x,y)  is  known,  an  answer 
still  may  be  obtained  by  the  use  of  double  Mellin  transform 
techniques.  To  date,  very  little  work  has  been  done  in  this  area  due 
to  the  difficulties  of  performing  the  double  transform  operations. 

Suppose,  however,  that  one  has  a  general  function  of  two 
variables  which  has  as  special  cases  all  of  the  bivariate  probability 
density  functions  in  some  group  of  interest.  If  a  solution  exists  for 
the  general  function  of  two  variables  by  application  of  the  double 
transform  techniques,  then  the  resulting  solution  covers  all  those 
problems  involving  the  special  cases.  This  is  the  motivation  for 
using  a  general  function  of  two  variables. 

Suppose  further,  that  the  application  of  the  double  transform 


techniques  to  the  general  function  of  two  variables  has  as  its 
solution  a  general  function  of  one  variable  which  has  as  special  cases 
all  the  univariate  probability  density  functions  of  some  interest 
group;  specifically,  the  H-function  of  one  variable.  If  this 
supposition  is  true,  then  one  has  a  powerful  technique  for  finding 
exact  distributions  of  products  and  quotients  of  dependent  random 
variables,  and  a  means  of  combining  these  distributions  with  other 
independent  univariate  distributions. 

The  primary  purpose  of  this  dissertation  is  to  develop  a 
general  technique,  presented  in  Chapter  5,  for  determining  the 
probability  density  function  and  the  cumulative  density  function  of 
the  random  variable 

Z  -  XpY^ 

where  X  and  Y  are  dependent  random  variables  with  joint  probability 
density  function  which  may  be  expressed  as  a  bivariate  H-function  and 
p  and  q  are  rational  constants.  The  general  function  of  two  variables 
known  as  the  bivariate  H-function  is  chosen  for  several  reasons. 

First,  the  bivariate  H-function  is  the  most  general  of  the  special 
functions  of  two  variables  and  includes  nearly  every  named  function  as 
a  special  case.  Second,  much  like  the  bivariate  normal  arises  as  a 
natural  extension  to  the  univariate  normal,  it  would  seem  that  the 
bivariate  H-function  distribution  should  be  a  natural  extension  to  the 
univariate  H-function  distribution.  This  is  indeed  the  case.  Chapter 
4  shows  that  many  of  the  properties  that  hold  for  univariate 


H-f unction  distributions  hold  also  for  the  bivariate  H-function 
distribution.  The  bivariate  H-function  distribution  has  the 
additional  feature  that  the  products,  quotients  and  rational  powers  of 
dependent  K-function  variates  are  reduced  to  univariate  H-function 
variates . 

In  developing  the  theory  above,  a  second  purpose  evolved  as  a 
natural  extension  to  the  primary  purpose  above.  A  technique  is 
presented  in  Chapter  5  for  determining  the  bivariate  probability 
density  function  for  the  random  variables 

n  p  n  q. 

Z  «  il  X.  ,  W  >  II  Y. 
i-1  1  i-1 

where  each  pair,  (X^Y^,  are  dependent  variates  with  a  given 
bivariate  density  function  which  is  expressible  as  a  bivariate 
H-function  and  X^,  i^j  and  Y^,  Y^  i#j  are  independent.  Exponents 
p^  and  q^  are  rational  constants.  It  is  shown  in  Chapter  5  that  such 
combinations  of  pairwise  independent  variates  result  in  the  dependent 
variates  Z  and  W  which  have  a  bivariate  density  function  which  is  also 
expressible  as  a  bivariate  H-function  distribution. 

In  the  course  of  developing  the  above  general  technique,  some 
secondary  purposes  became  evident.  One  is  the  attempt,  in  Chapter  3, 
to  relate  bivariate  H-functions  to  known  special  functions  and  to 
other  simpler  H-functions.  The  general  form  for  the  bivariate 
H-function  is  a  double  contour  integral  containing  products  and 


quotients  of  gamma  functions  and  is  not  readily  identified  by  this 
form.  The  properties  and  identities  given  in  Chapter  3  prove  useful 
in  evaluating  H-functions  in  later  Chapters. 

In  his  dissertation,  Cook  (5)  gives  a  convergence  proof  for 
the  univariate  H-function  distribution  and  a  readily  applied  technique 
for  inverting  the  H-function  using  residue  theory.  Chapter  6 
demonstrates  the  applicability  of  using  residue  theory  in  the 
bivariate  case  and  several  examples  are  given.  While  the  techniques 
given  are  generally  applicable,  a  bivariate  proof  similar  to  that 
given  by  Cook  would  still  be  required  to  develop  practical  guidelines 
for  when  left  half  plane  residues  versus  right  half  plane  residues 
should  be  summed  in  each  contour  integral  to  evaluate  a  given 
bivariate  H-function. 

To  study  and  develop  the  power  of  the  bivariate  H-function 
distribution,  three  new  bivariate  exponential  distributions  are 
developed  in  Appendix  C.  These  distributions  are  shown  to  be  special 
cases  of  the  bivariate  H-function  distribution  in  Chapter  4.  These 
distributions  indicate  the  versatility  of  the  bivariate  H-function 
distribution  in  the  contours  the  bivariate  H-function  can  undertake  as 
given  in  Appendix  D. 

Some  important  limitations  to  the  scope  of  this  dissertation 
must  be  stated.  For  instance,  techniques  are  presented  for 
determining  products,  quotients  and  rational  powers  of  dependent 
H-function  variates.  The  study  of  sums  and  differences  of  independent 


H-function  variates  has  been  established  through  the  work  of  Carter 
(3,4),  Eldred  (7),  and  Cook  (5).  The  extension  to  dependent 
H-function  variates  for  sums  and  differences  is  not  immediately 
accomplished.  To  do  so,  one  must  first  extend  Prasad's  theorems  (57) 
for  obtaining  the  Mellin  transform  from  its  Laplace  transform  and 
conversely  to  the  bivariate  case.  Then,  one  must  be  able  to  obtain 
the  Laplace  transform  of  the  bivariate  H-function.  Some  work  in  this 
area  has  already  been  accomplished  by  Goyal  (69).  Finally,  one  must 
determine  if  a  general  technique  for  using  double  Laplace  transforms 
exists  for  solving  for  the  density  function  of  a  random  variable  which 
is  the  sum  of  two  dependent  random  variables.  While  it  would  seem 
natural  that  such  a  theory  similar  to  that  developed  for  products  and 
quotients  of  dependent  random  variables  using  double  Mellin  transforms 
exists,  the  development  of  such  a  theory  would  present  a  formidable 
task. 

The  bivariate  H-function  is  not  defined  for  a  zero  or  negative 
real  value  of  its  arguments.  Therefore,  only  probability  density 
functions  that  are  defined  to  be  zero  for  nonpositive  arguments  are 
treated.  Techniques  for  finding  probability  density  functions  defined 
non-zero  for  both  positive  and  negative  arguments  are  handled  by 
dividing  such  functions  into  four  components,  one  for  each  quadrant. 
Such  techniques  are  presented  by  Fox  (10)  and  Subrahamian  (19).  These 
techniques  result  as  an  extension  to  work  done  by  Epstein  (9)  and 
Springer  and  Thompson  (18)  for  independent  variables. 


The  algebra  of  random  variables  is  a  vast  field  of  study  and 
the  study  of  the  algebra  of  dependent  random  variables  is  still  in  a 
relative  state  of  infancy*  Combining  the  advantages  of  a  general 
function  and  of  certain  properties  of  the  bivariate  H-function  and  its 
subsequent  reduction  to  a  univariate  H-function  for  products  and 
quotients  of  dependent  H-function  variates  is,  hopefully,  a  meaningful 
contribution  to  this  field  of  study. 


1.2  Literature  Survey 

Since  the  1920's,  considerable  attention  has  been  given  to  the 
derivation  of  probability  distributions  that  are  the  result  of  some 
algebraic  combination  of  random  variables  with  known  probability 
distributions.  Early  authors,  including  Aroian  (24),  Craig  (26),  and 
Craig  (27,28,29)  have  presented  detailed  discussions  on  sums, 
products,  and  quotients  of  independent  random  variables.  The  use  of 
Fourier  and  Laplace  transforms  extended  the  earlier  work  dealing  with 
sums  and  differences  of  random  variables.  Springer  provides  an 
excellent  discussion  and  bibliography  on  this  subject  in  his  book 


(17). 


The  problem  of  treating  products  and  quotients  of  random 


variables,  however,  was  limited  to  a  few  special  cases.  The  first 
practical  approach  for  dealing  with  products  and  quotients  of 
independent  variables  was  presented  by  Epstein  in  1948  (9).  Epstein 
used  the  Mellin  integral  transform  to  derive  the  probability  density 
functions  of  the  Student  t  and  Fisher  F  statistics  and  of  the  product 
of  two  standardized  normal  variates.  His  work  was  limited  to  two 
random  variables.  In  1966,  Springer  and  Thompson  (18)  extended  the 
work  of  Epstein  to  n  random  variables. 

The  most  significant  advances  in  the  algebra  of  random 
variables  came  in  1972  when  Carter  (3,4)  tied  together  the  physical 
science  work  on  H-functions  and  the  probability  work  on  Mellin 
integral  transforms  into  a  powerful  general  theory.  He  introduced  a 


/  * 


new  probability  distribution,  the  H-function  distribution,  which 
includes  as  special  cases,  ten  of  the  well  known  classical 
distributions  -  gamma,  exponential,  chi-square,  Weibull,  Rayleigh, 
Maxwell,  half-normal,  beta,  half-Cauchy,  and  general  hypergeometric. 
Carter  also  proved  that  products,  quotients,  and  rational  powers  of 
independent  H-function  variates  are  also  H-f unctions. 

In  1979,  Eldred  (7)  extended  the  work  of  Carter  by  developing 
a  computer  program  to  calculate  the  probability  density  function  of 
combinations  of  products,  quotients,  and  powers  of  H-function 
variates.  He  also  derived  the  H-functional  form  for  the  half-Student 
and  F  distributions. 

Cook  (5,6)  carried  this  work  even  further  to  provide  a  simpler 
method  for  calculating  the  H-function  resulting  from  some  combination 
of  H-function  variates.  He  also  developed  a  new  computer  program 
which  could  handle  sums  as  well  as  products,  quotients,  and  powers  of 
H-function  variates. 

Most  of  the  current  work  on  transform  and  H-function  theory 
for  the  algebra  of  independent  random  variables  is  given  in  books  by 
Springer  (17),  Mathai  and  Saxena  (14),  and  Giffin  (98)  and  in  papers 
by  Eldred  (7)  and  Cook  (5,6). 

Today,  the  H-function  techniques  are  powerful  enough  to  handle 
most  algebraic  combinations  of  independent  random  variables.  However, 
the  algebra  of  dependent  random  variables  has  received  little 
attention.  Much  of  this  has  been  due  to  the  inability  to  separate 


random  variables  residing  in  a  multivariate  density  function  due  to 
the  dependency  structure  of  the  density  function.  Indeed,  simply 
defining  a  bivariate  density  function  for  two  dependent  random 
variables  with  known  marginal  density  functions  has  been  a  major 
obstacle.  Except  for  the  bivariate  normal,  no  unique  bivariate 
density  function  can  be  derived  for  two  random  variables  with  a  given 
covariance  matrix  and  marginal  density  functions. 

In  the  1920's,  other  bivariate  distributions  were  constructed 
which  had  as  marginal  distributions  corresponding  well  known 
univariate  distributions  which  included  -  bivariate  Students  t, 
bivariate  beta,  and  Rhodes  distributions.  Little  more  was 
accomplished  until  the  development  of  a  bivariate  gamma  distribution 
in  1941. 

In  1960,  Gumbel  (37)  studied  a  bivariate  distribution  which 
has  exponential  margins,  but  no  meaningful  derivation  for  the 
distribution  is  known.  Marshall  and  Olkin  (48,49)  introduced  a 
bivariate  distribution  with  exponential  marginals  by  studying  a 
two-component  system  which  fails  to  function  after  a  shock  to  one  or 
both  components.  Further,  by  making  a  simple  variable  transformation 
they  were  able  to  express  a  meaningful  bivariate  Weibull  distribution 

In  a  recent  book,  Mardia  (47)  provides  an  excellent  summary 
and  bibliography  of  most  of  the  well  known  bivariate  distributions 
derived  up  through  1970.  Ord  (55)  provides  a  similar  summary  on 
families  of  frequency  distributions  which  include  certain  classes  of 
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bivariate  distributions.  A  summary  of  the  more  classical  bivariate 
distributions  given  by  these  two  authors  is  given  in  Appendix  A.  More 
recent  contributions  include  a  derivation  of  a  compound  gamma 
bivariate  by  Hutchinson  (42)  in  1981  and  a  new  class  of  bivariate 
logistic  distributions  by  Ali  and  Mikhail  (21)  in  1978.  Also,  current 
interests  have  extended  certain  known  bivariate  distributions  into  the 
complex  space.  Recent  articles  include  those  by  Brock  and  Krutchoff 
(25),  Giri  (36),  and  Saxena  (58). 

With  the  exception  of  the  bivariate  normal  distribution,  most 
of  the  work  in  multivariate  analysis  has  been  in  the  area  of 
characterizing  a  given  bivariate  distribution  either  through  a 
thorough  study  of  its  marginal  and  conditional  probability 
distributions,  or  through  a  study  of  its  characteristic  function. 

Such  analysis  is  reproduced  in  books  by  Springer  (17),  Anderson  (23), 
Feller  (34,35),  Mood  and  Graybill  (52),  and  Parzen  (56).  In  addition, 
recent  articles  on  the  subject  include  those  of  Lukacs  and  Beer  (46) 
and  Abrahams  and  Thomas  (20). 

Work  on  the  actual  distribution  of  algebraic  combinations  of 
dependent  random  variables  has  been  mostly  limited  to  the  bivariate 
normal  distribution.  In  his  book.  Springer  (17)  defines  the  bivariate 
normal  distribution  as  well  as  an  established  methodology  for  deriving 
the  sums  of  bivariate  normal  variates  using  double  Fourier  integral 
transform  techniques.  The  distribution  of  the  product  of  two 
dependent  normal  variates  was  studied  by  Aroian  (24)  in  1947,  and  the 


distribution  of  the  quotient  of  two  dependent  normal  variates  by  Craig 
(27)  in  1942.  Nicholson  (54)  used  a  geometrical  approach  to  study  the 
ratio  between  two  dependent  variables,  but  was  able  to  apply  his 
results  only  to  the  bivariate  normal. 

Along  with  new  interests  in  studying  and  developing  new  and 
meaningful  bivariate  distributions  came  revived  interests  in  studying 
products,  quotients,  and  sums  of  the  dependent  variates  of  these 
distributions.  Current  research  is  very  specific  in  nature  in  that 
techniques  used  for  one  bivariate  functional  do  not  apply  to  another 
bivariate  functional.  Indeed,  the  techniques  for  products  or 
quotients  may  not  even  be  the  same  for  a  given  bivariate  density. 
Current  research  in  this  area  include  Abrahams  and  Thomas  (20),  Alsina 
and  Bonet  (22),  Gupta  (38),  Lee,  Holland,  and  Flueck.  (45),  Mathai 
(50),  Tan  (60),  and  Wallgren  (61).  From  these  studies  it  becomes 
clear  that  a  more  general  theory  that  is  applicable  to  a  wide  variety 
of  cases  is  a  much  needed  tool  for  the  multivariate  analyst. 

In  1944,  Reed  (16)  defined  the  double  Mellin  integral 
transform  and  its  inversion  integral  with  associated  theorems  for 
each.  He  also  presented  a  theorem  for  deriving  the  double  Mellin 
transform  of  the  product  of  two  bivariate  functionals.  As  an  example 
of  his  theorems,  Reed  derived  the  double  Mellin  transform  identities 
for  Appell's  hyper geometric  functions  of  two  variables. 

Fox  (10)  provided  the  first  practical  method  of  handling 
products  of  bivariate  density  variates  by  extending  the  results 


derived  by  Reed  (16)  to  cases  of  statistical  distributions.  He  also 
included  some  work  on  quotients  of  bivariate  density  variates  as  well 
as  a  detailed  discussion  on  how  to  handle  functionals  that  resided  in 
other  than  the  first  quadrant.  Fox's  work  however,  was  limited  to 
finding  a  bivariate  density  function  that  resulted  from  products  or 
quotients  of  variates  of  two  bivariate  density  functions. 

Subrahamian  (19)  was  the  first  to  provide  significant  insight 
as  to  deriving  a  univariate  density  function  which  resulted  from  a 
product  or  quotient  of  two  dependent  variables  which  share  a  single 
bivariate  distribution.  He  used  the  results  of  Fox  (10)  and  combined 
them  with  earlier  results  on  independent  variates  by  Epstein  (9)  and 
Springer  and  Thompson  (18).  Subrahamian  applied  his  conclusions  to 
the  bivariate  normal  distribution  and  derived  by  this  alternative 
method  the  same  results  previously  presented  by  Aroian  (24)  and  Craig 
(27,28). 

As  a  result  of  the  labors  of  the  authors  listed  above,  today's 
statisticians  have  a  set  of  powerful  tools  in  transform  theory  for 
handling  certain  algebraic  combinations  of  independent  and  dependent 
variables.  However,  it  would  seem  a  natural  and  powerful  extension  to 
this  theory  if  one  could  apply  certain  results  of  Carter's  H-function 
analysis  to  bivariate  distributions  as  well.  A  brief  review  of  work 
on  H-functions  of  two  variables  follows. 

In  the  early  1970's,  Verma  (86,87),  Mittal  and  Gupta  (75),  and 
Goyal  (69)  extended  the  H-function  defined  by  Fox  to  a  generalized 


H-function  of  two  variables.  In  their  book,  Mathai  and  Saxena  (14) 
reproduce  a  formal  definition  of  the  H-function  of  two  variables  as 
well  as  some  important  properties  and  identities  for  the  H-function  of 
two  variables.  They  point  out  the  importance  of  this  function  arises 
from  the  fact  that  it  contains  as  special  cases  H-functions  of  one 
variable,  G-function  of  two  variables,  Whittaker  functions  of  two 
variables,  and  Appell's  functions  of  two  variables. 

The  majority  of  H-function  work  has  been  highly  theoretical 
and  generally  restricted  to  a  few  special  cases.  Most  of  the  articles 
on  the  subject  are  by  authors  from  India  and  are  published  in  foreign 
or  little  known  journals,  and  are  not  easily  accessible  to  the  U.S. 
researcher.  Almost  no  applications  are  given  in  the  literature  and 
the  few  given  are  for  physics  and  engineering.  Due  to  the  notation 
and  the  curse  of  dimensionality,  this  problem  is  particularly  true  for 
the  H-function  of  two  variables.  Thus  far,  the  limited  work  done  on 
bivariate  H-functions  has  been  in  the  area  of  extension  of  work  from 
G-functions  to  H-functions,  Argarwal  and  Singhal  (62),  and  on 
identifying  special  cases  of  the  bivariate  H-function,  Anandani  (64). 
Most  of  the  significant  work  done  on  the  bivariate  H-function  lies  in 
the  area  of  solving  differential  or  integral  equations.  Solutions  of 
dual  integral  equations  by  the  use  of  H-functions  can  be  seen  in  the 
works  of  Pathak  and  Prasad  (78)  and  Saxena  (82,83).  In  1972,  Mittal 
and  Gupta  (75)  used  a  generalized  function  of  two  variables  to  solve 
certain  classes  of  integral  equations. 


Along  with  a  more  general  and  subsequently  more  useful 
definition  of  the  bivariate  H-function,  Goyal  (69)  provides  some 
insights  on  the  applicability  of  taking  the  Laplace  transform  of  the 
bivariate  H-function.  Other  work  in  this  area  includes  Prasad  and 
Maurya  (79)  and  DeAnguio  and  Kalla  (66). 


1.3  Bivariate  Probability  Theory 
1.3.1  Definitions: 


(35:66-74; 56: 354-365; 52: 82-98, 198-215; 23:60-66;43:4-l2) 

Let  X  and  Y  be  jointly  distributed  random  variables  having  the 
bivariate  density  function  denoted  by  f^  Y(x,y)  and  the  bivariate 
cumulative  distribution  function  by  F^  Y(x,y).  Then 

FXY(x,y)  «  P(X  <  x  ,  Y  <  y)  (1.1) 

If  we  are  interested  only  in  the  cumulative  distribution  of  X, 
then  it  is  apparent  from  (1.1)  that 

P(X  <  x)  -  FXY(x,-)  (1-2) 

Therefore,  Fx(x)  ■  Fx  y(x*")  defines  the  cumulative 
distribution  function  for  X  and  its  associated  density  function  is 
given  by 


fx(x)  *  /  fXtY<x»y)dy  (1*3) 

Similarly,  the  density  function  for  y  is  given  by 

m 

fY(y)  * 1  fx,Y(x’y)dx  u-m 

The  density  functions  fx(x)  and  fY(y)  are  known  as  the  marginal 
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density  functions  for  the  joint  density  function  fY  v(x,y). 

X,  Y 


2 

The  expectation  and  variance  ox  of  x,  if  they  exist,  are 


given  by 


and 


“  E(x)  ■  /  /  xfx  Y(x,y)dxdy 

— <>o  — oo  * 


00  00 

°x2  “  Var<x>  *  I  /  (x-wx)2  fXY(x,y)dxdy 

Similar  identities  are  given  for  the  expectation  uy  and 


(1.5) 


(1.6) 


variance  a  1  of  Y. 

y 


The  expectation  u  or  first  product  moment  of  X  and  Y,  is 


given  by 


Jxy  -  E(xy)  »  /  /  xyfXY(x,y)dxdy 


(1.7) 


Equations  (1.5)  and  (1.6)  can  be  written  in  a  more  general 
form  for  higher  ordered  moments  of  fXY(x,y).  Let  <*nl  n2  be  the  nl,n2 
ordered  moment  for  fx  y(x,y),  then  anl  n2  is  given  by 


anl,n2  “  E(xnlyn2) 


/  /  xn^yn2f X^y(x,y)dxdy 


(1.8) 


Then,  o10  -  mx,  a0>1  -  Wy,  and  <*lfl  - 


xy' 


1.3.2  Properties  of  Moments:  (56:354-365) 

Let  g(x,y)  be  a  function  of  X  and  Y  where  X  and  Y  are  jointly 
continuous  with  a  joint  density  function  fx  y(x,y).  The  expected 
value  of  a  function  of  two  real  variables,  E[g(x,y)],  is  defined  as 

ao  ao 

E[g(x,y)]  -  /  /  g(x,y)fXY(x,y)dxdy  (1.12) 

«* oo  ao 

From  this  definition,  the  following  linearity  property  for 
expectations  of  jointly  distributed  random  variables  is  derived. 
Theorem  1.1:  If  X  and  Y  are  real  random  variables  which  are  jointly 
distributed  by  fx  y(x,y),  and  if  X  and  Y  have  finite  expectations  E(x) 
and  E(y),  then  the  sum  X+Y  has  a  finite  expectation  given  as 

E(x  +  y)  -  E(x)  +  E(y)  (1.13) 

A  similar  relation  may  be  found  for  finding  the  variance  of 
two  jointly  distributed  random  variables. 

Theorem  1.2:  If  X  and  Y  are  real  random  variables  which  are  jointly 
distributed  by  fx  y(*»y)»  and  if  X  and  Y  have  finite  variances  Var(x) 
and  Var(y)  and  a  finite  covariance  Cov(x,y),  then  the  sum  X+Y  has  a 
finite  variance  given  as 

Var(x+y)  ■  Var(x)  +  Var(y)  +  2Cov(x,y)  (1.14) 

1.3.3  Moment  Generating  Functions:  (56:354-365;52:200-204) 

The  joint  moment-generating  function  for  a  probability  density 
function  of  two  variables  fx  Y(x,y)  is  given  by 


Cov(x.y) 


a 


FX-E(x),Y-E(y)(0,0) 

Proceeding  in  a  similar  manner  will  yield  any  degree  of 
moments  of  higher  order. 

Properties :  If  Z  *  X  +  Y,  where  X  and  Y  are  independent,  F^(t)  is  the 
moment  generating  function  of  X,  and  ?Y(t)  is  the  moment  generating 
function  of  Y,  then  the  moment  generating  function  of  #z(t)  is  given 
as 


Fz(t)  «  Fx(t)FY(t). 

Characteristic  Function;  In  cases  where  the  integral  in  (1.15)  does 
not  exist,  moments  may  still  be  found  in  an  identical  fashion  by  use 
of  the  characteristic  function.  The  characteristic  function  is 
identical  to  equation  (1.15)  with  t ^  and  t2  replaced  by  itj^  and  it2 
respectively. 


1.4  Transformations  of  Random  Variables 

Emphasis  in  this  area  is  on  the  use  of  integral  transforms  to 
obtain  probability  density  functions  for  certain  transformations  of 
random  variables.  First,  a  review  of  some  related  probability 
concepts  should  be  made. 

A  one-to-one  transformation  h(x)  from  a  set  S  into  a  set  T 
means  that  for  each  y,  an  element  of  T,  there  exists  one  and  only  one 
x,  an  element  of  S,  such  that  h(x)  ■  y.  When  a  function  h(x)  is  a 
one-to-one  transformation  from  a  set  S  to  a  set  T,  then  the  inverse 
transformation  h~*(y),  from  T  onto  S,  exists  and  h-1[h(x)]  -  x. 

Stating  that  a  set  S  in  the  set  of  positivity  for  a  transformation 
h(x)  means  that  S  is  the  set  of  values,  x,  for  which  h(x)  is  positive. 

Two  random  variables  X  and  Y  are  independent  if  their  joint 
probability  density  function  f^  y(x,y)  equals  the  product  of  their 
marginal  density  functions  f^(x)  and  fy(y).  This  means  that  any 
variation  in  X  will  in  no  way  affect  the  outcome  of  Y,  or  vice  versa. 
Theorem  1.3:  Let  X  be  a  random  variable  with  continuous  probability 
function  f^(x)  ant*  suppose  that  y  ■  h(x)  is  a  one-to-one 
transformation  from  S,  the  set  of  positivity  of  fy(x),  onto  T,  the 
image  of  S  under  h(x) .  If  h“*(y)  is  differentiable  and  its  derivative 
is  continuous  on  T,  then  the  probability  density  function  of  Y  may  be 
given  as 


fY(y) 


f(h_1(y)]  2.  h_i(y) 
dy 


,  y  e  T 


.else 


Theorem  1.4;  Let  X  *  (X^t  X2,  ...  Xk)  be  a  set  of  k  random  variables 
having  the  joint  continuous  probability  density  function  fx(x).  Let 
T  ■  h(x)  *{h^()c),  t^Cx).  **•  hk(x>>  b®  a  set  °f  relations  forming  a 
one-to-one  transformation  from  S,  the  k-dimensional  set  of  positivity 
of  fx,  onto  T,  the  k-dimensional  image  of  S  under  h(ic).  The  inverse 
transformation  exists,  X"1  -  h'^y)  -  (g^),  g2(y),  •••  gk(z.)}.  If 
the  partial  derivatives  of  h-*(£)  exist  and  are  continuous. 


8i j  -  *yi  {gi(yi.  y2»  •••  yk» 


i,j  *  1,2,  ••  .k 


then  the  joint  probability  density  function  of  Y  is  given  by 


*2^*  8k(Z>>|j| 

fy(Z)  "j 


,y  e  T 


.else 


where  J  is  the  Jacobian,  the  determinant  of  first  partial  derivatives 


Using  Theorem  1.4,  one  can  find  Che  distributions  for  the  sum, 
product,  difference,  and  quotient  of  two  jointly  distributed  random 
variables . 

Example  1.1:  Suppose  the  probability  density  function  of  Z  ■  X/Y  Is 
desired.  Let  W  »  Y  so  that  X  *  ZW  and  Y  *  W  and 

3x/3z  3x/3w  w  z 

J  ■  ■  ■  w 

3y/3z  3y/3w  0  l 

By  Theorem  1.4,  for  the  approprate  ranges  of  z  and  w, 

fZ,W^z,w^  "  fX,Y*zw,w^  lJ  I 

The  marginal  distribution  of  Z  *  X/Y  is  found  by  integrating  the  above 
joint  distribution  f^  y(z,w)  over  the  proper  range  of  w: 

00 

fz(z)  -  /  fXY(zy,y)ydy 

Notice  that  if  W  ■  X  so  that  Y  ■  W/Z  then  the  determinant  of  the 


Jacobian  equals  -w/z  .  Then 


fz(z)  ■  /  (x/z^)fx  y(x,x/z)dx 

The  point  here  is  that  the  two  forms  for  calculating  fz(z)  are 
not  identical  by  a  simple  interchange  of  the  x  and  y  variables.  This 
is  the  only  case  where  the  symmetry  does  not  hold.  While  both  forms 
are  certainly  valid,  the  first  form  is  the  one  most  commonly  seen  in 
the  literature  and  is  usually  the  most  easily  applied. 

Using  Theorem  1.4  similarly  to  find  the  distributions  for  the 
difference,  product,  and  sums  of  two  random  variables  gives  the 
following  theorem  (41). 

Theorem  1.5:  If  X  and  Y  are  jointly  continuous  random  variables  with 
probability  density  function  f x  y(x,y),  then 

(1)  the  probability  density  function  of  the  random  variable 
Z  ■  X  +  Y  is  given  by 

00  00 

fZ(z)  *  /  fX,Y<Z-y*y)dy  ■  ^  fx>y(x,z-x)dx 

If  X  and  Y  are  independent,  then 

00  00 

fz(z)  -  ;  fx(z-y)fy(y)dy  ”  1  fx(x)fY(z-x)dx 
—00  —00 

(2)  the  probability  density  function  of  the  random  variable 


Z  *  X  -  Y  is  given  by 


f Z(Z)  ■  /  fx  Y(z+y,y)dy  »  /  fx  Y(x,z-»-x)dx 


If  X  and  Y  are  independent,  then 


fz(z)  *  /  fx(z+y)fY(y)dy  ■  /  fx(x)fy(z4x)dx 

(3)  the  probability  density  function  of  the  random  variable 
Z  ■  XY  is  given  by 


fz(z)  -  /  |y~1|fX  Y(z/y.y)dy  ”  /  |x_1jfx  Y(x,z/x)dx 

If  X  and  Y  are  independent,  then 

00  00 

fz(z)  -  /  |  y  1 |fx(z/y)fY(y)dy  «  /  | x”1 | fx(x)fY(z/x)dx 

00  111100 


(4)  the  probability  density  function  of  the  random  variable 
Z  -  X/Y  is  given  by 

00  00 

fz(z)  *  |  y  |^X,Y^zy •yWy  ■  ^  |  ~2  I^X.Y^’z-^* 

If  X  and  Y  are  independent,  then 


fz(«)  -  /  |y  |fx(zy)fY(y)dy  -  /  |  i2j  fx(x)fY(|.)dx. 

'■oo  2 


Oy  ■  1  and 


This  is  equivalent  to  Craig's  (28)  derivation  if  » 

U  *  Li  *  0. 

x  y 

Example  1.3:  Now  consider  Morgensterm' s  bivariate  uniform  density 
given  by 

fx>Y(x,y)  -  1  +  p(2x  -  l)(2y  -  1)  0  <  x,y  <  1 

Let  Z  *  XY.  Then 

00 

fz(z)  -  /  | y-1| fXfY(z/y»y)dy 

The  lower  limit  of  integration  of  y  is  determined  from  the  relation 
x  -  z/y,  0  <x  <  1,  0  £  y  £  1 

Since  x  £  1,  the  lower  limit  on  y  is  z.  Furthermore,  since  x  £  1  and 
y  £  1,  the  maximum  value  of  z  *  xy  is  1.  Then 

1 

fz(z)  -  ;  (l/y) [1  +  P(2z/y  -  l)(2y  -  l)dy 

z 

-  [p(4z  +  1)  +  1 ]ln(l/z)  +  4p(z  -  1)  ,  0  <  z  £  1 

If  p  »  0,  then  X  and  Y  are  independent  and  fz(z)  reduces  to 
fz(z)  -  -ln(z)  ,  0  <  z  £  1 

This  is  the  same  form  as  that  derived  by  Springer  (17:91-94)  for  the 
case  of  the  product  of  two  independent  identically  distributed  uniform 


variates . 


Example  1.4:  Suppose  the  probability  density  function  of  Z  ■  XY  is 
desired  where  X  and  Y  are  jointly  distributed  by  Kellogg-Barnes  III 
distribution  given  as 


fX,Y^x,y) 


ga2  vc  -(ax  +  Bxy) 
r(c) 


x,y>0 

a,S>0  ,  c>2 


from  Theorem  1.5,  case  (3),  the  distribution  for  Z  is  given  by 


fz(z)  *  /  -j-  fX  Y(x.z/x)dx 


.  &2L-  e_8z  /  xc_1 

r(c)  Jo  x 


-<» 


dx 


*  ge“®z  z>0 

which  is  the  univariate  exponential  distribution  with  parameter  1/g. 
Example  1.5:  Consider  the  Kellogg-Barnes  II  distribution  given  by 

f„  v<x,y>  -  BaZ  .-<«  +  Sy/‘>  x,y>0 

a,  B>0 


From  Theorem  1.5,  case  (4),  the  distribution  of  the  random  variable 
Z  «  X/Y  may  be  found  by 


f-(z)  -  f  yf„  Y(zy,y)dy 

-  3a2  e_6/z  /"  ye-012^  dy 

0 

-  (8/z2)e'8/z 

Now  suppose  Che  reverse  is  desired,  Chat  is,  Che  discribuClon 
of  Z  *  Y/X.  Theorem  1.5,  case  (4),  scill  applies  by  using  a  simple 
change  of  variables  and  fz(z)  is  given  by 

OO 

fz(z)  ■  /  xfx  y(x>zx)^x 

■  8a2  e-^2  /  xe”ax  dx 
0 

-  8e“8z 

While  Examples  1.4  and  1.5  are  relaCively  scraighc  forward. 
Examples  1.2  and  1.3  show  chac  Che  use  of  Theorem  1.5  for  produces  and 
quodencs  of  dependenc  variables  can  be  an  arduous  Cask..  For  products 
and  quotients  of  dependenc  variables  Che  task  can  be  simplified 
through  Che  use  of  Mellln  transform  techniques  as  will  be  shown  in  the 


next  chapter. 


CHAPTER  2 


Application  of  Integral  Transforms  to  Statistical  Analysis 
2.1  General  Remarks 

Section  1.4  of  Chapter  1  shows  that  using  convolution 
integrals  for  the  transformation  of  variables  to  find  the  distribution 
for  the  sum,  difference,  product,  or  quotient  of  two  random  variates 
can  be  a  difficult  task.  This  chapter  outlines  techniques  to  simplify 
the  problem  by  utilizing  integral  transform  techniques.  A  review  of 
integral  transforms  and  the  associated  techniques  for  finding 
distributions  of  algebraic  combinations  of  independent  variates  is 
presented  followed  by  a  discussion  of  the  extension  to  products  and 
quotients  of  dependent  variates  using  double  Mellin  transform 
techniques. 

Since  this  dissertation  is  devoted  to  products  and  quotients 
of  dependent  variates  using  double  Mellin  transform  techniques,  the 
double  Mellin  transform  is  developed  more  completely.  Theorems 
governing  its  use  as  given  by  Fox  (10)  and  Reed  (16)  are  presented. 
Extensions  to  univariate  Mellin  transform  properties  are  presented  for 
the  double  Mellin  transform.  These  properties  prove  useful  to  Mellin 
transform  manipulations  in  later  chapters. 

Finally,  theorems  are  given  for  the  distribution  of  products, 
quotients,  and  rational  powers  of  two  dependent  variates.  Theorems 
are  also  presented  for  products,  quotients,  and  powers  of  variates 
from  two  bivariate  distributions  which  are  pairwise  independent. 


2.2  Integral  Transforms 


2.2.1  Double  Fourier  Transform:  (17:67-75;106:76-79) 

A  real  function  of  two  variables  fx  Y(x,y),  where  each 
variable  is  defined  over  the  whole  real  line,  is  doubly  Fourier 
transformable  if  the  integral 

-  -  iklX  +  ik2y 

f  ■  |fx  Y(x,y)|e  dxdy 

converges  for  some  real  value  for  and  k2»  Then, 


ci  »c 


1  ’  2 


{fX,Y(x*y)}  m  f  * 


4tlx 


it2y 


fX,Y(x‘y) 


dxdy 


(2.1 


is  the  double  Fourier  transform  of  fY  v(x,y).  F..  (fY  v(x,y)}  is 

a,9  i  * 

called  the  bivariate  characteristic  function  of  fY  v(x,y),  and 

A,  I 

e*Clx  +  it2y  ca^2.ed  the  kernel.  The  inversion  integral  is  given  by 


h+i-  k+i«  -it.x  -  it2y 

fX,Y(x,y)  *  2v  *  f  e  F(t1,t2)dt1dt2  (2.2) 

h-i«  x-i» 


2.2.2  Double  Laplace  Transform:  (l06:221-228;34:452-458) 

A  real  function  of  two  variables  f^  Y(x,y),  defined  everywhere 
for  x  >  0  and  y  >  0,  with  x  and  y  real,  is  doubly  Laplace 
transformable  if  the  integral 


y 


-  -  -kxx  -  k,y  .  , 

/  /  e  I  fx^Y(x,y)|dxdy 


converges  for  some  real  values  of  k^  and  k2»  Then 

•  •  -rjx  -  r2y 

Lr  -  tfX  Y<x*y>>  m  *  f  e  fX  v(x.y)dxdy 

rl  »r2  *•*  0  0  ’ 

is  the  double  Laplace  transform  of  fx  Y(x,y),  where  rt  and  r2  are 
complex  variables. 

From  this  definition  it  can  be  shown  that 


(2.3) 


^1  r^(fX  Y(ax,by)}  -  (ab)  L^/a.rj/b) 


(2.4) 


The  double  Laplace  transform  has  been  used  in  the  past  for 
functions  of  x  and  y  where  x  and  y  are  independent.  Only  limited 
consideration  has  been  given  to  dependent  functions  of  x  and  y. 

2.2.3  Double  Mellin  Transform:  (17 : 151-156; 10; 16; 19) 

A  real  function  of  two  variables  fx  Y(x,y),  defined  everywhere 
for  x  >  0  and  y  >  0,  with  x  and  y  real,  has  been  defined  by  Reed  (13) 
to  have  the  double  Mellin  integral  transform  given  as 


•  •  s.-l  8,-1 

m(8i»82)  ■  /  /  x  y  fx  Y(x,y)dxdy 
0  0  ’ 


and  its  inverse  as 


(2.5) 


V.V.V.V - 


fX,Y^x,y^ 


1 


(2.6 


h+i«  k+i»  -s  ^  -S2 

f  f  x  y  M(s. ,s7)ds1ds. 

(2ui)‘  h-i-  k-i- 

The  conditions  for  which  (2.5)  and  (2.6)  are  valid  are  stated 

by  the  following  theorems.  The  proofs  are  given  by  Reed  (16)  and  Fox 

(10). 

Theorem  2.1 :  If 

(i)  M(s^,s2)  is  a  regular  function  of  both  variables  s^,s2  in 

the  strips  a  <  Sj  <  b,  c  <  s2  <  d 

(ii)  in  these  strips  M(slfs2)  -  0(  jst j “m)0( Js2 j~n)  for  some 

m  >  0,  n  >  0,  as  |SjJ  and  |s2|  tend  to  infinity  independently  of  each 

other ; 

(iii)  a  <  h  <  b  and  c  <  k  <  d 

00  00 

(iv)  f  f  j  M(Sp82)  |  |dsjJ  |ds2|  exists  when  taken  along 

—00  —OO 

any  lines  parallel  to  the  imaginary  axis  in  the  strips  defined  in  (i) 

(v)  f^  y(x*Y)  defined  by  equation  (2.6) 

then 


-  -  sri  s2-l 

^81*s2>  "  X  y  fx,Y(**y)dxdy 


is  true. 


Theorem  2.2:  Let  X  denote  a  part  of  the  complex  x  plane  which  is 
bounded  by  two  lines  through  the  origin  and  which  includes  the  whole 
of  the  positive  real  axis  from  0  to  +*.  Let  Y  denote  a  similar  region 
in  the  complex  y  plane.  If  with  x  in  X  and  y  in  Y  the  following 
conditions  are  satisfied: 

(i)  there  exists  two  real  numbers  h  and  k  such  that 
x  yk^X  y(*,y)  is  a  regular  function  of  both  x  and  y; 

(ii)  *hykfXtY<x,y)  “  0( ) logx | ”m)0( |l°8y I ~n) ♦  m  >  0,  n  >  0,  as  x 
and  y  tend  to  infinity  independently 

(ill)  //jx11  yk  fx  y(x,y) j |dx | jdy|  exists,  when  taken  along  any 
lines  in  the  X  and  Y  regions 

(iv)  M(Sj ^2)  is  defined  by  equation  (2.5) 
then  equation  (2.6)  is  true 

These  two  theorems  give  an  exact  analogue  to  the  single  Mellin 
Integral  transform  theorems. 

2.2.4  Mellin  Transform  Properties 

By  making  the  change  of  variables  x*u/a  and  y»v/b  in  the 
defining  integral  (2.5), 

m  co  g  g 

MS  .  {  f(ax,by)  }  -  /  /  f(u,v)u  1  v  ds.ds, 
l’s2  00 

-  a  Sl  b  82  M  {  f(u,v)  }  (2.7) 

91*s2 

Since  multiplying  f(x,y)  by  xay^  merely  results  in  changing  s^ 


to  Sj+a  and  S2  to  S2+b,  we  also  have  the  following  relation: 


M  {  xayb  f(x,y)  >  *  M  +a  ..  {  f(x,y)  } 

3^,82  ®  1  '  b 


(2 


For  a,b  >  0,  making  the  change  of  variables  x*u^a  and  y«v^b  yields 


Ms  «  «<*' 
*l»s2 


.  •  •  s./a-l/a  s,/b-l/b  .  . 

a,yb)}  ■  /  /  f(u,v)u  1  v  2  ((ab)*1u1/a_1v1/b  *dudi 


0  0 


.  00  00  Sl/a-l  s2/b-l 

(ab)_1  /  /  f(u,v)u  v  dudv 


0  0 


<ab>'1  ’V.vb1  f(u>v>  1 


(2.' 


The  Mellin  transform  also  has  certain  unique  properties  for 
derivatives  and  integrals  of  functions.  By  definition 


Ms  <  fx  >  *  /  /  xSl  1  y82  1  [  |  f(x,y)  ]  dxdy 

S1  »s2  3X  0  0  3x 


If  [ (3/3x)f (x,y) ]  is  continuous  for  y  constant  and  0  <  y  <  *  then  the 
equation  above  can  be  written  (89; 92 : 179-180) 


s2“l  "  Si-1 


Ms  s  (  fx  f(x.y)  }  -  /  y  4  /  x  1  j  f(x,y)  dxdy 

sl»®2  00 


integrating  the  inner  integrand  by  parts  yields 


f(x,y)  -  0 


37 


then  the  following  hold. 


Ms  a  t  Tx  f<x»y)  J  *  “(s2_1)^s  „  J  f(x,y)  > 

s1,Si  <*x  ‘  Sl»s2  1 


(2.12) 


M  {  /  f(x,v)dv  }  -  s2"  M  {  f(x,y)  } 

*l*s2  y  *1»s2T1 


(2.13) 


Higher  derivatives  can  be  dealt  with  in  a  similar  fashion. 
Applying  (2.11)  and  (2.13)  iteratively,  the  following  property  for  a 
double  integral  can  be  written. 


M  {  /  /  f(u,v)dudv  )  -  (s.s,)"1  M  {  f(x,y)  } 


a  -  i  .  ' 

sl,s2  x  y 


12'  Sj+l^+l 


(2.14) 


Equation  (2.14)  above  is  a  useful  property  in  that  it  may  now 
be  used  in  developing  the  cumulative  of  the  H-function  distribution. 
The  properties  developed  above  are  summarized  below. 


M«  ,  <  )  -  a  1  b  2  H  (  ) 


*V  t(X'y)  1  ‘  M.l+.,,2+b<  «*•»>  1 

f<*a-yb)  1  *  <abrl  M.l/a,.2/b(  f(,-y>  >  ;  *•’>  2  0 


"sj.Sjl  'j  )  -  («1»2)_1  M,1+it,2+i(  t 


2.2.5  Mellin  Transform  of  Appell's  Functions:  ( 16 ; 8 : 232 ) 

Reed  (16)  used  the  double  Mellin  Integral  transform  to  obtain 


Che  following  transforms  for  Appells  hypergeometric  functions  of  two 
variables : 


00  -  s1-i  s 2 “l 

fix  y  F. (a ,b ,b' ;c ; -x,-y)dxdy 
0  0 


r(c)r(s1)r(s2)r(a-s1-s2)r(b-s1)r(b,-s2) 


r(a)r(b)r(b’)r(c-8l-s2) 


when  0  <  Re(s1+s2)  <  Re(a),  0  <  Re(Sl)  <  Re(b),  0  <  Re(s2) 


•  *  sri  s2-l 

fix  y  F_(a,b,b' ;c ,c ' ; -x,-y)dxdy 
0  0 


r(s1)r(s2)r(c)r(c,)r(a-s1-s2)r(b-sl)r(b,-s2) 


r(a)r<b)r(b'  ^(c-sprcc'-sj) 


when  0  <  Re(Sj+s2)  <  Re(a),  0  <  Re(s^)  <  Re(b),  0  <  Re(s2) 


•  •  Sj^-1  s  2~1 

fix  y  F3(a,a',b,b';c;-x,-y)dx  dy 
0  0 


r(s1)r(s2)r(c)r(a-s1)r(a,-s2)r(b-s1)r(b,-s2) 

r(a)r(a’)r(b)r(b,)r(c-s1-s2) 


(2.1 

<  Re(b') 

(2.1 

<  Re(b’) 

(2.1 


when  0  <  Re(8l)  <  min(Re(a),  Re(b)),  0  <  Re(s2)  <  min(Re(a' ) ,Re(b' )) 


2.3  Integral  Transforms  for  Independent  Variates 

From  the  examples  It  can  be  seen  that  using  Theorem  1.5  can  be 
a  difficult  process.  Integral  transforms  can  help  simplify  the 
process.  For  the  case  of  X  and  Y  independent,  the  following  formulas 
have  been  most  helpful  in  solving  algebraic  combinations  of  random 
variables. 

Let  Ft,  Lp,  and  Mg  represent  the  Fourier,  Laplace,  and  Mellin 
integral  transforms  of  one  variable  respectively.  Then  special 
properties  of  these  transforms  are 


Ft{fx(x)}  Ft(gY(y)}  •  Ft{  /  fx(x)gY(y-x)dx  } 


Lrtfx(x)}  Lr(gY(y)}  -  Lr {  /  fx(x)gY(y-x)dx  } 

00 

Mg{fx(x)}  Mg{gY(y)}  -  Mg{  /  x”1  fx(x)gy(y/x)dx  } 

Combining  these  formulas  with  Theorem  1.5,  the  following  results  can 
be  derived: 

(1)  the  probability  density  function  of  the  random  variable 
Z  -  X+Y  is  given  by 

fz(z)  -  Pl“l(  Ft(fx(x)]  Ft[fY(y)]  } 
or 

fz(z)  -  Lr[fx(x)J  Lr[fy(y)]  }  x.y.z  >0 

where  *  and  are  the  inverse  Fourier  and  Laplace  transforms  of 

one  variable  respectively. 

(2)  the  probability  density  function  of  the  random  variable 
Z  *  X-Y  is  given  by 

fz(z)  -  Fj"l{  Ft[fx(x)]  Ft[fY(-y) ]  } 

where  F^”*  is  the  inverse  Fourier  transform  of  one  variable. 

(3)  the  probability  density  function  of  the  random  variable 
Z  ■  XY  is  given  by 


fz(z)  -  M#[fx(x)J  Mg[fY(y>]  }  x.y.z  >  0 


is  the  inverse  Mellin  transform  of  one  variable 


■>  myw 7 


where  M^ 

(4)  the  probability  density  function  of  the  random  variable 
Z  ■  X/Y  is  given  by 

fz(z)  «  Mg[fx(x)]  M2_g[fY(y)J  }  x,y,z  >  0 

where  M^-1  is  the  inverse  Mellin  transform  of  one  variable. 

A  distinct  advantage  to  transform  techniques  is  that  these 
above  formulas  can  be  easily  extended  to  more  than  two  variables. 
However,  these  formulas  are  restricted  to  cases  where  X  and  Y  are 
independent . 

2.4  Mellin  Transforms  for  Dependent  Variates 

In  1957,  Fox  (10)  applied  the  double  Mellin  integral  transform 
to  the  theory  of  bivariate  statistics  and  derived  the  following 
conclusions  for  relations  of  two  bivariate  distributions. 

The  expectation  of  $(x,y),  E($(x,y)],  is  defined  by 
(31:260-265) 


E[4>(x,y)]  -  /  /  $(x,y)f(x,y)  dxdy 

0  0 


If  <p(x,y)  -  x 


Sj^-l  s2“l 


,  then  E($(x,y))  is  the  double  Mellin  transform 


definition  given  by  (2.5).  Let  fx  ^  (x^y^  and  fx  Y  (x2,y2)  be  two 
bivariate  density  functions  having  the  double  Mellin  transforms 
^(8^,82)  and  M2(s1,s2)  respectively,  and  where  x^,  x2,  y^,  y2  >  0. 
Further,  assume  X^,Y^  are  Independent  of  X2,Y2«  Then 


El^x^y^^Xj.yj)]  *  E(4'1(x1,y1)]  E[4>2(x2,y2)] 


since  X^y^  is  pairwise  independent  of  X2,Y2*  If  4li(x^,y1)  ■ 
srl  s2-l  s^-1  s2“l 

xt  yx  and  <j>(x2,y2)  -  x2  y2  .  then 

sl~l  s2~l 

E[  (xlX2)  (Yj^)  1  *  Ml(8l»82)  M2(81»s2) 

The  joint  probability  density  function  of  (Xjx2,yiy2)  can  now  be 

Sj-1  s2-l 

solved  by  substituting  E[(x^x2)  ^1X2^  1  i^to  the  right  hand 

double  integral  of  (2.6): 

^Z,w(z,w)  *  M2  ^  I  M^(s^,S2)M2(Si ts2)  ]  (2.19) 

where  Z  ■  XjX2,  W  -  Y^Y2,  ant^  M2~*  *s  the  Mellin  inversion 

transform  as  defined  by  (2.6).  Continuing  in  this  fashion.  Fox  showed 
the  following  to  be  true. 

fZ,W<z’w)  "  m2_1  l  W1(2-s1,s2)  ]  (2.20) 

for  Z  -  1/Xt  and  W  -  Yx 

fZ,W<z’w)  "  M2_1  1  Ml(2-sl,2-s2)  ]  (2.21) 


for  Z  -  1/X,  and  W  ■  1/Y, 


Subrahamian  (19)  combined  the  work  of  Fox  (10)  to  derive 
results  for  cases  of  products  and  quotients  of  two  dependent 
variables.  Let  X  and  Y  be  dependent  random  variables  with  probability 
density  function  f^  Y(x*y)  which  is  positive  in  the  first  quadrant  and 
zero  elsewhere.  Further,  suppose  that  the  double  Mellin  transform  of 
fXY(x,y)  exists  and  is  given  by  M(s1,S2).  If  <Kx,y)  ■  xs-1ys-1  then 

OP  OP 

E[  xs-1ys_1  ]  -  /  /  xs_lyS_1  f„  v(x,y)  dxdy 
0  0  A,T 

*  Ms,s{  fX, Y(x’y>  } 

Let  Z  *  XY,  then 

E[  xs_1ys_1  ]  -  E[  zs_l  J 

Substituting  E[  z3-*  ]  into  the  inversion  integral  for  the  univariate 
Mellin  transform,  the  density  function  for  Z  can  be  found. 
Specifically, 

fz(z)  -  Mt“l[  M(s,s)  }  (2.23) 

where  Z  ■  XY  and  M^-*  is  the  univariate  Mellin  transform  inversion 


operator 


Similarly,  Che  probability  density  function  of  the  random 


variable  Z  ■  X/Y  is  given  by 

f2(z)  -  M(s,2-s)  1 


(2.24) 


The  work  accomplished  by  Fox  and  Subrahamian  can  be  extended 
to  n  pairs  of  dependent  variables  which  are  pairwise  independent. 

This  work,  shown  in  the  following  theorems,  is  similar  to  that  done  by 
Carter  (3)  for  independent  variables.  Theorems  extending 
Subrahamian' s  work  to  raise  dependent  variates  to  rational  powers  are 
also  given. 

Theorem  2.3:  If  X  and  Y  are  dependent  continuous  random  variables 
with  a  bivariate  probability  density  function  fx  Y(x,y),  x»Y  >  0,  then 
the  bivariate  probability  density  function  of  Z  ■  Xa  and  W  *  Y^,  a,b 
rational,  is  given  by 


M_-1  M 


2  as  -a+l,bs,-b+ll  X,Y 


{  fY  Y(x,y)  }  z,w  >  0 


fZ,W(z*w) 


(2.25) 


otherwise 


Proof  of  Theorem  2.3 

The  Mellin  transform  of  the  probability  density  function 
^Z,w(z*w)  *s  8i-ven  by 


] 


s 

-  E(  z 


I"1  s2-1 
w 


where  E  is  the  expected  value  operator.  From  the  definition  of  Z  and 
W,  this  bacoaes 


sl~l  S2~l 

Mg  {  fz  W(z,w)  }  -  E[  (xa)  (yb)  ] 

Sl»s2  6,tl 

*  •  as^-a  bs2~b 

"  Vo  fX,Y<x»y>x  y  dslds2 

u  'J 

•  •  (asj-a+O-l  (bs2-b+l)-l 

*  Vo  fx»Y<x,y)x  7  dSld82 

"  Mas1-a+l,bs2-b+lt  fX,Y(x,y)  } 

Froa  the  inverse  Mellin  transform,  (2.25)  follows. 

Theorem  2 .4:  If  X  and  Y  are  dependent  continuous  random  variables 
with  a  bivariate  probability  density  function  fx  Y(x,y),  x,y  >  0,  then 
the  probability  density  function  of  the  random  variable  Z  *  XaY**,  a,b 
rational,  is  given  by 


fz(z) 


M1  Mas-a+l,bs-b+l{  fX,Y^x*y^  *  z  >  0 

0  otherwise 


(2.26) 


Proof  of  Theorem  2.4 


Let  U  -  Xa  and  V  »  Yb.  From  Theorem  2.3 

fu>vU,v)  -  M2  1  [  MaSi_a+1>bs2_b+1 
from  which  it  follows  that 


t  fx,Y(x,y)  *  1 


\,52(  fu,v(u''')  >  -  Masl-a+l,bs2-l>+l1  fx,Y<x'>'>  1 


Since  Z  -  UV,  applying  Equation  (2.24)  yields 


fz(z)  -  Mj^”1  [  Mg  g{  fu  v(u,v)  }  ] 

"  M1  Mas-a+l,bs-b+l{  fX,Y(x»y)  }  1 

Example  2.1:  (19; 17 : 154 -156) 

Consider  the  bivariate  standard  normal  distribution.  For 
x,y  0,  the  double  Mellin  transform  is 

i  "  *  ®i'l  S2~l  r  o  o  "] 

Mv  Y(s1,s2)  -  — ;==s-  f  f  X  y  exp  I  - i—  (x  -pxy-yz)  i  dxdy 

1  L  211/l-p2  0  0  1 2(l-p2)  J 


,  [(s1+s2+1)/2-1) 

r(s1)r(s2)(i-pz) _ 


(si+s,+l)/2 

2  r(i/2)r{(s1+s2+i)/2] 


3  2  3i~f82+1 
2  ’  2 


From  Theorem  2.4,  for  Z  ■  X/Y,  a  -  1  and  b  ■  -l.  The  Mellin  transform 
of  f  (*)  is  then  given  by 


,  r(3)r(2-3) 

2T(l/2)r(3/2) 


i 


.2  z2  .  J  • 
*  2  ’  2  * 


1 


The  transform  above  is  valid  for  z  >  0  only,  but  by  symmetry, 
the  inverse  for  z  <  0  may  be  derived  also.  Subrahamian  completed  the 
inversion  to  give 


MO 


A  -  p2 


ir(l  -  2pz  +  z  ) 


This  is  exactly  the  form  derived  in  Example  1.2  for  a  -  o  «  1,  and 

x  y 

P  a  LI  a  0. 

x  My 

Example  2.2;  Consider  again,  Morgenstera* s  bivariate  uniform  density. 

fX,Y<x»y>  “  1  +  P(2X  -  D(2y  >1)  0  <  x,y  <  l 

The  double  Mellin  transform  of  this  density  is  given  by 

My  y<s1’s2>  "  /"/"  xSl  1  y82  1  [l  +  P(2X  '  1)(2y  "  l^d*dy 
*  0  0 

1  P(«!  -  1)(92  ~  1) 

sl82  sl92^sl  +  1)(®2  + 


If  one  is  interested  in  the  density  function  fz(z),  where  Z  »  XY,  then 
from  Theorem  2.4  a  ■  b  ■  1.  The  Mellin  transform  of  fz(z)  is  then 


sz(s  +  1)Z 


s 


L 


Ths  density  function  of  fz(z)  *8  found  by  performing  the  inverse 
operation  of  Mz(s). 


fz(z) 


1_ 

2U 


h+i®  r  -s 
h-i®  L  sz 


Prs  (s  -  d2|Hc 

s2(s  +  l)2  J 


The  density  of  fz(z)  may  be  found  by  performing  the  inversion  integral 
directly,  or  by  summing  the  residues  of  two  terms  at  s  *  0  for  term 
one  and  s  -  0,  s  *  -1  for  term  two. 

For  term  one,  *  residue  at  s  -  0. 


R1  "  ds  Z 


-s 


ln(l/z)  0  <  z  <  1 


For  term  two,  ■  residue  at  s  ■  0. 


.  L  Pz'3(3  -  D2 


ds 


(s  + 1  y 


■  -4p  +  pln(l/z)  0  <  z  <  1 


For  term  two,  ■  residue  at  s  ■  -1. 

“8  2 

R  d  pz — (9  -  1)  I  -  4pz  +  4pzln(l/z)  0  <  z  <  1 

3  ds  s2  |  s  -  -1 

Since  all  residues  are  valid  for  0  <  z  £  1,  fz(z)  is  equal  to 

the  sum  of  R^,  R^ ,  and  R^.  Summing  these  terms  and  rearranging,  fz(z) 
is  then  given  by 


W,  this  becomes 


fZ,W<*.«>  1  “ 
and,  since  XlfXj  and  YltYj 

M3l,S2t  fZ,w(z,w)  >  - 


n  a .  81_1  a(  s2*1 

Et  a  (x,  )  (y,  )  1 

i-i 


independent  for  all  1  f  j, 


n  a  j  s  ^  -a  ^  b  ^  s  2  b  ^ 


*1  E[ 


1-1 


a^s^-aj  b^S2~b^ 


/  /  Mx^.y^x,  y,  ‘dx  dy. 

1-1  0  0 


1M  "i‘l«1-»l+l'bls2-bi+1<  f‘(*l,Tl>  ’ 


Applying  the  inverse  Mellln  transform,  (2.27)  follows. 

Example  2.3:  Consider  the  Kellogg-Barnes  I  distribution  given  by 


fXfY<*.y> 


4a0+1  ,Jl  ^  2*8  „-a(x2  +  y2) 

i.r<fl+l)  Cx  +7  ’  6 


x,y  >  o 


2  2 

Letting  r  -  x  +y  ,  the  Mellin  transform  is  given  by 


m<s1,s2) 


g+1  •  ir/2 

4a  r  r 

*TW> 


/  /  -~ar 
0  0 


S1+S9+2B-1  S1-1  S7-1 

r  L  c  cos0  1  sin0  L  d9dr 


a0+1r(s ,/2)r(s2/2)  -  Sl+s2+2e-2  2 

"  itr($+i)r(sl/2+s2/2)  ;0  r  e  2rdr 

Letting  u  -  r^  and  du  -  2rdr,  then 


<rT1T<s  /2)r(s,/2)  -  3,/2+3,/2+e-i 

(  l‘s2>  “  itr(B+l)r(s1/2+s2/2)  ;0  u  e 


ar(31/2)r(s2/2)r(B+s1/2+s2/2) 

s, /2+s?/2 

ur(8+l)r(s1/2+s2/2)a 


Now  suppose  the  bivariate  probability  density  function  of 

z  “  Xl/X2»  W  »  Yj/Y2  is  desired  where  X^Y^  and  X2»Y2  are  distributed 
according  to  the  Kellogg-Barnes  I  distribution  given  above.  Suppose 

further  that  ^(x^.y^)  has  parameters  «  1,  8^  ■  0,  and  f2(x2,y2) 
has  parameters  a2  ■  l,  82  ■  8.  Substituting  the  respective  parameter 
values  into  the  Mellin  transform  for  the  Kellogg-Barnes  I  distribution 
and  from  Theorem  2.5,  a^  ■  b^  ■  1,  a2  *  b2  *  -1,  the  probability 
density  function  of  Z,  W  is  given  by 

fz>w(z,w)  -  M2_1[^MslfS2t  *!<*!. Yl>  >  m2-s1,2-s2{  f2<x2»y2>  }] 


_x  rr(s1/2)r(s2/2)r(i-s1/2)r(i-s2/2)r(B+2-s1/2-s2/2)  I 
2  L  "2r(B+l)r(2-Sl/2-s2/2)  J 


r<8+2)  „  -i[  r(2)r(s1/2)r(s2/2)r(i-s1/2)r(i-s2/2)r(8+2-s1/2-s2/2)i 
,2r(8>1)  m2  L  r(8+2>r(i)r(i)r(2-s1/2-s2/2)  J 


Using  property  (2.9)  for  a  ■  b  »  2,  then  the  density  function  for  Z,W 


above  can  be  rewritten  as 
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f  (  2  -  *(3+i)  M  -i  T r(2)r(s1)r(32)r(i-81)r(i-82)r(B+2-»1-s2>i 

z‘mCz  ,2  2  L  r(S+2)r(i)r(i)r(2-srs2)  J 


Using  Equation  (2.15)  for  b  ■  b'  *1,  a  ■  8+2,  and  c  *  2,  the  inverse 
may  be  found  directly  and  is  given  as 

fZ,W<z’w>  *  ■■■" t1  ~  F1(3+2;1,1;2;-22,-w2) 


where  is  Appell's  hypergeometric  function  of  two  variables  as 


defined  in  Appendix  B.  The  series  converges  when  |x|  <  1  and  |y|  <  1. 
For  other  values  of  x  and  y,  the  function  can  be  evaluated  by  the 
usual  methods  of  analytic  continuation. 

Except  for  the  normalization  constant,  the  results  given  above 
are  identical  to  the  results  derived  by  Fox  (10).  The  distributions 
used  by  Fox  were  defined  over  the  range  -*•  <  x,y  <  •. 

By  applying  Theorem  2.5  followed  by  Theorem  2.4  or  Equation 
(2.23),  a  general  theorem  for  finding  the  distribution  of  a  random 
variable  which  is  the  product,  ratio,  or  power  of  an  arbitrary  number 
of  bivariate  random  variables  can  be  proved. 

Theorem  2.6;  If  X^,Y^;  X2,Y2;  ...;  Xn,Yn  are  n  pairwise  independent 
continuous  random  variables,  X^Yi  dependent  for  all  i,  with 

probability  density  functions  f  1(xl,y2),f2(x2,y2) . WV* 

Xi’^i  ^  ®  i“l,2,...,n  ,  then  the  probability  density  function  of  the 
random  variable 


f°r  ai,l>i  rational,  i*l,2,...,n,  is  given  by 


[n 

i-i  M»i«-«1+i.b18-bi+i{  fi<xi*yi>  z  >  0 


M*)  - 


(2.28 


otherwise 


The  proof  of  Theorem  2.6  follows  directly  from  application  of 
(2.23)  to  Theorem  2.5. 

The  advantage  of  these  theorems  is  that  all  the  theorems 
developed  by  Carter  (3)  are  special  cases  of  these  theorems. 
Therefore,  these  theorems  provide  the  techniques  for  finding 
distributions  of  random  variables  which  result  from  products, 
quotients,  and  powers  of  an  arbitrary  number  of  independent  and 
dependent  random  variables.  The  only  restriction  is  that  any  one 
random  variable  is  correlated  to  at  most  only  one  other  random 


variable. 


CHAPTER  3 


The  H-f unction 

3. 1  General  Remarks 

Tne  H-function  was  first  introduced  by  Fox  in  1961  as  a 
symmetric  Fourier  kernel  to  the  G-function  of  Meijer  and  was  used 
extensively  in  physics  and  engineering.  Carter  (3)  demonstrated  the 
importance  of  this  function  in  statistical  applications  when  viewed  as 
a  probability  distribution.  The  reasons  for  this  importance  are 
two-fold.  First,  the  H-function  is  the  most  general  special  function, 
containing  most  of  the  other  special  functions  as  special  cases. 

Thus,  anything  accomplished  with  the  general  form  for  the  H-function 
is  valid  for  all  special  cases.  This  allows  the  user  to  solve  a 
problem  for  a  large  class  of  functions  with  a  single  derivation. 

The  second  advantage  to  H-functions  is  readily  seen  in  the 
following  sections.  The  properties  of  the  H-function  are  such  that 
they  are  reduced  to  simple  adjustments  of  given  parameters.  The 
simple  parameter  changes  needed  to  find  the  Mellin  transforms  or  the 
derivatives  of  an  H-function  are  trivial  compared  to  performing  these 
same  operations  for  various  special  cases.  Indeed,  the  derivative  of 
an  H-function  is  another  H-function. 

Carter  (3)  used  these  properties  to  show  that  products, 
quotients,  and  rational  powers  of  Independent  H-function  variates 
yield  a  random  variable  which  also  follows  an  H-function  distribution. 
These  results  provide  a  robust  method  for  determining  algebraic 


combinations  of  independent  random  variables  and  strong  motivation  for 
extension  of  this  theory  to  bivariate  distributions. 

A  G-function  of  two  variables  was  given  by  Sharma  (34)  in  an 
attempt  to  generalize  classes  of  functions  of  two  variables.  Since 
the  Il-function  given  by  Fox  is  not  a  special  case  of  this  function, 
several  workers  have  extended  the  univariate  H-function  and  called  it 
an  H-function  of  two  variables.  The  H-function  of  two  variables 
contains  as  special  cases  most  of  the  known  functions  of  one  and  two 
variables,  Appell's  functions,  G-function  of  two  variables,  Whittaker 
functions  of  two  variables,  H-function  of  one  variable,  product  of  two 
H-f unctions,  etc. 

Chapter  4  presents  the  bivariate  H-function  distribution,  a 
bivariate  probability  function,  expressed  in  terms  of  an  H-function 
times  an  appropriate  constant.  Some  of  the  classical  bivariate 
distributions  are  shown  to  be  special  cases  of  the  bivariate 
H-function  distribution.  Methods  for  solving  resultant  distributions 
derived  from  products  and  quotients  of  dependent  H-function  variates 
are  presented.  The  important  result  is  that  such  combinations  result 
in  a  distribution  that  is  an  H-function  distribution  of  one  variable. 
Finally,  methods  for  handling  products  and  quotients  of  dependent 

H-function  variates  from  two  bivariate  H-function  distributions  which 
are  pairwise  independent  are  given.  Such  combinations  result  in 


distributions  which  are  also  bivariate  H-function  distributions. 


3.2  Definitions.  ( 5: 32; 7: 93; 3: 35; 14: 25; 69: 37; 12) 

Although  there  are  slight  variations  in  the  definition  of  the 
H-function,  this  paper  uses  the  H-function  of  one  variable  given  by 
Cook  (5:32)  and  a  slight  modification  to  the  definition  of  the 
bivariate  H-function  given  by  Goyal  (69:37).  The  H-function  of  one 
variable  may  be  defined  by 


H(z)  =*  H  [  z  :  (b.,0.)  ;  (+.,♦.)  ] 
P,Q  J  J  J  J 


a  r(*  S)  a  ra-a.-o.s) 

1  ,  j-i  J  J  j-i  2  3  -s  . 

_  /  Z _ _ _ _  z  ds 

2*  i  C  P  Q 

a  r (0,40.*)  .1  ro-*.-*,s) 

j»N+l  J  J  j*M+l  J  J 


where  C  is  a  contour  in  the  complex  s-plane  running  from  w-i<»  to  w+i« 
The  following  assumptions  are  made. 

(i)  M,  N,  P,  Q  are  integers  such  that  0<M<Q  and  OOKP 

(ii)  parameters  are  complex  numbers  and  0j,  j  are 

positive  real  numbers 

(iii)  empty  products  are  defined  to  be  equal  to  unity 

(iv)  all  poles  of  r($j+$js)  lie  to  the  left  of  C,  and  all  poles 
of  r(l-bj-Ojs)  lie  to  the  right 

The  bivariate  H-function  is  defined  by 
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M^,N^,M2,N2,M3,N3  x  (a^,A^)  ;  (C|,C^) 
H[x,y]  =  H 

» ^2  * ^2  * ^ 3  * ^3  y 

(b1,Bi)  ;  (di,D1) 


1  s  1  s2 
■■ —  2  ^  ^  *^ ( ~s 3 )x2 ( — s2 )x3( _s 3~s2 )x  y  ds^d^  (3.2) 

(2xi)^  Cx  C2 


where  an  empty  product  is  denoted  by  unity,  is  a  contour  in  the 
complex  s^  plane  running  from  h-i®  to  h-H®,  C2  is  a  contour  in  the 
complex  Sj  plane  running  from  w-i«  to  w+i®,  and 


"1  "1 

.1  T(bi-Bis1)  il  r(l-aj.+A.s1) 

*l(-*i)  "  1  1 

p_ 

11  r(ai-Als1)  .1  r(l-bi+Bis1) 


2  n2 

•i  r(di-D1s2)  ,i  r(i-c1+c1s2) 

x2^_s2^  m  1  1 _ 

^2  ^2  ” 
a  r(c1-cis2)  a  r(i-d1+D1s2) 


“3 

a  r(ei-Ei(s1-t-s2))  a  r(i-fi+Ft(s1+s2)) 
x3(-s1-s2)  -  1 _ 1 - 

^3  Q3 

a  r(fi-F1(s1+s2))  il  r(l-et+Ei(s1+s2)) 


(3.3) 


(3.4) 


(3.5) 


The  following  assumptions  are  made 


(i)  M^,  ,  P^,  Qit  i*l,2,3  are  non-negative  integers  such  that 

°lNi£Pi*  °iMi<Qi»  i-1.2.3 

(ii)  parameters  a^t  b^,  c^,  dt,  e^  are  real  or  complex  and 
parameters  A,,  Bj,  Cj_,  Dlt  Eit  Ft  are  real  positive  numbers 

(iiia)  poles  of  r(bi-Bisi)»  (i-l.-.M^),  I’(ei-Ei(s1+s2)) , 
(i=l...M^)  He  to  the  right  of  and  poles  of  r(l-a^+A^s^), 

(i-1-..Nj^),  r(l-fi+F1(s1+s2),  (i-l...N3)  lie  to  the  left 

(iva)  poles  of  IXd^-D^Sj),  (i*i...M2),  r(ei-Ei(s1+s2),  (i*l..«M3) 
lie  to  the  right  of  Cj  and  poles  of  T( 1-Cj+Cj.s2) »  (i*l-..N2), 
r(l-fi+Fi(s1+s2)),  (i*l...N3)  lie  to  the  left. 

If  — s ^  is  substituted  for  s^  and  -s2  for  s2  and  recognizing 

that 

a  b 

/  f(-s)d(-s)  *  /  f(-s)ds 
b  a 

and  perform  that  operation  twice,  then  H[x,y]  may  be  redefined  as 

.  -si  -S2 

H[x,y]  -  — - — ,  /  I  x1(s1)x2(s2)x3(s1+S2)x  7  dslds2  (3*6) 

(2xi)z  Cj.  c2 

where 


■I  r(b1+Bta1)  it  rCl-aj-AjSj^) 


r (a.+A.s.) 

n 

r(i-bt- 

-BiSl) 

l+1 

ml+i 

»2 

r ( d . +D . s  ^ ) 

n2 

n 

ii 

r(i-c.' 

“C< So) 

it  Iacj+c^)  it  r(i-d{-D{s,) 


f  (ej+E^(sj+s2))  *I  r(l_f£-F£(s3+s2)) 


jI  rff^Fj^Csj+s,))  it  r(l-ei-Ei(s1+s2)) 

N-j+1  M^+l 

Assumptions  (i)  and  (ii)  remain  unchanged  and  assumptions  (iiia)  and 
(iva)  are  changed  to 

(iiib)  poles  of  r<b1+Bis1),  (i-l..^),  r(e1+Ei(s1+s2>) , 
(i*l...Mj)  lie  to  the  left  of  and  poles  of  TCl-a^-A^s^), 

(i«l...N1),  r(l-f1-Fi(s1+s2)),  (i-l...N3)  lie  to  the  right. 

( ivb)  poles  of  r(d1+Dis2),  (i-l...M2),  I^ej+E^Sj+s^), 

(i»l...M^)  lie  to  the  left  of  C2  and  poles  of  TCl-c^-C^Sj), 

(1-1. ..N2),  r(l-f1-Fi(sl+s2)),  (i-l...N3)  lie  to  the  right. 

This  second  form  of  the  definition  has  an  advantage  in  that  it 
is  of  the  form  of  a  double  Mellin  integral  transform  inversion 
integral.  Form  (3.6)  of  the  H-function  definition  i9  used  hereafter. 


because  of  the  direct  relation  for  the  Mellin  transform. 

Form  (3.6)  of  the  H-function  definition  is  useful  in  that  its 
Mellin  transform  is  x^(s^)x2(s2>x-j(s^+S2)»  It  Is  clear  from  this 
definition  that  the  bivariate  H-function  can  represent  functions  whose 
Mellin  transforms  have  terms  of  the  form  x3(S]+s2)  or  X3(-s^-S2)» 
Practically  speaking,  it  is  very  likely  that  a  great  many  functions 
will  have  Mellin  transforms  with  terms  of  the  form  x^(s^-S2)  or 
x3^s2”sl^*  *s  Possible  t0  modify  definition  (3.6)  to  accomodate 
such  functions. 

Using  definition  (3.6),  if  -s^  is  substituted  for  s^,  then  an 
equivalent  representation  for  (3.6)  is 

H[x,y]  =•  — - — j  ^  *  xl(-s1)X2(s2)x-j(-s1+s2)x  1  y  2  ds1ds2  (3.10) 

(2iti)z  ci  c2 

Similarly,  if  -S2  is  substituted  for  S2  in  (3.6) 

1  ~S1  s2 

H[x,y]  -  — i — x  f  f  x. (s, )x,(-s»)x^(s1-s2)x  y  ds.ds,  (3.11) 

(2iti)2  CL  C2 

where  the  poles  of  the  x2  and  the  X3  terras  are  interchanged  about  the 
Cj  axis  in  the  S2  plane  for  (3.11)  and  the  poles  of  the  x^  and  the  X3 
terms  are  interchanged  about  the  axis  in  the  s^  plane  for  (3.10). 

From  this  it  is  clear  that  the  mathematics  of  the  H-function 
can  handle  functions  whose  Mellin  transforms  have  terms  of  the  form 
*3(81 —3 2 )  or  X3(S2~Si)«  Definition  (3.6)  can  now  be  modified  to 
accomodate  such  functions. 


Definition:  A  bivariate  functional  whose  Mellin  transform  is  given  by 
x^(s^)x2(-S2)x;j(s^-S2)  has  the  bivariate  H-function  given  by 


H[x,y] 


Ml,M1,M2,N2,M3,N3 

2H 

^l’Ql’^2  ’^2 » ^*3 » Q3 


<ei,Ei) 

x  (a^A^  ;  (cj.Cj.) 
y  Cfi*Fi> 

<bi,Bt)  ;  (d^D^ 


1  ~S1  s2 

”  - •-  *  f  X1(s1)x2(s2)x3(si+s2)x  y  ds.ds2  (3.12) 

(2iri)2  Cx  c2 

where  x1(s1),  X2(s2),  and  X3(sj+s2)  are  defined  by  (3.7),  (3.8),  and 

(3.9).  Assumptions  (i),  (ii),  (iiib),  and  (ivb)  hold. 

Henceforth,  jH[x,y]  shall  be  used  to  denote  definition  (3.6) 

and  2H[x,y]  to  denote  definition  (3.12).  The  relationship  between 

s2 

jH[x,y]  and  2H[x,y]  is  readily  derived  by  observing  that  y  may  be 
”s2 

rewritten  as  (1/y)  in  (3.12).  Comparing  (3.12)  with  this 
substitution  to  (3.6),  it  is  readily  seen  that  the  following 
relationship  holds. 

2H[x,y|  -  jHfx.l/yl  (3.13) 

3.3  Properties:  (14: 24-2 5 ; 69 : 3 9-4 1 ; 12) 

Making  a  change  of  variable  substitution  in  (3.2)  yields  the 


following  identities  which  are  useful  in  bivariate  manipulations: 


V  *.  V 


*.*  V  O 


**  4 


H[x,y] 


Ml,N1,M2,N2,M3,N3 

1 

P1  *  Ql » p2  *  ^2  ’  p3  *  ^3 


(ei»El> 

1/x  (a^A^  ;  (cj^.Ci) 

1/y 

(bi,Bi)  ;  (di,Di) 


^1»P1*Q2’P2'Q3»P3  F 

(l-aj.Aj)  ;  (l-c^.C^) 


For  k  >  0,  the  following  is  true. 


<ei,kEi> 

M1,N1,M2,N2>M3,N3  xk  (^.kA^  ;  (Ci.kq) 

k*  h 

P1  »Q^ » P2  *^2  ’  P3»^3  y 

<bi.kBi>  ;  Cdi,kDt) 


(e1.E1) 

MpN^,M2,N2,M3,N3  x  I  (a^.A^)  ;  (c^.C^) 


pl»Ql»p2»Q2»p3»Q3 


(fl.pi) 

(b4  »B, )  ;  (d, ,D» ) 


(3.14) 


(3.15) 


VtV»rAW>' 


Making  use  of  the  Mallin  transform  property  M  {  xmynf 

sl,s2  X 

Ms1+m,s2+n^  fX,Y^x,y^  the  property  is  obtained 

(ei.Ei> 

m  n  mi»N^.M2,N2,M2,N^  x  (a^.A^)  ;  (c^.C^) 

x  y  XH 

^1 ’^1  * P2 *^2 ,P3 *^3  y  (^i»Fi) 

(bi.Bi)  ;  '*,,0^ 


C  ®  j+E^(m+n )  9  ) 

M1»n1*m2*n2«m3»n3  x  <ai+Aim»Ai)  5  (ci+C1n,C1) 
f,l’^l’^2*^2,^3,^3  y  (F  j+F^(nH-n)  ,F^) 

(b1+Btm,Bi)  ;  (dj+Djn.D^ 


m  n  u  M1,!,1’M2,N2,M3’N3  x  (ai*Ai)  5  (ci*Ci) 

x  y 

PltQl.p2«Q2»p3»Q3  y  ^fi»Fi^ 

(b^,B^)  j  (d^,D^) 


(ei+Ei(m-n),Ei) 

^1  ’ ^1  ’^2 ’ ^2 * ^3 ’ ^3  x  >  (cj_Cjn,Cj) 

Pl»Ql»p2»Q2,P3*^3  y  (f i+F^Cm-n) ,F^) 

(bl+Bitn,Bi)  ;  (d^-D£n,Di) 


3.4  Mellin  Transform:  (3 : 37; 5: 35; 7 : 102 ; 12) 


Definition  (3.1)  of  the  H-function  is  exactly  that  of  a  Mellin 
transform  inversion  integral  so  that  the  Mellin  transform  of  the 
H-function  is  directly  given  as 


Mg(H(cz)} 


M  N 

*1  r  ( 1»  s )  H  r(l-0.-0.s) 

c~s  j-1  1  j»l _ 1  1 


P  Q 

*i  r(9.+Q.s)  .i 

i-N+1  i-M+1 


r(i-4>i-*is) 


(3.18 


Form  (3-6)  of  the  bivariate  H-function  definition  is  exactly 
that  of  a  double  Mellin  transform  inversion  integral,  so  that  the 
Mellin  integral  transform  of  jH{x,y]  is  directly  given  as 

-si  -s  2 

Ms1»s2  *  lHfglx,82yl  *  *  x1(s1)x2(s2)x3(s1+s2)g1  g2  (3.19 

where  x^(s^),  x2(s2),  and  x3(Sj+s2)  are  defined  by  equations  (3.7), 
(3.8),  and  (3.9)  respectively. 

Under  definition  (3.12)  for  2H[x,y]  and  assuming  convergence 
of  the  integral  in  the  definition,  the  Mellin  transform  can  be  found 

by  interpreting  the  bivariate  H-function  as  the  inverse  Mellin 

“S1  ~a2 

transform  of  the  coefficients  on  x  y  •  Then 

l  ,  "S1  -s2 

- - J  *  xl(s1)x2(-s2)x3(s1-s2)(g1x)  (82y>  dslds2 

(2*i)Z  C.  C- 


2HtSix,g2y 1 


saq 


J  r<e1+Ei(s1-s2))  il'  rd-fj-F^Sj-Sj)) 
x3 ( 3 1 ~s 2 )  "  1  l _ 

11  r(fi+Fi(s1-s2))  il  r(l-e  -Ei(s1-s2)) 

M3+l 


(3.20) 


Using  the  definition  of  the  Mellin  transform,  one  can  express 
2H[g1x,g2y)  in  the  form 


-1  ~si  ~s2 

2Hf8ix.g2yl  "  M2  t  xi(si)x2(-s2)x3(s1-s2)g1  g2  ] 


where  M2  is  the  inverse  operation  for  the  doable  Mellin  transform  as 
defined  by  Fox  (10).  It  follows  that 


Mslfs2{2Hl8lx-82yl}  *  x1(sl)x2(-s2)x3(s1-s2)g1  Slg2  82 
3.5  Special  Cases:  ( 14: 26-28; 16: 569-571 ; 12) 


(3.21 


From  the  results  of  Reed  (16)  on  doable  Mellin  transforms  of 
Appel Is  functions,  the  following  special  cases  for  the  H-function  of 
two  variables  may  be  derived. 


(1^,1) 

1,1, l, 1,0,1  x  I  (l-b,l)  ;  (l-b’,1) 


1,1, 1,1, 1,1 


U-a,l) 

(0,1)  ;  (0,1) 


r(a)r(b)r(b’)F1(a.b.b,:c:-x.-y) 

r  (c) 


(3.22) 


Kampe*  de  Feriet's  function. 


1 ,  B,  1,  B,  0,  A  x 

1H 

B,D+1,B,DH,A,C  y 


(l-Ci.l) 

(1-b^l)  ;  (1-b.M) 
(0,l),(l-dlfi)  ;  (0,l),(l-d1' ,1) 


A  B 

■i  r(at)  :i  r(b1)r(b1I ) 

=»  i  l 

C  D 

-i  r<c .)  n  rcd^rcd^) 

i  l 


ai;-*-aA 
bi»bi' ; • • *bB*bB' 


ci ; 


i»  •  •  •'-c 


' * • • •^d,^d' 


-x,-y 


(3.26) 


The  identity  given  by  Mathai  and  Saxena  is  of  a  slightly 
different  form  due  to  the  fact  that  they  used  the  definition  of  the 
form  of  (3.2).  The  identity  above  results  from  the  definition  for  the 
H-function  of  two  variables  given  by  (3.6). 

If  M-j“N3=*P3*Q3“0,  the  H-function  of  two  variables  breaks  up 
into  a  product  of  two  H-functions. 


m1,n1,m2,n2,o,o  x 
lH 

Pl,Ql»p2»Q2»°'°  y 


(af,Af)  ;  (c^,C^) 


v 


a-"!  r 

H  x 


(ai»Ai>  j  hM2*N2 


(bi»®i)  -1  p2»$2 


(ci,Ci) 

(di,VL) 


(3.27) 


-  H1(x)*H2(y) 

Here,  H^(x)  and  H2(y)  are  univariate  H-functions  as  defined  by 
(3.1).  Similar  results  hold  for  2R{x,y]. 


^l«Nl,M2jN2,OfO  x  (a^,A^)  ;  (c^,C^) 


pi.Qi.p2.Q2.o.°  y 


mi»ni  r  i  <ai»Ai^i  n2*m2  r  i  (i_di*Di) 

H  x  *  H  y  I 

Pl»Ql  _  (k^*®i)_  Q2*P2  - 


(3.28) 


From  (3.27)  and  (3.28),  it  is  clear  that  for 
i*2,3,  the  bivariate  H-function  reduces  to  a  univariate  H-function  as 
defined  by  (3.1).  It  then  stands  to  reason  that  distributional 
analysis  of  products  and  quotients  of  H-function  variates  in  the 
univariate  domain  is  a  special  case  of  distributional  analysis  of 
H-functions  of  higher  order.  This  fact  shall  be  demonstrated  in 
Chapter  4. 

Example  3.1:  Consider  once  again  Morgensterm' s  bivariate  uniform 
distribution.  The  double  Mellin  transform  was  given  in  Example  2.2  as 


1  p(sx  -  l)(s2  -  1) 

MX,Y<S1’S2>  -  —  +  Sls2(Sl  +  "l ) (s 2  +  1) 


Noting  that  (s  +  n) 


?(s  +  n  +  1) 
f(s  +  n) 


MX  y^sl,s2^  may  ^e  rewritCen  *n  terras  of  gamma  functions. 

r(s1)r<s2)  pr (s ^  )r(s2 )r( s j )r( s2 ) 

mx,y(si,s2)  ■  r(s1+i)r(s2+r)+  r(s1+2)r(s2+2)r(s1-i)r(s2-i) 


While  this  form  is  not  immediately  expressible  as  a  single  H-function 
of  two  variables,  it  may  be  redefined  as  a  sum  of  two  H-f unctions  of 
two  variables. 


fx,Y<x»y>  *  Hit  x»y  1  +  U2r  x,y  i 


where 


Hlt  x,y  ]  *  jll 


1,0, l, 0,0,0  x  (1,1)  ;  (1,1) 


1,1, 1,1, 0,0  y 


(0,1)  ;  (0,1) 


h2(  x.y  1  *  Pi» 


2, 0,2, 0,0,0  x  (2,1), (-1,1)  ;  (2,1), (-1,1) 


2, 2, 2, 2, 0,0  y 


(0,1), (0,1)  ;  (0,1), (0,1) 


In  this  instance,  it  must  be  remembered  that  H^[  x,y  ]  and 
Hjf  x,y  ]  are  not  H-functional  representations  of  density  functions, 
but  are  two  general  H-functions  of  two  variables,  the  sum  of  which  is 
a  bivariate  density  representation. 

While  it  is  desirable  to  formulate  fx  Y(x,y)  as  a  single 
H-f unction  of  two  variables,  the  representation  above  is  still 
valuable  in  that  given  Z  «  XY  or  Z  =  X/Y,  f2(z)  may  be  represented  as 
the  sum  of  two  H-functions  of  one  variable  which  can  be  derived  from 
HjJ  x,y  ]  and  x,y  ].  Numerical  inversion  of  HjJ  z  ]  and 
iljt  z  ]  can  be  accomplished  by  methods  presented  by  Eldred  and  Cook. 
Numerical  evaluation  of  fg(z)  can  then  be  derived  by  appropriately 
summing  the  inversions  of  H^[  z  ]  and  H2(  z  ]. 


CHAPTER  4 


The  H-function  Distribution 

4.1  General  Remarks 

In  this  chapter  a  new  bivariate  probability  density  function 
based  on  the  H-function  of  two  variables  is  introduced.  The  new 
distribution,  called  the  bivariate  H-function  distribution,  includes 
as  special  cases  many  of  the  more  common  bivariate  distributions  -  t! 
bivariate  gamma,  the  bivariate  beta,  and  the  bivariate  Cauchy.  Threi 
new  bivariate  distributions,  called  the  Kellogg-Barnes  distributions 
are  also  shown  to  be  special  cases  of  the  bivariate  H-function 
distribution.  Also,  by  extension,  all  the  univariate  H-function 
distributions  are  special  cases. 

Formulas  for  finding  the  moments  of  the  bivariate  H-function 
distribution  and  the  normalizing  constant  are  given.  The  cumulative 
of  the  bivariate  H-function  distribution  is  shown  to  be  another 
H-function,  a  property  that  is  not  shared  by  other  bivariate 
distributions.  Examples  are  given. 

4.2  Definitions : 

Definition:  Consider  a  random  variable  Z  with  probability 
density  function  given  by 

C k  H(cz) ,  cz  e  S 

fZ(z)  -  )  (*• 

(_  0,  otherwise 

where  H(cz)  represents  a  univariate  H-function  as  defined  in  section 
3.2,  k  and  c  are  real  constants  such  that 


and  S  is  a  subset  of  the  positive  real  values  u  for  which  H(u)  is 
convergent.  The  random  variable  Z  will  then  be  called  an  H-f unction 


variate  or  a  random  variable  with  an  H-function  distribution 
(3:4l;7:103;5:84;l7:20Q). 

Definition;  Consider  the  random  variables  X,  Y  with  joint 
probability  density  function  given  by 

{krH[giX,g2y]  ,  Sjx  e  Sx  ,  g2y  e  S2 

(4.2) 

0  ,  otherwise 

where  rH[gjx,g2]  represents  an  H-function  as  defined  in  section  3.2, 

k,  g^,  and  g2  are  real  constants  such  that 
00  00 

/  /  fx  Y(x,y)  dxdy  •  1  , 

0  0  *,Y 

and  S^,S2  are  subsets  of  the  positive  real  values  u,v  for  which  H(u,v) 
is  convergent.  The  random  variables  X,  Y  will  then  be  called 
dependent  H-function  variates  or  random  variables  with  a  bivariate 
H-function  distribution. 

4.3  Special  Cases 

Since  the  univariate  H-function  distribution  may  be  expressed 
as  a  special  case  of  (4.2),  all  of  the  classical  univariate 
non-negative  probability  distributions  studied  by  Carter  (1972), 

Eldred  (1979),  and  Cook  (1951)  are  special  cases  of  the  bivariate 


H-function  distribution.  In  addition,  some  of  the  classical 
non-negative  bivariate  distributions  may  be  expressed  in  the  form 
(4.2).  Converting  a  probability  density  function  into  its  H-function 
form  is  accomplished  by  taking  the  Mellin  transform  of  the  density 
function  and  arranging  the  transform  such  that  it  is  products  and 
quotients  of  gamma  functions.  The  H-function  form  may  then  be 
identified  by  taking  the  inverse  Mellin  transform. 

( i)  McKay 1 s  bivariate  gamma  distribution 

p+q  ,  , 

fX  Y(x’y>  *  xP  (y*)q  e“ay  y>x>0 

r(p)T(q)  a,p,q>0 

Let  c  *  (aP+q)/(r(p)r(q)).  Taking  the  Mellin  transform  and 
integrating  first  with  respect  to  x  and  then  with  respect  to  y  yields 

*  y  —i  s  — i 

M(s,,s2)  «  cl  I  xp_1  (y-x)q_1  e-ay  xSl  y  2  dxdy 
0  0 

00  St-!  y  Si+p-2  . 

*  c/  y  e  ay  [/  x  (y-x)q  1  dx]  dy 

0  0 

Using  the  Mellin  transform  property  (2.8),  the  transform  above  becomes 

*  So-1  y  ,  , 

m(s1»s2)  -  c /  y  e-ay  [/  xz_1  (y-x)q_1  dx]  dy 

0  0  I  z«s1+p-i 

Using  (15:16,  if  2.20)  to  evaluate  the  inner  integral,  the  transform 


can  now  be  written 


74 


cr(q)r(Sj+p-l)  08  sl-t-s2+P+-q-3  _ay 


M(S1’32>  •  iT,1+^-l)  y 


e  dy 


Realizing  that  the  integral  is  just  the  gamma  function  and  replacing 


the  value  for  c  yields 


2  s  ^  s  £ 


H(»1,s2> 


r(p-l+s1)r(p+q-2+s1+s2) 

r(p)T(p+q-l+s1) 


fx  Y(x,y)  is  returned  by  taking  the  inverse  Mellin  transform  given  by 


(2.6). 


-  ffo  —l" 


2  i  rCp-l+s, )r(p+q-2+s,+s7)  -s,  -s, 

5 _  _  f  T  /  av^  r 


rcp+q-l+s^ 


(ax)  (ay)  ds  j. 


(p+q-2,1) 


2  1,0, 0,0, 1,0  ax  (p+q-l,l)  ; 


1, 1,0, 0,0, 1  ay 


(4.  )) 


(P-1,1)  ; 


(ii)  The  bivariate  beta  distribution 


r(p1+p2+p3)  pi-l  p2"l  P3-l 

r(Pl)r(p2)?(p-  >  *  y 


x,y>0  ,  x+y<l  ,  p1,p2,P3>0. 


Let  c  -  r(p1+p2+p3)/r(p1)r(p2)r(p3).  Taking  the  Mellin  transform  and 


integrating  first  with  respect  to  y  and  then  with  respect  to  x  yields 


L  1-X  p  -1  p--l  p  -1  g  -1  g  -1 

M(s1,s2)  -  cf  f  x  y  (1-x-y)  X  y  dydx 


1  (Sl+Pl-1)-1  (.2+P2-l)-l  P3-1  ,  , 

cj  x  [/  y  (1-x-y)  dy]  dx 

0  0 


Using  (15:16  #  2.20)  to  evaluate  the  inner  integral  and  substituting 
the  value  of  c,  the  transform  can  now  be  written 

r(p1+P2+P3)r(s2+P2-l)  1  (sj+pj-l)-!  s2+p2+p3~2 

M(s,  S_)  «  - -  /  x  (1-x)  dx 

2  r(pi)r(p2)r(s2+P2+P3-i)  o 


Using  (15:16  //  2.20)  once  again  to  evaluate  the  integral  yields 

r(Pi+P2+P3)r(sl+Pi"1)r(s2+P2~1) 


1  z  r(p1)r(p2)r(sl+s2+p1+p2+p3-2) 

The  ll-f unction  form  is  found  by  using  the  Mellin  transform  identity 
for  the  H-function  as  defined  by  (3.18)  or  (3.20). 

fx,Y<x-y)  *  m2_1[  M(s1.92>  i 


r(Pj+P2+P3)  1,0, 1,0, 0,0  x  ; 

r(P1)r(P2)  1  o,i,o, 1,1,0  y  (Pl+p2+p3-2,i) 

(Pi-i.i)  ;  (p2-i,n. 


(iii)  The  quarter  Cauchy  distribution 


fx>Y(x,y)  -  (c2  +  x2  +  y2)‘3/2 


x,y>0  ,  c>0 


2  2  2 

Letting  r  ■  x  +y  ,  x  ■  cos9,  y  ■  sin9,  and  dxdy  *  rdrdQ,  the  Mellin 
transform  can  be  written  as 


-  */2  Si -1  s?-l  _  _  -/9 

M( s.  ,s,)  -  —  /  /  (rcos9 )  (rsinQ )  (cz+rz)"J/ z  rdQdr 

1  "  00 


*  8 i ‘+'8 2-i  «  9  i/j  ^/2  Si- 1  s,-l 

.f£.  /  r  1  L  (c2+r2)  3/2  [/  cose  1  sin9  1  dG]dr 
11  0  0 


_  °°  Si+Sj  —  l  ~  n  ,  in 

B(s,/2,s,/2)  /  r  (c  +rZ)  iU  dr 

w  o 


where  B(u,v)  is  the  beta  function  and  is  given  by 


Using  (15:15  #  2.19)  to  evaluate  the  integral  and  Mellin  transform 
property  (2.8),  the  transform  becomes 


Sj+S2 


M(*i ,82) 


cr(s1/2)r(s2/2)  c  r(s1/2+s2/2)rC3/2-*i/2“82/2) 
2*r(.1-f<2)  c3r<3/2) 
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-1, 


fX,Y^x,y^  "  M2  *  m(si.s2)  1 


1,0,1, 0,1,0 


.H 


*  1  0,1,0, 1.1,1 


(0,1/2) 


/a  x 
/a  y 


(0,1/2) 

(0,1/2)  ;  (0,1/2) 


(4.6) 


(v)  Kellogg-Barnes  Type  II  distribution 


r  -  .  0  2  -ax  -  0y/x 

fX,Y(x*y)  "  Sa  e 


x,y>0 

a,0>O 


The  Mellln  transform  Is  given  as 


M(s  ,  )  .  B«2  /*  y'2"1  [/’  X*1’1  e-ax-0y/’t  dx]  dy 
0  0 


Using  (95:313  #  17)  to  evaluate  the  Inner  Integral  yields 


Si  /2 


-  0«2  y"2”1  2<Sy/a)  1  ‘  V  2<“6>”1/2  1  d” 


2  81/2  "  82"*"®1^— 1  i/2 

-  2Bc/(B/a)  /  y  K  [  2(a0y)1/z  1  dy 

0  91 


Sj/2  •  s^+2s2“l 


-  48a‘(B/a)  /  u  K  [  2/aB  u  J  du 

0  81 


where  Ky  is  the  modified  Bessel  function  as  defined  by  Erdelyi 


(95:371).  Using  (95:331  #  26)  to  complete  the  integration  yields 


2  ”sl”s2  -s2 


m(s1,s2)  -  6c*  o  6  r(Sl+s2)r(s2) 


Inverting  and  using  the  H-function  definition  (3.6)  the  H-function 


form  for  f^  y(x,y)  is  obtained. 


fv  f  ^(8i  |8o)  ] 


(0,1) 


8a2  XH 


0,0, 1,0, 1,0  ax 


0,0,0, 1,0,1  a0y 


5  (0,1) 


(vi)  Kellogg-Barnes  Type  III  distribution 


'*.*<«•*>  -  {&>  *c  * B:,y  • 


x,y>0 


a,8>0  ,  c>2 


The  Mellin  transform  is  given  by 


M(s1»82)  "  f^c)  f  x 


!®c  /  *8l+c  1  „-°x  ti  82  1 


e  x  [/  y  ‘  e  dy]  dx 
0 


„  I  “3  2 
ac8  2 


r(c)  r(82)  ;n  x 


*  8,-S-^C-l 

1  4.  ^ax  . 

X  e  dx 


-So  -s 


■  r(s2)r(c+8l-s2)(8/a)  o 


2  81 


(4.7) 


kl 
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Inverting  and  using  definition  (3.12)  or  using  the  Mellin  transform 
property  (3.20)  for  the  H-function,  the  H-function  form  for  fx  ^(x,y) 
is  obtained. 

™  ^2  ^  ^sl*s2^  1 

(c,l) 

-  ;  (l.D 


4.4  Moments  of  the  H-function  Distribution: 

Carter  (3)  showed  that  the  moments  of  the  univariate 
H-function  distribution  can  be  found  by  taking  the  appropriate 
derivatives  of  the  characteristic  function  of  the  univariate 
H-function  distribution.  To  use  this  approach  for  the  bivariate 
H-function  distribution,  the  double  Fourier  transform  of  the  bivariate 
H-function  must  first  be  shown  to  exist.  To  do  this  an  extension  to 
Prasad's  theorems  (57)  must  be  developed.  Fortunately,  while  the 
above  is  considered  beyond  the  scope  of  this  dissertation,  there  is  a 
simpler  method  for  obtaining  the  moment,  of  the  bivariate  H-function 
distribution. 

Using  the  notation  of  section  1.3,  the  nl,n2  ordered 
noncentral  moment  for  fx^Y(*.y).  anl,n2*  is  deflned  by 


»  -  jrpjk  ■ 


-_w  ■_> 


\v  7  VV  ' 


W  V/J'/ 


81 


anl,n2 


E[  x 


nlyn2 


] 


»  /  /  xnl  yn2  fx  Y(x,y)dxdy 
0  0 


where  E  is  the  expected  value  operator.  From  the  definition  of  the 
Mellin  transform,  it  is  clear  that 

si-1  S2~l 

Ms  a  t  fx  v(x.y)  *  ■  Ef  x  y  ) 

for  distributions  defined  for  x,y  >  0.  The  nl,n2  ordered  moment  for 
f x  Y(x,y)  may  then  be  obtained  from  the  Mellin  transform  of  the 
probability  density  function.  Specifically, 

anl,n2  *  Mslfs2{  fX,Y^y)  } 

Then  from  the  Mellin  transform  property  (3.18)  of  the  bivariate 
H-function,  the  nl,n2  ordered  moment  for  k^H{g^x,g2y)  is  given  by 

anl,n2  "  nl+T"  nS+1  *i<nl+l)x2<n2+1>x3<nl+ll2+2)  (4’9> 

*1  g2 

where  x^(u),  x2(v),  and  Xj(u+v)  are  defined  by  (3.7),  (3.8),  and  (3.9) 
respectively. 

Similarly,  from  the  Mellin  transform  property  (3.20)  of  the 


s. *nl+l 
s,-n2+l 


bivariate  H-function,  the  nl,n2  ordered  moment  for  k2Hl8lx»82y)  is 
given  by 


(4.10) 


anl,n2  "  ~nl+T  n2-Hl  xiCnl+1)*2(-n2_1)x3(nl~n2) 

81  g2 

where  x^u),  x2(-v),  and  x3(u-v)  are  defined  by  (3.7),  (3.4),  and 
(3.19)  respectively. 

Following  the  procedures  outlined  in  section  1.3,  (4.9)  and 

o 

(4.10)  above  can  be  used  to  find  u  ,  p  ,  y  ,  o  ,  and  o  .  From 

x*  y*  xy*  x  y 

these  values,  the  covariance  and  correlation  for  X  and  Y  may  be  found. 
Example  4.1:  From  (4.3),  McKay's  bivariate  gamma  distribution  can  be 
represented  as  a  bivariate  H-function  distribution. 

(p+q-2,1) 

(p+q-1,1)  ;  — 


(P-1,1)  ; 

Using  (4.9),  may  be  found  by  setting  nl  ■  1  and  n2  «  0. 

,2  r<p+i)1Wn 
aJ  r(p)r(p+q+l) 


■  p/a 


Similarly, 


a3  r(p)r(p+q) 


-  ( p+q ) / a 


Using  the  identities  o  ■  a„  _  -  a,  «  and  a  *  a_  -  -  an  ,  ,  the 

x  2,0  1,0  y  0,2  0,1 

variances  may  be  found. 

o  2  ,  a2  U21 -DllSHSm  -  (p/a)2 

X  a4  r(p)r(p+q+2) 

»  P(P+1)  _  £2 

a2  a2 

*  p/a2 

Similarly, 

a  2  .  a 2  r(p)r(p-*i+2)  _  (p+q )2 
y  a4  r(p)r(p+q)  a2 

,  (p+q+l)(p+q)  _  (p+q)2 

a2  a2 

*  (p+q)/a2 

Using  (1.10),  the  covariance  may  be  found  by 
cov(x,y )  -  o1>x  -  o1|Q  aQ>1 

,  a2  r(p+l)r(p+q+2)  _  p(p+q) 
a4  r(p)r(p+q+l)  a2 

■  p/a2 


Finally,  the  correlation  is  equal  to  the  covariance  divided  by  the 
standard  deviations  for  x  and  y. 
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*.  \ 


p(x,y) 


cov(x,y )/o  o 

x  y 


£  a_  a_ 

&  /p  •'p+q 


■  / p//p+q 


The  moments  computed  above  agree  with  those  given  by  Mardia 
(47:89)  for  the  bivariate  gamma  distribution.  The  moments  for  all  six 
of  the  special  cases  of  the  H-function  distribution  given  in  section 
4.3  were  computed  using  (4.9)  or  (4.10)  and  agreed  with  those  given  by 
Mardia  or  computed  using  the  equations  given  in  section  4.4. 

4 . 5  Cumulative  Distribution  Functions: 

4.5.1  Univariate  Cumulative  Distribution 

The  cumulative  distribution  function  Hc(z)  of  a  probability 
density  function  H(z)  is  defined  as 


H  (z)  -  /  H(u)du 
C  0 


Using  a  well-known  Mellin  transform  relation  from  Erdelyi 
h95:307),  Eldred  (7:139)  derived  an  expression  for  Hc(z).  Cook 
(5:103)  improved  on  this  expression  and  derived  the  form 


H  (z)  -  1  -  (k/c)  H 

c  f 


M+l.N  r  (e  +G.,0.),(1,1); 

cz  J 

P+i.Q+iL  (O.lM^j.^j) 


(4.11) 


1 
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Cook  then  showed  a  second  formulation  for  the  cumulative 

1  distribution  and  proved  that  the  cumulative  distribution  function  for 

i 

an  H-function  distribution  is  another  H-function.  Using  the  Laplace 
transform.  Cook  derived  the  following: 


» 

f 

f 


HCU> 


f  M,N+1 

,  (1,1), (6  +0  ,0  );  ' 
cz  1  J  J  J 

\  (k/c)  H 

P+1, Q+l 

(t>j+»j,»j),(°,1)  . 

all  <  1,  J-1,...,M 


M+l  ,N 

.  (~k/c)  H 

V  p+i *o+i 


cz 


(9j+0j,Qj),(1,1);" 

(O.lMtyWj.tj) 

if  any  j/* j  >  1,  j-1 


.,M 

(^•12) 


i 

i 


Hc(z)  can  also  be  computed  using  the  Mellin  transform  property 
for  integrals  given  by  Oberhettinger  (15:12)  and  Sneddon  (106:269). 


2  i 

M  {  /  f (u)du  >  -  -7  M  ,{  f(*)  ) 

s  o  s  »+i 

Letting  -1/s  ■  r(-s)/r(l-s)  results  in  (4.12)  case  I  and  letting  -1/s 
■  -r(s)/r(s+l)  results  in  (4.12)  case  II. 

Comparing  (4.11)  to  (4.12)  case  II,  it  appears  at  first  glance 
that  one  of  the  two  formulations  must  be  in  error.  The  Laplace 
transform  introduced  by  Carter  (3)  can  result  in  one  or  more  of  the 
poles  associated  with  the  ($j»*j)  overlapping  with  the  poles 


associated  with  the  new  (0,1)  term  in  the  numerator.  Under  such 


conditions,  no  contour  exists  to  properly  separate  the  poles*  To 
correct  this  problem,  Cook  (5:81)  eliminates  the  overlap  by  using  an 
equivalent  expression  for  the  T(l-s)  term  introduced  in  the 
development  of  the  Laplace  transform  of  the  H-function.  The 
replacement  of  this  equivalent  identity  results  in  (4.12)  case  II. 

Assume  _>  1  for  some  j.  Then  (4.12)  case  II  is  valid 

and  (4.11)  is  not  valid.  Further,  (4.11)  is  only  invalid  if  the  poles 
associated  with  r(l-0j-0j-O^s)  overlap  with  the  poles  associated  with 
the  new  term  T(s).  However,  these  poles  will  only  overlap  if 
(9j+0j)  1  for  some  j,  j  ■  1,...,N.  Now  consider  the  poles  of  the 

density  function.  If  this  condition  holds,  then  some  of  the  poles  of 
r($j+$js),  (for  -♦j/*j),  will  occur  at  values  of  s  >  1  and  some  of  the 
poles  of  r(l-0j-OjS),  (for  6j+Gj  >  1),  will  occur  at  values  of  s  £  1. 
By  definition  (3.1)  such  an  overlap  of  poles  is  not  allowable. 
Therefore,  the  condition  of  will  never  occur  in  H-function 

distributions  and  the  cumulative  distribution  is  given  by  (4.12)  case 
I  only.  These  results  are  summarized  in  the  following  theorem. 

THEOREM  4.1:  If  kH(cz)  is  an  H-function  probability  density  function 
as  defined  by  (4.1),  then  <  1,  j  ■  1,...,M. 


From  the  discussion  above,  it  would  seem  that  it  should  also 
be  true,  for  the  same  reasons,  that  <  !•  While  a  check  of  all 

the  special  cases  listed  by  Cook,  (6:85-87),  support  this  idea,  it  is 
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noC  possible  Co  go  from  (A. 11)  to  (A. 12)  and  prove  chat  this  is  true 
in  general. 

It  should  also  be  pointed  out  that  the  improved  Laplace 
transform  for  the  H-function  as  given  by  Cook  (6:82)  is  not  in  error. 
Specifically,  it  does  not  necessarily  hold  that  <  1  , 

j  ■■  1,...,N,  for  the  H-function  in  general.  It  is  only  true  if  the 
H-function  is  a  density  function. 

A. 5.2  Bivariate  Cumulative  Distribution 

The  cumulative  distribution  function  Hc(x,y]  associated  with 
the  probability  density  function,  H[x,y]  is  defined  by 

x  y 

Hc[x,y]  *  /  /  H[u,v]  dudv 
0  0 


H.[x,y]  can  be  obtained  by  direct  integration  or  through  the 
use  of  the  Hellin  transform  of  H[x,y}.  The  latter  method  is  most 
often  preferable  as  it  avoids  the  necessity  of  evaluating  H[x,y]  in 
order  to  derive  Hc[x,y]. 

The  use  of  the  preceding  procedures  in  evaluating  Hc[x,y]  is 
made  possible  by  noting  from  equation  (2.1A)  that 


M  { 
sl»s2 


1-Hc[x,y]  } 


(sl82>  Ms1+1,s2+1{  Htx’y>  } 


Inverting  the  equation  above  and  subtracting  one  from  both  sides 
yields 


»cIx.yl  ■  1  -  «2  'I  <»1*2>  «,1+i,.2+l 


(  H[x,y 1  )  ) 


Consider  first  jH[x,y].  Substituting  the  Mellin  transform  for  ^HJx.y] 
given  by  (3.18)  into  the  equation  above  and  replacing  s^  by  Sj+1  and 
S2  by  Sj+l  yields 


^[x.yl  -  1  -  M 


-1, 


-81 


~Sr 


®1®2S1S2 


x1(81+l)x2(82+2)x3(s1+S2+2)g1  g2  ] 


where  x^(u),  x2(v),  and  Xj(u+v)  are  defined  by  (3.7),  (3.8),  and  (3.9) 
respectively.  Recognizing  that 


^  r(sj)r(s2) 

sj72  "  r(s1+i)r(s2+i) 

Substituting  this  identity  into  the  inversion  above  yields 
xHJx.y]  -  1  - 

r(sj)r(s2)  -83  -82 

r(s1+l)r(s2+l)  xi<si+1>*2<82+l>x3<sl+s2+2>8l  g2 


Completing  the  inversion  and  using  the  bivariate  H-function  definition 
(3.6),  the  cumulative  H-function  form  may  be  found. 


where  K  -  k/gjg2 


(4.13) 


Using  property  (3.16),  m  ■  n  *  1,  equation  (4.13)  can  be  written  as 
^[x.y]  -  1  - 

(ej.E^ 

(aj*Aj)>(0,l)  i  (tj iCj ),(0,1) 
(fi*Fi) 

(-l,l),(b1,Bi)  ;  (-l,l)t(di,Di) 
where  K  ■  kxy  (4.14) 

Using  the  sane  procedure  outlined  above,  a  similar 


M^I.Nj.Mj+I.Nj^.Nj  glx 

^1^1 ’Ql^l ’ ^2^1 ’ ’ ^3 1 ^3 


representation  can  be  obtained  for  the  cunulative  distribution 
function  for  2H[x,y]. 


,H  fx.yj  -  1  - 


<VEi> 


glx  (ai+Ai,Ai),(l,l);(l,l),(ci-C1,Ci) 
tQj^l  •  P2^"l  »Q2"^1  ’  F3  *^3  (fi»^i) 

(0,l),(bi+Bi,Bi);(di-Di,  DjMO.l) 


where  K  -  k/g1g2 


(A. 15) 


Using  property  (3.17),  m  ■  n  ■  1,  equation  (4.15)  can  be  written  as 
2Hc(x,y]  -  1  - 

(ej.Ej) 

M1+l,N1,M2,N2+l,M3,N3  glx  ^ai*A1),(0, 1)  ;  (2, 1) , (Cj ,Cl ) 

2  P1+1,Q1+1,P2+1.Q2+1,P3.Q3  g2y  (fi*Fi) 

(-l,l),(bl,B1)  ;  (dj.DjMl.l) 


where  K  -  kxy 


(4.16) 


As  indicated  by  (4.14)  and  (4.16),  rHc(x,yJ  can  be  found  at 
the  same  time  as  rH(x,y]  by  using  the  calculations  for  the  residues  of 
rH[x,y].  For  jHfx.yJ,  multiply  each  residue  in  the  s2  plane  by 

l/(82-l),  and  then  add  the  pole  s2*l  (or  increasing  by  1  the  order  of 

an  existing  pole  at  s2*l).  For  2H[x,y],  multiply  each  residue  in  the 

s2  plane  by  l/(-l-s2),  and  then  add  the  pole  s2— 1  (or  increasing  by  1 

the  order  of  an  existing  pole  at  s2— 1).  For  both  1H[x,y)  and 
2H[x,y],  multiply  each  residue  in  the  s^  plane  by  l/(s^-l),  and  then 
add  the  pole  s^«l. 


Another  formula  for  the  cumulative  distribution  function  of  an 


H-function  probability  density  can  be  derived  using  the  Mellin 
transform. 


THEOREM  4.2;  The  bivariate  cumulative  distribution  function  for  a 
bivariate  H-function  probability  density  function  is  a  bivariate 
H-function. 


Extending  the  Mellin  transform  for  integrals,  ( 15 : 12 ; 106 :269) , 
to  two  variables  and  letting  -1/Sj  -  I*(-s,  )/r(l-Sj)  and  -l/s2  ■ 
r(-s2)/r(l“s2)  yields  the  bivariate  Mellin  transform  property 


x  y 

M  {/  /  f(u,v)dudv} 
ls2  0  0 


r<-s1)r(-s2) 

r(i-Sl)r(i-s2)  Msl+i,s2+itf(x*y)} 


(4.1 


This  property  may  also  be  obtained  as  a  simple  perturbation  of 
property  (2.14). 

Applying  (4.17)  and  (3.18)  to  jHfx.y]  yields 


xHcfx.y]  ■ 


kr(-s1)r(-s2) 

g1g2r(i-s1)r(i-s2) 


-si  -s2 

x1(s1+l)x2(s2+l)x3(s1+s2+2)g1  g2 


(e^+2Ej > ) 


KjH 


m1,n1+i,m2,n2+i,m3,n3 
P  L+1 ,  Q1-»-l ,  P2+l ,  Q2+l ,  P3 ,  Q3 


®lx 

g2y 


(l,l).(a1'*A1,A1);(ltl)t(c1+CltC1) 

(fi+2Fi‘Fi> 

(bi+B1,B1),(0,l);(di+D1,Di),(0,l) 


where  K  *  k/g^g2 


(4.i; 


^Hc[x,y]  may  also  be  obtained  by  reversing  the  order  of 
integration  of  the  cumulative  distribution  function  integral  and  the 
bivariate  H-function  contour  integral.  Let 


1  _S1  ~s2 
H[x,y]  -  k  - II  (~)(g.x)  (g2y)  dsjdsj 

(2ni)2 

where  //  represents  the  bivariate  H-function  contour  integrals  as 
defined  by  (3.6)  and  (-)  represents  the  gamma  products  in  the 
bivariate  H-function  definition  (3.6)  which  do  not  depend  on  the 
variables  x  and  y.  Assuming  convergence  of  the  H-function,  then 


x  y  i  ~si  ~s2 

,H  [x,y ]  -  /  /  k - ,  II  (~)(g,u)  (g2v)  ds.ds2dudv 

1  c  00  (2xi)^ 


i  x  -8.  y  -s2 

-  k -i — „  II  (-)[  /  (glu)  /  (g2v)  dvdu  ]  ds1ds2 


(2ni) 


-s 


1, 


-s  - 


,  (g,x)  (g2y) 

k  - -  fI  /!-.  W1  —  N~ 


(2xi) 


(l-s1)(l-s2) 


(4.1 


replacing  [ ( 1  — s ^ ( 1  — s 2 )  1  by  ^ <  1  — s ^ ) T ( 1  — s 2 > / T ( 2 — s  1 ) T ( 2 — s 2 )  and  using 

(3.6),  (4.19)  may  be  expressed  as 


iHc(x.y)  - 

M1,N1+l,M2,N2+l,M3iN3 


KtH 


Px+1 ,Q1+1 , P2+l , Q2+l , P3 , Q3 


gxx 

g2y 


(ej.Ej) 

(Otl)t(ai,A1)  ;  (0,1), (c^C*) 
<fi’Fi> 

;  (di,D1),(-l,l) 


where  K  «  kxy 


(4.20) 


Using  property  (3.16),  m  *  n  ■  1,  (4.20)  may  be  expressed  as  (4.18). 

Using  similar  procedures,  the  cumulative  distribution  function 
2Hc[x,y]  for  2H[x,y]  may  be  found  and  is  also  given  as  a  bivariate 
H-function.  Using  property  (4.17)  and  the  Mellin  transform  for 
2Ht*»y]»  (3*20),  2Hc[x,y]  is  given  as 


2Hc[*»y]  ■ 


kr(-s1)r(-s2) 

8l82r<1“sl)r(1"s2) 


~S1  -s2 

x1(s1+l)x2(-s2-l)x3(s1-s2)g1  g2 


Replacing  s2  by  -s2  and  completing  the  inversion  by  use  of  (3.12),  the 
cumulative  distribution  function  for  2H[x,y]  may  be  found  and  is  given 


2Hc[*.y] 


(e^) 


k2h 


m1,n1+i,m2+i,n2,m3,n3 
P^+1 , Q^+l t P2+i , q2+i ,p3,q3 


«lx 

g2y 


(l,l),(a1+A1,Ai);(ci-Ci,Ci),(l,l) 
(fi,Fi > 

(b1+B1,Bi),(0,l);(0,l),(di-Di,D1) 


where  K 


k/gig2 


(4. 2 


Using  property  (3.17),  (4.21)  may  be  expressed  as 


2Hc[x.y]  - 


P^+l ,Q^+1 ,  P2+l  ,Q2+1 ,  P3 ,  Q3 


8l* 

g2y 


^ei ,Fi ) 

(0,l),(ai,Ai)  ;  (ci,Ci),(2,l) 

(fi.Fi) 

(bi,Bl),(-l,l)  ;  (l,l),(di,Di) 


where  K  »  kxy 


(4.2 


Bivariate  Gamma  Cumulative  Distribution:  Applying  (4.18)  to  the 
bivariate  gamma  probability  density  function  given  by  (4.3), 


iHc[x,y] 


ap+q 

r(p)T(q) 


x  y 

;  / 
n  n 


,P_1 


(v-u)**”*  e”av  dudv 


( p+<? .  i ) 


l 

r<p)  i 


1,1,0, 1,1,0 

H 

2, 2, 1,1, 0,1 


ax 

ay 


(l,l),(p+q,l)  ;  (1,1) 


(p, 1) ,(0,1)  ;  (0,1) 


y>x>0  ,  a,p,q>0 


(4.23 


Bivariate  Beta  Cumulative  Distribution:  Applying  (4.18)  to  the 
bivariate  beta  probability  density  function  given  by  (4.4), 


iHc[x.y] 


r(P!+P2+P3)  rx/  Pi*1  P2*1 

r(p1>r(p2)r(p3)  ;o;o  u 


P3-1 

(1-u-v)  dudv 


r(Px+P2+P3) 

r(p1)r(p2) 


1,1, 1,1, 0,0  X 

1H 

1,2, 1,2, 1,0  y 


(1,1)  ;  (1,1) 

(Pl+P2+P3»l) 

(P1,1),(0,1)  ;  (p2,l),(0,l) 


x,y>0,  x+y<l,  PL,P2,P3>0  (4.24 


Quarter  Cauchy  Cumulative  Distribution:  Applying  (4.18)  to  the 
quarter  cauchy  probability  density  function  given  by  (4.5), 


iH„lx.y]  “ 


dudv 


Hc[x,y] 


6a2  /  /  e‘(au  +  Sv/u) 
o  o 


■  H 


0,1, 1,1, 1,0 
1,1, 1,2,0, 1 


ax 

aSy 


(2,1) 

(l.D  ;  (i,i) 


(0,1)  ;  (1,1), (0,1) 


x,y>0,  a,6>0 


(4. 


ellogg-Barnes  III  Cumulative  Distribution:  Applying  (4.21)  to  the 
Kellogg -Barnes  III  probability  density  function  given  by  (4.8), 


1HC [x»y 1 


/V  uc  e"<au  +  Buv>  dudv 

r(c)  o  0 


0,1, 1,1, 1,0 


r(c)  2  1,1,2, 1,0,1 


ax 

iy 

a 


(c,l) 

(1,1)  ;  (0,1), (1,1) 

(0,1)  ;  (0,1) 


x,y>0,  a,$>0,  c>2 


(*■ 


4.6  Evaluation  of  the  Bivariate  H-function  Constant 


Carter  (3),  Eldred  (7),  and  Springer  (18)  provide  special 


cases  and  a  definition  of  the  univariate  H-function  distribution . 
However,  it  was  Cook  (5)  who  gave  a  methodology  for  determining  the 
normalization  constant  k  in  definition  (4.1).  A  similar  procedure  can 
be  applied  to  the  bivariate  H-function  distribution  constant. 

One  approach  to  finding  k  is  to  investigate  rHc[x,y]  for  large 
x  and  y,  since 

lim  rHc(x,y]  *  lim  k(H-function  given  by  (4.18)  or  (4.21))  *  1. 

x-»  ®  x-*°* 

y  — *  oo  y— »oo 

That  is,  if  k^H^x.y]  is  a  proposed  H-function  probability  density, 
use  (4.18)  or  (4.21)  to  find  the  associated  H-function  for  the 
cumulative  distribution  function,  which  for  large  x  and  y  will  be  1/k. 

While  the  above  method  is  feasible,  it  yields  a  numerical 
approximation.  An  exact  method  may  be  found  by  equating  the  right 
sides  of  (4.14)  and  (4.20),  for  ^Hc[x»y],  which  immediately  yields: 


(ei ) 


m1+i,n1,m2+i,n2,m3,n3 

P^l  ,Q1+1  ,P2+1  ,Q2+1 ,  P3  ,Q3 


®lx 

g2y 


+  1H 


m1,n1+i,m2,n2+i,m3,n3 
Pj^+1  ,Q^+1 ,  P2+l  »Q2+1 » ^3  *^3 


8lx 

g2y 


(ai *At ) * (0, 1 )  ;  (e1,C1)1(0>l) 
<fi’Fi> 

(-1,1), (bi,B.)  ;  (-l,l),(d1,Dl) 

^ei,E1) 

(0,l),(ai,Ai)  ;  (0,1), (cj,^) 
<fi’Fi> 

(VBiM-i.i)  ;  (d1,D1),(-i,i) 


;/ 


(2nl )' 


(*-)  + 


(S1-1)(S2-1)  '  (l-81)(l-82) 


(-) 


ds1ds2  (4.29) 


where  SS  represents  the  appropriate  contour  integrals  as  defined  in 

“S1  ~s2 

(3.6)  and  (-)  *  x^Sj )x2(s2)x3(s j+s2)(gjx)  (g2y)  also  defined  in 

(3.6) . 

Compare  the  residues  of  the  two  bivariate  H-functions  in 
(4.29).  Inverting  with  respect  to  s2,  each  RHP  residue  of  the  first 
bivariate  H-function  has  a  matching  RHP  residue  of  the  second 
bivariate  H-function  that  is  exactly  equal  but  opposite  in  sign, 
except  at  the  residue  at  ■  1.  Similarly,  each  LHP  residue  of  che 
second  has  a  matching  residue  of  the  first  that  is  exactly  equal  but 
opposite  in  sign,  except  at  the  residue  at  s2  -  1.  Similar  results 
hold  true  for  inverting  with  respect  to  s^.  Therefore,  whether  (4.29) 


is  inverted  by  LHPCs^),  LHP(s2);  LHP(s^),  RHPCsj);  RHP(s^),  LHP(S2>; 
or  RHP(s^),  RHP(s2).  it  reduces  to  only  one  term  on  the  right  side: 

—  »  (+LHP  LHP  residue  at  s,«s,«l) 

Kxy  S.  L  L 


(-LHP_  RHP„ 

S.  S  M 

residue 

at 

S^Sj-l) 

1  2 

(-RHP  LHP 

S1  ®2 

residue 

at 

s1-s2-l) 

(+RHP  RHP 

residue 

at 

8^82-1) 

If  the  bivariate  probability  density  function  jHfx.y]  has  no 
pole  at  Sj  *  1  or  »  1(  then  the  bivariate  cumulative  distribution 
function  jH.fx.y]  has  a  pole  of  order  1  at  Sj  »  1  and  s2  *  1  and 
(4.29)  reduces  to: 


M1  N1  M2 

.  n  r(b.+B. )  n  r(i-a,-A. )  n  r(d,+D. ) 

L  *  1  l  l 


^1  ^1  ^2 

n  r (a .+A. )  n  r(i-b,-B.)  n  r^+c,) 

Nl+1  Mx+1  N2+l 


N,  M.  N, 

n  r<i-c  -c  >  a  r(e.+2E. )  n  r(i-f.-2F. ) 

x  1 _  1  1 _ _1 _  (4.3 

^  giX82y 

ir  r(i-d  -d.)  n  r(f.+2F.)  n  r(i-ei-2E.) 

m2+i  1  1  n3+i  m3+i 


Solving  for  k  and  noting  the  Mellin  transform  property  (3.18) 


%  . 


IV 


V.\ 

v 


»S;V 

R$ 


m 


i.-> 


m 


k: 


h\ 


K;. 


r-  .  - 
v> 


/,  - 
>, '. 
y. 


:<• 

o»\- 


K'i 


V 
%  w  ■ 


k  “  (  1/MS  c  {  i  >  ) 

sl’s2  1 


Si-1 

s2-i 


The  above  result  is  summarized  in  the  following  theorem. 
Theorem  4.3:  If  ^H(x,y]  is  a  bivariate  H-function  probability  density 
defined  by  (4.2)  for  r  -  1,  such  that  — <  1,  i  ■  1,...,M^, 

-d1/Di  <  1,  i  -  1,...,M2,  (ai+Ai)  <  1,  i  -  1 , . . . ,N^ ,  (Cj+Cj)  <  1, 
i  »  1,...,N2,  (which  implies  that  ^H[x,y]  has  no  pole  at  s^  *  1, 
s2  -  1),  then 


-  (  1/M  {  jH[x,y]  }  ) 

S1  ,s2 


V1 


*2*1 


gl*2 


1  ^1  P2 
n  r (a.+A. )  n  r(i-b1-B.)  n  iXc.+c^) 

mi+i  1  ”  •* 


nl+i 


N2+l 


M.  N,  M2 

n  r (b.+B. )  a  r(i-a.-A. )  n  r(d.+D.) 
1  1  1  1  1 


Q3 


a  r(i-d1-Di)  n  r(f1+2F1)  n  r(i-ei-2Ei) 


m2+i 


n3+i 


m3+i 


(4.31) 


n2  m3  n3 

n  r(i-c  -c  )  n  r(e.+2E, )  n  r(i-f.-2F. ) 
1  1  1  1 


Example  4.2 :  Consider  the  probability  density  function  given  by 


/  x  i  ~ax  -  By/x 

7v  VI  m  V  A  y 


V  .  v>0 
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k  1 


H 


0,0, 1,0, 1,0 
0,0, 0,1, 0,1 


ax 

aBy 


from  Theorem  4.3,  k  is  given  as 


k  -  (a)(aB) 


r(l)T(2) 


a28 


(0,1) 


;  (0,1) 


which  agrees  with  the  definition  of  the  Kellogg-Barnes  II  probability 
density  function. 

The  five  special  cases  of  the  bivariate  H-function  probability 
density  functions  given  in  section  4.3  all  meet  the  conditions  of 
Theorem  4.3  and  their  constants  agree  with  (4.31) 

The  normalizing  constant  for  jHfx.y]  may  be  found  in  the  same 
fashion  by  equating  the  right  sides  of  (4.16)  and  (4.22), 


IV 

tvfl 


r. 

*  v  s 
AV.V', 


Mj+I.Nj^.Nj+I^.Nj  glx  (a^AjMO.l)  ;  (2,  IMCj.Cj) 

i 

P1+1,Q1+1,P2+1,Q2+1,P3iQ3  g2y  (fi»Fi> 

(-l,l),(bi,Bi)  ;  (dj.DjMl.l) 


(ej.Ej) 

gXX 

(0,l),(ai,A1)  ;  (c 

g2y 

(fi’Fi> 

*^1+1 ,M2+1 , N2 tM3 , N3  g^x  ( 0, 1 ) , (a^ , )  ;  (c^,Cj), (2,1) 

Pl+1  »Qi"*'l  1  P2+l  ,Q2+1  tP^  i  Q3  (f^»Pj) 

;  (l,l),(di,Di) 

■  —  0  11  ,  W  , - x  <~>  +  ,T  .  wix-  \  dsids,  (4.32) 

(  2t»  i ) 2  L(81"1)(_1*82>  (1-81)(1+s2)  J  1  2 


where  //  represents  the  appropriate  contour  integrals  as  defined  in 

~8l  s2 

(3.6)  and  (~)  -  x1(s1)x2(s2)x3(s1+s2)(g1x)  (g2y)  also  defined  in 

(3.6) . 

Compare  the  residues  of  the  two  H-functions  in  (4.32). 
Inverting  with  respect  to  s2,  each  RHP  residue  of  the  second 
H-functlon  has  a  matching  RHP  residue  of  the  first  H-function  that  is 
exactly  equal  but  opposite  in  sign,  except  at  the  residue  s2  »  -1. 
Similarly,  each  LHP  residue  of  the  first  has  a  matching  residue  of  the 
second  that  is  exactly  equal  but  opposite  in  sign,  except  at  the 
residue  at  s2  •  -1.  Similar  results  hold  true  for  inverting  with 
respect  to  s^  except  that  extra  pole  is  at  s^  ■  1. 
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If  the  probability  density  function  2H[x,y]  has  no  pole  at 
■  1  and  S2  ■  -1,  then  the  cumulative  distribution  function  2Hctx*y] 
has  a  pole  of  order  1  at  Sj  ■  1  and  s2  ■  -1  and  (4.32)  reduces  to 


M1  N1  M2 

l  n  r(b1+B1)  n  ra-a^Ap  n  r^-op 

kxy  *  J - 5 - j - 

n1  r (a .+a. )  n1  r(i-b,-B, )  n2  r(C.-c. ) 

Nj+l  1  Mj+1  N2+l 


N  M  N3 

n  r<i-c  +c  )  n  r(e.)  n  r<i-f  > 

l  ill 


n  r<i-d  +d .)  n  r(f.)  it  r(i-e.) 

•  1  ^  ^  kT  l  4  ^  X  .  4  ^ 


g!xg2y 


m2+i 


n3+i 


m3+i 


(4.33) 


solving  for  k  and  noting  the  Mellin  transform  property  (3.21)  yields 


(  1/Ms  s  {  2H*x’yJ  }  > 
Sl‘s2  l 


•lm\ 

*2-1 


The  above  result  is  summarized  in  the  following  theorem. 
Theorem  4.4:  If  2H(x,y]  is  a  bivariate  H-function  probability  density 
defined  by  (4.2)  for  r  -  2,  such  that  -bj/Bj^  <  1,  i  ■  1,...,!^, 
d1/Di  >  l,  i  -  1,...,M2,  (a1+Ai)  <  1,  i  -  1,...,N1,  (cj-Cj)  <  1, 
i  ■  1,...,N2,  (which  implies  that  2Hlx»y]  bas  no  pole  at  ■  1 
or  s2  «  -1),  then 


k  "  (  1/MS  s  {  2Htx-yi  }  > 

Sl*s2  L  s.-l 


P1  Q1  P2 

n  r(a1+A1)  n  r(i-b1-B1)  n  r^-Cj) 


m,  n,  m2 

n  r(b  +b  >  n  r(i-a,-A. )  n  r(d1-D, ) 
1  i  1 


Q2  p3  Q3 

n  r<i-d  +d.>  n  r(f.)  n  ro-e.) 

<2+1  N^+l  M^+l 

N,  M~  N3 

n  r(i-c  +c.)  n  r(e.)  it  ro-f  ) 

i  i  1  l 


(4.34) 


Example  4.3:  Consider  the  bivariate  probability  density  function 


given  by 


fXfY(x,y)  -  k(ax)c  e-(“x  + 


x,y>0 


(c,l) 


0,0, 0,1, 1,0  ox  - ;  (1,1) 

k  2H  S 

0,0, 1,0, 0,1  £  y  - 

a 


Conditions  of  Theorem  4.4  are  satisfied,  and  k  is  given  by 
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t* 


k  •  <!>"”  mki  • 


which  agrees  with  constant  for  the  definition  of  the  Kellogg-Barnes  II 
distribution. 

Another  way  to  arrive  at  Theorems  (4.3)  and  (4.4)  is  to 

consider  the  zero-zero  moment  of  a  probability  density  function 

k..H[x,y],  x,y  >  0,  and  M  {  _H[x,y]  }  has  no  pole  at  s,  »  1  or 
r  s  *  S2  r 

S2  *  1,  for  r  »  1,  or  at  Sj  ■  1  or  Sj  *  -1,  for  r  *  2: 

00  OO 

E[x°y°]  -  /  /  kH[x,y]  dxdy  -  1 
0  0 


“  Ms  o  t  kH[x,y]  } 
S1  ,s2 


s.-04-l 

S2*0+l 


V 


§ 

$ 

r-y- 

'• 


tv.' 

.V.-. 

,v. 

.v.* 

I 


>y 


k  Ms  .  {  H(x,y}  } 
l’s2 


sr; 

s2-l 


solving  for  k  yields 


k  -  (  1/M  {  H[x,y ]  }  ) 

l,s2 


v; 

82-l 


from  which  Theorems  4.3  and  4.4  immediately  follow. 


CHAPTER  5 


Transformations  of  H-function  Variates 
5. 1  General  Remarks 

This  chapter  has  two  distinct  parts.  In  the  first  part, 
various  combinations  of  products,  quotients,  and  powers  of  dependent 
H-function  variates  are  examined.  Theorems  are  presented  to  show  that 
the  product  or  quotient  of  two  dependent  H-function  variates  is  an 
H-function  variate.  Powers  of  dependent  H-function  variates  also 
result  in  H-function  variates.  This  result  has  significance  since  if 
the  joint  density  function  of  two  dependent  variates  can  be 
represented  as  a  bivariate  H-function,  then  the  probability  density 
function  of  a  random  variable  of  the  form  Z  ■  X1^  is  given  by  a 
univariate  H-function  and  may  be  inverted  using  the  inversion  routine 
given  by  Cook  (3). 

The  second  part  of  the  chapter  examines  various  combinations 
of  H-function  variates  from  two  or  more  bivariate  Il-function 
probability  density  functions.  Theorems  are  given  for  products  and 
quotients  of  variates  from  bivariate  H-function  distributions  which 
are  pairwise  independent.  These  results  are  also  significant  in  that 
by  combining  these  results  with  the  results  from  part  one,  a  powerful 
general  theory  for  finding  rational  combinations  of  mixtures  of 
dependent  and  independent  random  variables  now  exists. 
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5.2  Transformations  of  Dependent  H-function  Variates 


5.2.1  Transformations  of  ^H[x,y ]  Variates 

As  is  stated  in  the  following  theorem,  one  of  the  most 
significant  properties  of  the  bivariate  H-function  distribution  is 
that  the  probability  distribution  of  products  or  quotients  of  its 
variates  reduces  to  a  univariate  H-function.  This  is  similar  to 
results  derived  by  Carter  (3)  for  independent  H-function  variates,  but 
this  theorem  shows  the  result  holds  for  dependent  H-function  variates 
as  well.  It  is  well  known  that  such  a  property  is  not  common  among 
the  classical  distributions. 

Theorem  5.1:  If  X  and  Y  are  H-function  variates  with  joint 
probability  density  function  fx  ^^.y)  where 

f  kiHl8ix»g2y)  *  *.y  >  0 

f x  v^’^  - 

’  [_  0  ,  otherwise 

then  the  probability  density  function  of  the  random  variable 

Z-  XpY^ 

p,q  rational,  is  given  by 


fZ<*> 


M  N 
H 

P  Q 


2  P2  <*Z 
81  *2  z 


'2  -I 


,  z>0 


(5.1) 


0 


otherwise 


where  for 


Case  I_  p>0  ,  q>0 

K-k,  M-ZM1,  N-IN.,  P-ZPj,  Q*ZQ1 ,  1=1, 2, 3  , 

*  (a^+Aj-A^p,A^p) ,  1*1. • «Np  (c^+C^-C^q ,Cjq ) ,  i—l...^, 

(f1+2Fi-Fi(p+q),Fi(p+q)),  1-1... Nj, 

p  ,  A^  p)  |  1  *N  (c  ^■fC^  -C^q  , C ^q  )  ,  i  =N2+  1 .  • . 

(fi+2Fi-Fi(p+q),Fi(p+q)),  i-N3+l...P3 
n2  *  (bi+Bi_BiP»BiP)»  i*l***M1,  (di+Di-Diq,Diq), 
(e1+2Ei-Ei(p+q),Ei(p+q)),  1-1. ..M-j, 

( "B^p, B^ p)i  1  *M ^+1 . .  »Qp  ( “D^q ,  D^q  )  ,  1  =M2^'  1 .  • . 
(e1+2Ej.-Ei(p+q),E1(p-hi)),  1-M3+1..  .Q3 
Case  II  p>0  ,  q<0  ,  |p|>jqj 

K-k,  N-N1+M2+N3,  P-P1+Q2+P3,  Q*Q1+P2-H33, 

*  (a^+A^-A^p,A^p),  1*1 « • .Np  ( 1  — d^ —D^+D^q ,  D^q  ) ,  1*1. ..M2 
(fi+2Fi-F1(p+q),Fi(p+q)),  1-1. ..N3 

(ai+Ai-AiP»AiP)»  l-N^+l...Pj,  (  l-d^-D^+D^q ,  -D^q  ) ,  i-M2+l 
(  f  i+2Fi  "Fi  <  P^  >  *  Fi  <P-^  >  >  •  i-»3+l  •  •  •  f3 


n2  *  (b^+B^~B^p ,  B^p ) ,  1*1. ..Mj,  ( l—c^“ C^+C^q » ~C^q  )»  i—l*...^ 
(e ^+2Ej -E^ (p+q ) ( E ^ ( p+q )),  1*1... M3, 

(b^+B^-B^p , B^p ) ,  i-M^+1. . .Q^,  ( 1  — c ^-C^+C^q , “C^q ) ,  1—N2+1 


(ei+2Ei-Ei(p+q),Ei(p+q)),  l-M3+l...Q3, 

Case  III  p>0  ,  q<0  ,  |p|<|qj 

K-k,  M-M1+N2+N3,  N-Nl4M24M3,  P-P. +Q2+Q3 ,  Q-Q.+Pj+P^, 


0^  *  (a^+A^-A^p.A^p) ,  iBl#  « »Nj,  (1“^j“Dj+D^(|i*D^(|),  iml*##M2 

(l-ei-2E1+Ei(p+q),-E1(p+q)),  i-1. •  .M-j 

^al+Ai_Aip,Aip^’  i“Ni+1,**pi»  (1“di_Di+Di<l»~D£cl)»  i"M2+l...Q2 
(l-e1-2Ei+E1(p+q),-E1(p+q)),  i-M-j+1..  .Q3 

^2  *  i-l...Mj,  ( 1  ~c ^ “C^+C^q , “C^q ) ,  i— l...N2, 

C l“f 1“2Fi+F1(p+q ) , -F1 (p4q ) ) ,  i-l...N3 

^bi+Bi“Bip,Bip^’  ^  1"ci"Ci+Cicl»“Citl)»  i*N2+l...P2 

(1'fi“2Fi+Fi(p+<l)*~Fi<p+q))*  i"N3+1*,'P3 
Case  IV  p>0  ,  q<0  ,  |p|»|q| 


K  -  k 


m3  n3 

n  r(e  +2E.)  n  rci-f.^Fj) 
1  1 _ 

F3  ^3 

n  r(f1+2F1)  n  rci-e^Ej) 


M-M1+N2,  N-N^,  P-P^,  Q-Q1+P2, 

*  (a^+Aj*A^P)Ajp) |  i-1. • «N^,  (l“d^—D^+D^q,~D£q),  i*l«*»M2> 
(ai+Ai-Aip,Aip),  i-N^+1.. .  P^,  (l-dj-Dj+D^.-Djq),  i-M2+l., 
^2  *  (b^+Bj~B^p | B^p ) ,  i— l...N2, 

(b  j+Bj  ~ B^p ,  Bjp  ) ,  i-M^+l. .  .Q^,  (l“Cj“C^+Cjqi*Cjq)i  i*N2+l»i 
Case  p<0  ,  q<0 

K-k,  M-ZN, ,  N-ZM. ,  P-ZQ, .  Q-ZP, ,  i-1, 2, 3  , 


«  (l-bj-Bj+BjP.-Bjp),  (l-di-Di+Diq,-Dtq),  i-l...M2, 


(l-e^E^E^p-KO.-E^p+q)),  1-1... H3 
(l-b1-B1+B1p,*Bip),  i-Mj+1. . .Q1# 

( l~d D^+D^q , — D^q ) ,  i“M2+l...Q2, 
(l-e1-2Ei+Ei(p+q),-Ei(p+q)),  i«M3+l...Q3 

n2  *  ( l~a^—A^+A^p t *A^p ) y  i**l . . . N| t  ( 1  “C  ^ “C^+C^q , —C^q ) t  i*l«»*N2# 

(1_fi“2Fi+Fi(p'Hl)'"Fl(p+q))*  1*1-*‘N3 
( 1"3^**A^+A^p  1  “A^p  )  |  1  *N  ^1  •  •  •  P  3 1 

(l^c1-Ci+Ciq,-Ciq),  i-N2+l...P2, 

( 1-f i-2Fi+Fi (p+q ) , -Fi (p4q ) ) ,  i-N3+l . . . P3 


Proof  for  Theorem  5.1  From  Theorem  2.4 


-1, 


fz(z)  -  Mj  [  Mps-p+i,qs-q+l*  fx,Y(x,y)  }  1 


-1 


M1  1  Mps-p+l,qs-q+l{  kiHtgl*’g2yl  }  * 


and  from  Equation  (3.19) 
-1 


f7(z)  »  M, 


k-g  1  ”*Ug 2  Vx1(u)x2(v)*3(u+v) 


u*ps-p+l 
v-qs-q+1  J 


h+i- 


fz(*)  -  k8lP”lg2q~1(2nl)  x1(u)x2(v)x3(u4v) 


(8lP82q*)'S  ds  (5.2) 


U-pS-p+1 

v-qs-q+l 


from  the  univariate  H-function  probability  density  function 


*  •  m.\  s-a  1  *l- 


definition,  (4.1),  and  the  univariate  H-function  definition,  (3.1) 
cases  1,  II,  III,  IV,  and  V  immediately  follow. 


Special  Cases 

Two  special  cases  of  this  theorem  that  are  of  particular 
Interest  are  a  simple  product  and  quotient  of  two  dependent  variables. 
The  distribution  of  Z  «  XY:  From  Theorem  5.1,  p  *  q  “  1,  (5.1) 
reduces  to 


fz(z) 


M  N 


kH 


P  Q  L 


®1®2J 


0  otherwise 


,  z>0 


where 

M-EM^  N 

-ZNlf  P-ZPj,  Q-EQj , 

i-i 

,2,3  , 

nl  ” 

(aj ) » 

i*l • • |  ) » 

i-i 

•  •  •  n2  1 

9  • • • N^|  | ) y 

i- 

Mj+1.. .P1 

1 

(ci »^i ) , 

i-N2+l...P2>  (f ± , 2F 

i>« 

i»N3+l.. 

*P3 

n2  * 

(bj.Bj), 

i-l 

•  •  •  M2  y 

(e1,2Ei) 

9  i*l • • *M3> (b^ » Bj ) 9 

i-M 

j^+1 ...  , 

(^i »^i ) » 

i*M2+l. • *Q2>  ) 2E 

i), 

i-M3+l..< 

'Q3 

(5.3) 


Equation  (5.3)  above  may  also  be  found  by  applying  Equation 
(2.23)  to  the  Mellin  transform,  (3.19),  of  the  bivariate  H-function 
definition.  Recombining  like  terms  and  inverting  by  (3.6),  Equation 
(5.3)  may  be  obtained  directly. 


» 
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(5.4)  may  be  obtained  directly. 

For  the  case  (5.4)  reduces  to  the  univariate 

H-functlon  probability  density  function  which  is  the  result  of  the 
ratio  of  two  independent  H-function  variates.  For  M^»N2*P2«Q3»0, 

(5.4)  agrees  with  the  form  given  by  Carter,  (3:62). 

Example  5.1:  Consider  the  Kellogg-Barnes  II  distribution  given  by 


fX,Y(x*y) 


exp( -ax  -  By/x)  , 


x,y>0 

a,B>0 


0,0,1, 0,1,0 

-  Sa  ,H 

0,0,0, 1,0,1 


ax 

aBy 


(0,1) 


;  (0,1) 


The  distribution  of  the  random  variable  Z®X/Y  may  be  found  with 
Theorem  5.1,  p  ■  1,  q  ■  -1,  or  from  (5.4)  directly,  and  is  given  by 


fz(z) 


Ba2 

77 


0  1  !  . 

r(2)H  ~  2 

i  o  L 8 


(-1,1) 


l  l_ 
B  2-ni 


/  r(2-s)(|  z)"8  ds 

h-i®  e 


(5.5 


Using  Cook's  convergence  theorem  for  the  H-function,  equation 

(5.5)  may  be  evaluated  using  the  sum  of  residues  in  the  right  half 
plane.  Poles  of  the  integrand  from  T(2-s)  occur  at  s®2,3,4, 


•  •  • 


Evaluating  the  residues  at  these  points  and  summing  yields 


r  residues  -  "  (-d'Ws)1 

s2  i-0  1! 


Recognizing  this  as  the  power  series  for  the  exponential,  the 
distribution  of  fj,(z)  is  then  given  by 

f_(z)  ■  4  ^  residues 
L  B 


8s*2 


which  agrees  with  the  results  obtained  using  the  Mellin  convolution 
integral  in  Example  1.5. 

Example  5.2:  Suppose  the  distribution  for  the  random  variable  Z»Y/X 
is  desired  for  X  and  Y  distributed  by  the  Kellogg-Barnes  II 
distribution  given  above.  Clearly,  Theorem  5.1  won't  hold  in  this 
case  since  it  does  not  allow  for  cases  where  p<0  and  q>0.  However,  if 
the  Mellin  transform  is  used  to  map  x  into  &2  and  y  into  s^  and  an 
appropriate  H-function  identity  found,  then  Theorem  5.1,  case  iv, 
holds.  The  H-function  representation  for  the  Kellogg-Barnes  II 
distribution  with  such  a  mapping  is  given  by 


Applying  Theorem  5.1,  case  iv,  the  distribution  for  the  random 
variable  Z-Y/X  is 


f7(z)  r<2)  H 


1  0 


0  1 


Bz 


(0,1)  J 


h+i— 

/  r(s)(6z)‘s  ds 

h-i« 


(5.6) 


The  distribution  f^(z)  can  be  found  by  evaluating  and  summing 
the  residues  of  (5.6)  in  the  left  half  plane.  The  integrand  of  (5.6) 
has  poles  at  s*0,-l,-2, . . .  Evaluating  and  summing  these  residues 
yields 


Z  residues  -  Z  (-l/O**)1 
1-0  1! 

Realizing  that  this  is  the  power  series  for  an  exponential,  the 
distribution  fjj(z)  may  be  given  as 

fg(z)  -62  residues 

-  (Se-Sz 

The  distribution  f^(z)  may  also  be  determined  directly  from  the 
H-function  form  of  Equation  (5.6).  Using  the  special  cases  given  by 


Carter  (3),  the  exponential  distribution  is  given  directly  by  (5.6). 
The  solution  is  verified  by  use  of  the  Mellin  convolution  integral  in 
Example  1.5. 


Example  5.3:  Consider  the  probability  density  function  given  by 


fX,Y(x*y)  “ 


1+8/2 


xS  fi-a(x2  +  y2) 


x,y  >  o 

a, 8  >  0 


Viewing  the  density  function  above  it  should  be  realized  that 
the  function  can  be  expressed  in  the  form  fx(x)*fy(y)»  Restated,  X 
and  Y  are  Independent.  While  theorems  given  by  Carter  (3)  may  be  used 
to  find  the  distribution  of  Z  -  XY,  it  is  illustrative  to  treat  the 
density  function  above  as  a  joint  density  and  use  Theorem  5.1  to  find 
the  density  function  for  the  random  variable  Z  »  XY. 

Applying  (2.5),  followed  by  (3.6),  to  the  joint  probability 
density  function  above,  the  H-function  representation  may  be  found  and 
is  given  by 


From  Theorem  5.1,  case  1,  p  ■  q  ■  1,  the  distribution  of  Z  ■  XY  is 
given  by 


—  f  e"(u  +  (az/u>  >  du 

/¥  o 


which  agrees  with  the  solution  above  using  a  bivariate  H-function 
probability  density  function  and  Theorem  5.1. 

5.2.2  Transformations  of  ^Hfx.y]  Variates 

Theorem  5.1  shows  that  the  probability  density  function  of 
products  or  quotients  of  jH[x»y]  variates  is  a  univariate  H-function. 
The  following  theorem  shows  this  property  is  shared  by  2H[x,y) 
variates . 

Theorem  5.2:  If  X  and  Y  are  H-function  variates  with  joint 
probability  density  function  f^  y(x,y)  where 


fx,y(x,y) 


k2Hl8ix»82yJ  *  x»y  >  0 

0  ,  otherwise 


then  the  probability  density  function  of  the  random  variable 


z-  yptf 


p,q  rational,  is  given  by 


f7(z) 


,  ,  M  N  .  n. 

Kg,15'1  gS~l  H  g.PgA 


otherwise 


,  z>0 


(5.9) 


where  for 

Case  1^  p>0  ,  q>0,  p>q 

JOk,  N-N14M2+N3,  P-PJ+Q2+P3,  Q-Q3+P2-K33 

nl  *  (ai+Ai“AiP‘AiP)»  i"l»»*M2, 

(f j-F^p+F^q.F^p-F^q),  1*1. ..N3, 

( a i +Af -AipjAfp),  i«Nj+l...P^,  ( l"d j+D^-D^q ,D^q ) ,  i **M2+1 • • • Q2 » 

( f^-Fjp+Fj.q,F^p-F^q  ) ,  i*N3+l . .  .P3 

n2  -  (bj+Bj^p.Bjp),  i-l...Mlf  (l^ci+Ci”C1q»Ciq),  i-l...N2« 
*ei_Elp+Eiq,Eip“Eiq^*  i-1,“M3» 

(b1+B1-Bip,Bip)>  1-Ml+1...Q1,  (l-Ci+Cj-Cjq.Cjq),  1-N2+1...P2, 
(ej-Ejp+E^q.Ejp-Ejq),  i«M3+l. .  .Q3 
Case  II  p>0,  q>0,  p<q 

K-k,  M-M1+N2+N3,  N-N14M2+M3,  P-P1-KJ2+Q3*  Q“Q1+P2+P3 
nl  *  (ai+Ai“AiP,Aip),  i*l»  ♦  »Nj,  (l-dj+Dj-Dj^q.Djq),  i»l...M2, 

(l-ej+E^p-E^-EjP+Ejq),  i»l..  .M3, 

(a^+A^-A^pjA^p) ,  i“N^+l. .  .P^,  (l”d£+Dj-Djq,Djq),  i*M2+l...Q2, 
( 1-ej^+Ejp-Ejq ,  -EjP+Ejq ) ,  i-Mj+l . .  .Q3 
h2  *  (b^+Bj ~B^p , B^p ) ,  ( 1-c £+C^—Cjq ,C^q ) ,  i*l...N2, 

(l-f1+F1p-F1q,-F1p+Fiq),  i-I.. .N3, 

(b^-fBj_B^p,  B^p) ,  i**Mj+l. « .Qj,  (l^c^+Cj^C^q.C^q),  I*N2+1...P2, 
(l-f1+Fip-Fiq,-Fip+F1q),  I-N3+I.. .P3 


1 


i 


\  \  V  \\V*  V 


Case  III  p>0,  q>0,  p*q 


n3  r<e  >  n3  ra-^) 

K,k| - i - 

ii3  r ( f  >  n3  r(i-ei) 

n3+i  m3+i 


m*m1+n2,  n=n1+m2,  p-p1-k)2,  q^q^Pj 

ni  "  (a^+A^-A^p ,A^p) ,  (l-dj+D^^-Djq.Dj^q),  i"*l...M2, 

^ai+Al_AiP’AiP^*  i*Nj+l. . •Pj,  (l'dj+Dj-Djq.Dj^q),  i=M2+l...Q2 
n2  -  (bi+Bi~Bip,B1p),  (l-c^+Cj-Cjq.Cjq),  1-1. ..N2, 

(bi+Bi_Bip,Bip)’  i“Mi+1***Ql»  (l-c1+ci-ci‘l»ci<l)»  i-N2+l...P2 
Case  IV  p>0,  q<0 

K-k,  M»EM1 ,  N*ENj ,  P-ZPj,  Q-EQj ,  1*1,2, 3 

nl  *  (aj+Ai~Aip,A1p),  1*1. .«Nj,  (Cj-Cj+Cjq.-C^q),  1*1. ..N2, 
(f^-FjP+F^q,FjP_Fjq ),  1*1. ..Nj, 

(ai+Ai“AiP»AiP)»  l*^i+l«  *  >Pj,  (Cj-Cj+Cjq.-Cjq),  i“N2+l...P2, 
^fi~Fip+Fiq»Fip-Fiq),  i*N3+l. . .p3 
n2  -  (bi+Bi-Bip,Bip),  1-1. ..Mj,  (dj-Dj+Djq.-Djq),  1-1. ..Mj, 
(e^-E^p+E^q , E^p-E^q ) ,  1*1. ..M3, 

(b1+Bi-B1p,Bip),  1-M1+1...Q1,  (dj-D^q.-Djq),  i»M2+l...Q2, 
(e^-Ejp+E^q ,  E^p*E^q ) ,  i*M3+l . .  .Q3 


Proof  for  Theorem  5.2  From  Theorem  2.4 


\TV^Tl?WV>  I 


/  ft. 


fz(z)  «  Mx  [  Mps_p+i ,qs-q+l^  fX,Y(x,y)  '  1 

“Ml"1{  Mps-p+l ,qs-q+l{  k2Ht8lx‘82yJ  }  1 


and  from  Equation  (3.21) 


fz(z)  -  Mj  1  j  kg1-ug2-vx1(u)x2(-v)x:J(u-v)j 


u*ps-p+l 

v*qs-q+l 


M!  1  kg1~ug2Vxi(u)x2(v)x3(w)| 


u-ps-p+1 

v*-qs+q-l 

w»(p-q)s-p+q 


fz(z)  -  k8lP”l82q_1(2xr)  f  x1(u)x2(v)x3(w)  (glpg2qz)  s  ds  (5.10) 

'  /  h~l* 

u»ps -p+1 

v“-qs+q-l 

w-(p~q)s-p+q 

from  the  univariate  H-function  probability  density  function 
definition,  (4.1),  and  the  univariate  H-function  definition,  (3.1), 
cases  I,  II,  III,  and  IV  immediately  follow. 

Special  Cases 

A  simple  product  and  quotient  of  2H[x,y]  variates  are 
presented  as  special  cases  of  Theorem  5.2. 

The  distribution  of  Z  ■  XY:  From  Theorem  5.2,  case  III,  p  ■  q  *  1, 


(5.9)  reduces  to 


1 1“  v 


1  H*  i*  n 


7*r7T 


123 


f7(*> 


where 


K  H 


M1+N2,N1+f!2 

p1+Q2,Q1+p2  L 


glg22 


otherwise 


M  N3 

h  r (e , )  n  r(i-f  ) 

.  1  1 

r2 - Ql 

II  r(f  )  n  r(i-e.) 

n3+i  m3+i 


'2  J 


,  z>0 


(5.11) 


*  (a^,A^),  1*1. . .N^,  (l~d^,Dj),  i*l...M2, 

(ai,Ai),  i*N1+l. .  .P^,  (l-df,D1),  i-M2+l...Q2 
n2  *  (bi,Bi),  i-l...Mlt  (1-c^Cj,),  1*1. ..N2, 

(bj.Bj),  1*M1+1 • • .Qj ,  (l-ci,Cl),  i-N2+l...P2 

Equation  (5.11)  can  also  be  derived  directly  by  applying 
(3.21)  and  (2.23)  to  the  bivariate  H-function  distribution  definition, 
(4.2)  and  the  H-function  definition  (3.12). 

The  distribution  of  Z  ■  X/Y:  From  Theorem  5.2,  case  IV,  p  *  1, 
q  *  -1,  (5.9)  reduces  to 


M*)  “ 


K  H 

l2 

g2 


M  N 


P  Q 


-2  -i 


,  z>0 


(5.12) 


otherwise 


K-k, 

,  M- 

N-El^ , 

P-£Pi . 

Q-EQi,  i 

L«l,2, 3 

(a^ 

. • • , 

<ci-2Ci 

•^i ) ,  i-1 • • • N2, 

<fi‘ 

’2Fi 

,2F1) 

,  i-1.. 

•N-j,  (a 

i.A^,  j 

L-Nj+1.. 

•Pl* 

(ci' 

■2Ci 

,c± ), 

i»N2+l 

. . . P2 , 

(fi-^, 

,2Fi),  i< 

-n3+i., 

.  .p3 

(br 

,Bt) 

,  1-1 

<di'2Di 

,Di ) ,  i-1 . . «M2, 

(V 

-2Ei 

,2Et) 

,  i-1.. 

.m3,  (b 

i.Bj),  J 

L*M1+1.. 

•Qj  , 

(di- 

-2Di 

.Di), 

i-M2+l 

. . . q2 , 

(e^Ej, 

,2^),  1' 

-m3+i  . , 

•  *Q3 

Equation  (5.12)  can  also  be  derived  directly  by  applying 
(3.21)  and  (2.24)  to  the  bivariate  H-function  probability  density 
definition  (4.2)  and  to  the  H-function  definition  (3.12). 

Example  5.4  Consider  the  Kellogg-Barnes  111  distribution  given  by 

fv  v(x,y)  -  .xc  exp ( -ox  -  Bxy)  ,  x,y>0 

r(c)  a,  (5>0  ,  c>2 

(c,l) 

«  - ;  (1,1) 

y  - 


Applying  Theorem  5.2,  case  III,  p-q-1,  or  applying  (5.11)  directly 
the  distribution  for  the  random  variable  Z-XY  is 


The  distribution  of  fz(z)  can  be  found  directly  from  (5.13)  by 
realizing  that  (5.13)  is  the  H-function  identity  for  the  exponential 
distribution  with  parameter  1/3  as  given  by  Carter  (3).  Then 

fz(z)  -  Be"02 

which  agrees  with  the  solution  obtained  in  Example  1.4  using  the 
convolution  integral. 

5.3  Transformations  for  Two  or  More  Bivariate  Distributions 

In  the  last  section  theorems  were  presented  deriving  the 
distribution  for  a  random  variable  resulting  from  some  combination  of 
dependent  H-function  variates  sharing  a  bivariate  H-function 
probability  density  function.  This  section  is  concerned  with 
procedures  for  finding  the  bivariate  probability  density  functions  of 
random  variables  which  are  the  products  or  ratios  of  other  bivariate 
random  variables. 


5.3.1  Tranformatlons  for  ^Hfx,y]  Variates 

5. 3. 1.1  The  Distribution  of  Products 

Theorem  5.3:  If  X^,Y^;  ***»  Xn,Yn  are  n  Pa*rwise  Independent 

random  variables  with  bivariate  probability  density  functions 
fl(xl‘yl>*  f2<x2»y2> . fn<xn»yn>  respectively,  where 


lH[g1jXj.g2jyj) 


Vyj >  0 


otherwise 


for  j*l,...,n,  then  the  bivariate  probability  density  function  of  the 
random  variables 


n  n 

z-n  X  ;  W-H  Y 
j-1  2  j-1 


is  given  by 


fZ,W<z-w> 


l(Ilkj>  1H 


^^1  j’^^l  j’^2  j’^^2  j»^*^3j*^^3j 


^^1  j»JQi  j»£p2j'Z(52j*ZP3j,S:Q3j 


(nglj)z 

(ng2j)w 


^2  *  ^3 


n,  ;  n 


5  *  '6 


z,w> 


0 


otherwise 


(5.14 


where  the  E's  and  H’s  are  indexed  over  j  and  go  from  1  to  n  and 


nl 

- 

(eij'Ei j>* 

i  *  1  •  •  •  M  j  j 

j-1 « « *n , 

(eij»Fij),  i-M^j+l. 

**Q3j. 

j-1. .  .n 

n2 

- 

(ai j'Ai j). 

j-1 • • *n , 

(aij.Aij),  i-Njj+1. 

’  *P1  j* 

j-1. .  .n 

n3 

- 

(Cij.cij). 

i-l • • • ^2 j  * 

j  * 1 . . .n , 

(cij,^ij),  i-®2j^"l» 

*  *P2 j* 

j-1. .  .n 

\ 

- 

(fij.Fij), 

i*l  •  •  »N<j  j» 

j-1* • *n. 

(fij,Fij>,  i-N3j+l. 

*  *P3  j* 

j-1. • *n 

n5 

at 

(bi j,Bij)* 

j— 1 « « *n , 

(bij,Etj)*  l—^l j^l* 

•*Qlj, 

j-1. • .n 

n6 

- 

(dij.Dij). 

l *1 ♦ • .M2 j  f 

j— 1. . *n, 

(d i j » Di j ) »  I" ^2 j+l * 

*  *Q2 j* 

j-1. • *n 

Proof  of  Theorem  5.3s  From  Theorem  2.5, 


ZfW(2,v)  -  Mf1  f  jfxj.yj)  >  ]  z,w  >  0 


and  from  the  Mellin  transform  of  the  bivariate  H-function,  (3.19),  the 
bivariate  density  of  Z  and  W  is  then  given  by 


fZ,w(z,w)  *  M2 


-1 


n 

[  n  {  k 

j-i 


-Si 


-s. 


j®lj  ®2j  X1 j^sl^x2 j^s2^x3 j(si+s2^  ^  1 


(  n  k.)  M 


-1 


j-1 


j‘ 


n 

(  n  { 

j-i 


-S, 


C1 j(si)x2j(s2^x3j^sl+s2^8l j  g2j 


n 


where  x]j(Sj),  X2j(s2^’  and  x3j(si+s2)  are  defined  by  (3.7),  (3.8), 
and  (3.9)  respectively.  Recombining  like  product  terms  and  using 
(3.6),  equation  (5.14)  may  be  written  directly. 

For  the  case  where  Mij"Ni j“pi j-Q* j"°,  1-2.3,  (or  for  i-1,3). 


j»l,...,n,  the  bivariate  H-functions  reduce  to  univariate  H-functions 
and  (5.14)  reduces  to  Carter's  Theorem  4-1,  (3:52),  for  the 
distribution  of  products  of  independent  H-function  variates. 

5. 3. 1.2  The  Distribution  of  Quotients 

Theorem  5.4:  If  X^,Y^;  X2,Y2  are  two  pairwise  independent  random 
variables  with  bivariate  probability  density  functions  ^(x^yi)  and 
^2^x2,y2^  respectively,  where 


f  kj  iHt8ijxj’*2/ji  •  xj»yj  > 0 

< 

[_  0  ,  otherwise 


for  j*l,2,  then  the  bivariate  probability  density  function  of  the 
random  variables 


Z  -  X^Xj  ;  W  - 


is  given  by 


kjkj  m1.N1,M2,N2,M3,N3 

i  H 


2  2  1 
12  22  ^l*Ql*^2’^2,*>3,^3 


fz,w(*.w) 


811 


®12 

®21 


g 


22 


n2  5  n3 
n4 

n5  5  % 


otherwise 


where  N^*N^j+M^2,  ^i,"^il"^i2,  an{^  ^i"^il^i2* 


z,w>< 


(5.15 


for  i-1,2,3 


il»Eil)* 

i*l . .  ,  (l*"f 

12'4Fi2»Fi2 

i 1  *  Eil 5 » 

i^M^^+l • • .Q^ , 

( 1-f i2~4Fi2 

il’Ail)' 

i-l...N11,  ( 1  — b 

12~2Bi2’Bi2 

il»AiP» 

i“N^^+l • • .Pj ^ , 

(1_bi2”2Bi2 

il.cil>. 

i-l...N21,  (1-d 

12-2Di2,Di2 

il»Cil)» 

i*^ ^1 . .  .P21 , 

(l-d12-2Di2 

il»Fil)* 

i*l. . .N^ ,  (1— e 

i2"4Ei2»Ei2 

il»Fil)* 

i"^31+l *  *  * F3 1  * 

( 1-ei2“4Ei2 

il»Bil)* 

i»l...M11,  ( 1-a 

i2~2Ai2,Ai2 

i  1  *  Bi  1  ^  * 

i-Mn+1. .  .Qn» 

(l~ai2“2Ai2 

il»Dil)» 

i«l. . .Mj^,  ( 1“C 

i2"2Ci2»Ci2 

il’0!!5’ 

*”M21+1 *  * ‘^21  * 

(l-ci2-2Ci2 

Proof  of  Theorem  5.4:  From  Theorem  2.5,  n*2,  a^bj-1,  and  a2*b2*-l» 
or  from  (2.22) 

fZ,W<z»w>  "  M2  1  [  Ms1,s2{  fl(xl,yl)  }  M2-s1,2-s2{  f2(x2»y2>  }  1 

and  from  Che  Mellin  transform  of  Che  bivariate  11-function,  (3.19),  the 
bivariate  density  of  Z  and  W  is  then  given  by 

—1  ~ ®i  ~®2 

fZ,W<z»w)  *  M2  ^  kl8ll  821  x11^8i)x21^82^x31^s1+s2^ 


s^-2  S2“2 

y  k2812  822  x12^”81^x22^"”s2^x32^~81”82^  1 


nT*1 


♦Tl 
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klk2  -1 

*  2  2  M2  ^  xll^sl^x21^s2^x31^s l"*"s2) 

g12  S22 

~S1  “s2 

X  xi2(2”si)x22^-s2^x32^-sl~s2^j>11^8l2^  ^®21^®22^  1 

where  x^f.),  x2j(*),  and  x3j(->,  j-1,2,  are  defined  by  (3.7),  (3.8), 
and  (3.9)  respectively.  Recombining  like  product  terms  and  using 
(3.6),  equation  (5.13)  is  obtained. 

For  the  case  where  Mi j“Nt j«Pt j-Q* j-0,  i-2,3,  j-1,2,  then  the 
bivariate  H-functions  reduce  to  univariate  H-functions,  and  (5.15) 
reduces  to  Carter's  Theorem  4-8,  (3:62),  for  quotients  of  independent 
H-function  variates. 

5. 3. 1.3  The  Distribution  of  Powers 

Because  of  the  restriction  in  the  bivariate  H-function 
definition  that  s^  and  Sj  must  have  identical  coefficients  in  the 
x3^sl+s2^  term*  a  general  H-function  of  the  form  jH[z,w]  for  powers  of 
H-function  variates  can  not  be  obtained.  However,  an  H-function  of 
the  form  ^H[z1^p,w^^] ,  where  p,q  are  rational  powers  of  x  and  y 
respectively,  is  possible.  For  the  special  case  where  p-q  an 
H-function  of  the  form  1H[z,w]  is  possible. 

Theorem  5.5:  If  X,Y  are  random  variables  with  bivariate  probability 
density  function  fx,Y^x,y^*  where 


fX>y(x,y) 


k  jHfg^.gjy]  ,  x,y  >  0 


0  ,  otherwise 


then  the  bivariate  probability  density  function  of  the  random 
variables 

Z  -  Xp  ;  W  -  Yq 
p,q  rational  and  p,q  >  0,  is  given  by 

fZ,W(z»vO  ■ 


(e1+(2-q-p)E1,E1) 

1/p 

*1* 

^ai+Ai“pAi,Ai^» ^ci+Ci-qCi,Ci^ 

a  „1/q 
g2 

(fi+(2-q-p)Fi,Fi) 

^bi+Bi~pBi » Bi ) * (di+Di-qDi ,Di ^ 

(^0  otherwise 

where  K  -  kg1p-1g2q-1/pq 

Proof  of  Theorem  5 » 5 :  From  Theorem  2.3,  letting  tj«s^,  t2*s2' 

fZiU(.,»)  •  «2  '[  Mptl-p+l,qt2-<|+ll  1  1 

*  »2_lf  V^l.qtj-q+l'  klHt8iX.82y]  J  ] 
and  from  equation  (3.19) 


fZ,W(z,vO  “  M2  kgx  g2  x1(u)x2(v)x3(u+v) 


u*pt j-p+1 
v*qt_-q+l J 


— t i  ~t  2 

fz>w(z,w)  "  K  —  2  x1<u)x2(v)x3(u+v)i  (81Pz)  (g2qw)  dtidt2 


<2iti) 


u»pt1-p+l 

v«qt2-q+l 


where  K  *  kg^'^g^-^  and  //  is  the  double  contour  integral  for  the 
inverse  Mellin  transform  as  defined  by  (2.6).  Let  s^*pt3,  s2*qt2,  and 
dSjds2*pqdt jdt2,  then 


-s./p  „  -s7/q 


fZ  W(2.W)  -  -  — 1— . Jf  X1(U)x2(v)x3(u+v)|(g1P2)  1  (g2qw)  2  dSlds2 


pq  (2ni) 


U-Sj-p+1 

I  v-s2-q+i 


-  - _//  x,(u)x2(v)x3(u+v)  (g,z1/p)  (g2w  /q)  ds^sj 

pq  (2»i)2 

U-Sj-p+1 

v-s2-q+l 


from  which  (5.16)  follows. 

For  the  case  where  p-q,  then  from  (3.15),  the  H-function  in 
equation  (5.16)  becomes 

(e ^+2E ^-2pEj , pE^ ) 

M^IN3,M2,N2,M3,N3  g]^z  (s^+A^-pA^ ,pA^ ) ; (c^+Cj-pC^ ,pC^ ) 

KlH  *.w>( 

pi»Qi,p2.Q2.p3.Q3  g2Pw  <  V2Fi-2pFi*PFi> 

( bi+Bi  -p Bj ,  p Bj  ) ;  ( d i+Di  -pDx ,  pDj^ ) 


,1 


7»V  >*»” 


where  K  - 


(5.17) 


If  (2.23)  is  applied  to  (5.17)  above,  the  distribution  of 
Z'-(ZW)P  may  be  found.  It  can  be  shown  that  such  a  procedure  will 
result  in  a  solution  identical  to  that  given  by  Theorem  5.1,  case  I, 
for  pmq. 

Example  5.5:  Consider  the  Kellogg-Barnes  I  distribution  given  by 


fX,Y(x,y) 


4a 


6+1 


irr(6+l) 


(x2  ♦  y2)6  + 


x,y  >  o 


(8,1/2) 


a 

u  T ( 8+1 ) 


1,0,1, 0,1,0 

1H 

0,1, 0,1, 1,1 


/a  x 
/a  y 


(0,1/2) 

(0,1/2)  ;  (0,1/2) 


Now  suppose  the  bivariate  probability  density  function  of 
Z  *  X^/Xj,  W  «■  Y^/Y2  is  desired  where  X^,Y^  and  X2»Y2  are  distributed 
according  to  the  Kellogg-Barnes  I  distribution  given  above.  Suppose 
further  that  f^(x^,y^)  has  parameters  ■  1,  8j_  ■  0,  and  f2(x2»Y2) 
has  parameters  a2  -  1,  02  -  0.  From  Theorem  5.4,  the  bivariate 
probability  density  function  of  Z,  W  is  given  by 


^  ^  r  *  .  «r.  «r.  j-.  ^v.  -'.  *  ,  --.v.  •-.  -v  “, 


(-1,1/2) 

1,1,1, 1,0,1  z  (0,1/2)  ;  (0,1/2) 

_  ~ 

TT2r(6+l)  1  1,1, 1,1, 1,1  W  (-1-0, 1/2) 

(0,1/2)  ;  (0,1/2) 


(5.18 


from  property  (3.15),  k  =  2,  Equation  (5.18)  above  roay  be  written  as 


(-1,1) 

,  1,1, 1,1, 0,1  z2  (0,1)  ;  (0,1) 

4  H 

u2r (8+i)  1  1,1, 1,1, 1,1  W2  (-1-0,1) 

(0,1)  ;  (0,1) 


(5.19 


Recognizing  that  (5.19)  is  the  H-function  form  for  one  of  Appell's 
hypergeometric  functions,  special  case  (3.22),  equation  (5.19)  may  be 
written  as 


f7  u(z,w)  -  -  r(0+2)r(l)r(l)F  (2+8,1, 1;2;-z2,-w2) 

nzr(0+i) 


4(-g^--)-  F1(2+8,1,1;2;-z2,-w2) 


where  is  Appell's  hypergeometric  function  of  two  variables  as 
defined  in  Appendix  B.  The  results  obtained  using  the  bivariate 
H-function  and  Theorem  5.4  agree  with  results  obtained  using  Mellin 
transform  techniques  as  outlined  by  Fox  (10)  and  illustrated  in 
Example  2.3. 
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Example  5.6:  Consider  Kellogg-Barnes  II  distribution  given  by 


fXfY(x,y) 


ga2  e-(ax  +  By/x) 


x,y  >  0 


(0,1) 


6a2  XH 


0,0, 1,0, 1,0  ax 
0,0, 0,1, 0,1  agy 


(5.20) 


*,  (0,1) 


From  (5.17),  the  bivariate  probability  density  function  of  Z-X1^2, 


W»Y1/Z  is  given  as 


0, 0, 1,0, 1,0  fa  2 


f  (z,w)  «  H 

’  0,0, 0,1, 0,1  /ag  w 


(1,1/2) 


(5.21) 


;  (1/2, 1/2) 


From  (3.15),  for  k«2,  (5.21)  can  represented  as 


(1,1) 


1/2  0,0, 1,0, 1,0  az 

fZ  W(z»«)  -  ™  , 

*  0,0,0, 1,0,1  a0w* 


(5.22) 


!  d/2,1) 


Equation  (5.22)  above  may  also  be  obtained  directly  from  Theorem  5.5, 
p»q*l/2.  Using  H-function  property  (3.16),  m»n»l/2,  (5.22)  may  now  be 
expressed  as 


fZ,W(2*v0  “  46a  zv 


0,0, 1,0, 1,0 


0,0,0, 1,0,1 


(0,1) 


az 


a(Sw * 


;  (o,i) 


(5.23 


Comparing  the  H-function  in  (5.23)  with  the  H-function  in  (5.20),  the 
distribution  for  Z,  W  may  be  written  as 


fZ,W(z,w) 


48a2  zw  f  e“(au  +  6v/u>  ] 


u»z. 


v»w 


46a2  4.  *-<“2  +  6<"/2)2> 


The  solution  obtained  above  may  also  be  found  using  Theorem 


1.4.  Let  Z*X1/2  and  W-Y1/2,  then  X-Z2,  Y-W2,  and  the  Jacobian  is 


3x/3z  3x/3w 

3y/3z  3y/3w 


2z  0 
0  2w 


4zw 


From  Theorem  1.4, 


fZ,W^z,w^  “  fX,Y*z  »w  H  J 


48a2  zw  e-(°z2  +  0(w/z)2) 


5.3.2  Transformations  for  ^H[xty ]  Variates 

This  section  presents  theorems  for  transformations  of  2H[x,y] 
variates.  Examples  for  testing  the  theorems  in  this  section  are 
difficult  to  derive,  and  unlike  the  last  section  cannot  be  checked  by 
comparing  the  results  given  here  to  those  obtained  for  independent 
H-function  variates  by  previous  authors.  This  does  not  mean  that  the 
following  theoretical  results  have  no  practical  applications.  Rather, 
the  reverse  may  be  true  -  this  may  be  the  only  practical  means  for 
obtaining  the  distribution  for  transformations  of  the  type  given  in 
this  section. 

5. 3. 2.1  The  Distribution  of  Products 

Theorem  5.6:  If  X^.Yjj  X2»Y2;  •*•;  xn»Yn  are  n  pairwise  independent 
random  variables  with  bivariate  probability  density  functions 

fl(xl,yl)»  f2(x2,y2)»  *•*»  fn(xn,yn5  respectively,  where 


kj  2Hfgljxj,82jyj!  ’ 


Xj,yj  > 


f j<XJ»yj)  " 


0  ,  otherwise 


for  j»l,...,n,  then  the  bivariate  probability  density  function  of  the 
random  variables 


n  n 

Z  -  il  X  ;  W  -  n  Y. 

j-1  3  j-1  J 


is  given  by 


'  .’j.  ■.<  v 


z,w>0 


(Hkj)  2H 


EMlj’ENlj*EM2j*EN2j*EM3j*EN3j 
EPlj*EQlj’EP2j,EQ2j‘EP3j‘EQ3j 


(ng!j)z 

(ng2j)« 


otherwise 


(5.24) 


where  the  I's  and  II' s  are  indexed  over  j  and  go  from  1  to  n  and 


nl 

* 

(eij,Ei j). 

i*l • • j  * 

j-l. . .n. 

(eij*Eij)» 

i«M3y4-l. 

*  *Q3 j* 

c 

• 

• 

• 

H 

*»“> 

n2 

m 

(ai j,Ai j) * 

i-i— "u. 

j”l • . *n , 

(ai j,Ai j>. 

i-Nx j+1. 

•*Plj* 

j*l. . .n 

n3 

m 

(Cij,cij), 

i-l...N2j, 

j*l • • . n , 

(c1;j,Cij), 

i-N2j+l. 

..p2j. 

j-1. . .n 

n4 

“ 

<fij,Fij>* 

i*l • • • N3  j  f 

j-l...n. 

<fij’Fij>’ 

i-N3j+l. 

’  *P3  j* 

j-1. . .n 

n5 

- 

(bij»Bij). 

i • 1 • • •Mljf 

j«l. . .n. 

(bij,Bi j)’ 

i-Mj^j+l. 

j-1. . .n 

n6 

* 

(di j,Di j)’ 

i*l  •  •  *  **2  j  f 

j*l • • .n , 

(dij,Di j). 

i«M2j+l. 

• *Q2 j* 

j-1. . .n 

Proof  of  Theorem  5.6:  From  Theorem  2.5, 


fZ,W<z’w> 


n 

[  il  M 

j-1 


®  1  |  S  r> 
1  *  l 


*  i 


z,w  >  0 


and  from  the  Mellin  transform  of  2H [x»y ] •  (3.21),  the  bivariate 
density  of  Z  and  W  is  then  given  by 


4 

1 

* 


\ 
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-1  n  -S1  ~s2 

^Z,W^z,w^  M  ^2  f  ^  ^  ^j®lj  82j  X1 j^sl^x2 j(~s2^x3j^sl_s2^  ^  1 

^  Xlj(8l)x2j<“82)x3j(8l“s2)8lj  g2j 

where  x^s^,  x2j(-s2),  and  x^^-Sj)  are  defined  by  (3.7),  (3.4), 
and  (3.20)  respectively.  Recombining  like  product  terms  and  using 
(3.12),  Equation  (5.24)  may  be  written  directly. 

For  the  case  where  Mi  j*Pi  j-Qj  j"0,  i-2,3,  (or  for  i»l,3), 
j»l,...,n,  the  bivariate  H-functions  reduce  to  univariate  H-functions 
and  (5.14)  reduces  to  Carter's  Theorem  4-1,  (3:52),  for  products  of 
Independent  H-function  variates. 

Comparing  Theorem  5.6  with  Theorem  5.3,  it  is  seen  that  the 
form  of  theorems  is  identical  except  that  the  first  deals  with  2H[x,y] 
variates  while  the  second  deals  with  jH[x,y]  variates.  It  is  clear 
that  a  more  general  theorem  can  be  established  for  both  ^H[x,y] 
variates  and  2H[x,y]  variates.  However,  this  property  does  not  hold 
for  |H[x,y]  variates  and  2H[x,y]  variates  when  transformations  of 
ratios  or  powers  of  the  variates  are  concerned. 

5. 3.2.2  The  Distribution  of  Quotients 

Theorem  5.7:  If  X^,Y^;  X2,Y2  are  two  pairwise  independent  random 
variables  with  bivariate  probability  density  functions  f^(xj,y^)  and 
^2^x2»^2^  fespectlvely*  where 


■ft 


ft 


&<$ 


k*. 


xj’yj  >  0 


fj(xj’yj} 


j  kj  2lIlgljxj'82jyjJ  » 
|_  0  ,  otherwise 


for  j-1,2,  then  the  bivariate  probability  density  function  of  the 
random  variables 


Z  -  Xj/Xj  ;  W  -  Y1  / Y. 


V  ‘2 


is  given  by 


klk2  mi.N1,M2,N2,M3,N3 


,H 


2  2  2 
®12  g22  P 


l’^l,P2*^2*P3,^3 


fZ,W<Z’w)  * 


gll 


g12 

g21 


g 


{  n- 


22 


n5  5  n6 


otherwise 


z,w> 


(5.25 


where  and  Qi“Qil+Pi2’  for  i"1»2*3 

and 


■  ^eil*^il^*  ^2 ,Pi2 ^ ’  i*l. • *N32 

^eil*^il)*  ^“^31"*’^- •  •  *^3i *  ^”^i2,Pi2^’  *“^32+*  *  *  *P32 
^2  ”  ^ail’^il^’  ( l”b£2-2B^2, Bj2) ,  i“l...M^2> 

Caii*Aii),  i*Nj^+l . .  •  ( l“l)j2”2 ®i2 ’ ®i2 ^ ’  *"M12+*  *  *  *^12 

h3  “  ^cil’^il^’  i“l* • *^21 »  d“di2+2Di2*  ^i2^*  i*l«.*M22 

^cil’^il^’  ^*^2i"^l* * *p2l>  (I“di2+^^i2*®i2^*  *”M22+* ••*^22 


ol 


(  f  i  1 1 

•Fil>’ 

i-1.. 

.N31 

,  (1-e 

12,Ei2^’  • • *M32 * 

(fil’ 

•Fil>’ 

i-N3i 

+  1  •  • 

,P31  * 

(l~e^2»Bi2^*  1^32^1'" 

(bii- 

,Bil), 

i-1 . . 

•  Mu 

,  (1-a 

i2~2Ai 2 , Aj 2 ) ,  i-1 • « 

*N12 

<bi  1  ’ 

.Bil), 

i-Mn 

+  1  •  •  1 

,Q11* 

( l~a^2~2Aj2 »^i2 ) »  1 

=N12 

(dil' 

,Dil > » 

i-1.. 

*M21 

,  (1-c 

i2+2Ci2,Ci2^’  i“1*’ 

*N22 

<dil 

,Dil), 

i*M21 

+  1  •  • 

•Q2i« 

( l~c £2^"2Ci2 ,Ci2 ) ,  i 

"N22 

Proof  of  Theorem  5.7:  From  Theorem  2.5,  n«2,  a^*b^«l,  and  a2mk2— 1» 
or  from  (2.22) 

fZ,W(z’w)  "  M2  1  1  Ms1,s2{  fl(xl*yl)  }  M2-Sl,2-s2{  f2(x2»y2)  }  1 

and  fr.j  the  Mellln  transform  of  2H[x,y],  (3.19),  the  bivariate 
density  of  Z  and  W  is  then  given  by 

-1  -S1  ~s2 

FZ,W^z,w^  *  M2  ^  ^1®11  ®21  *11 ^ 8 1 ^x21 ^ — s2 ^x31 ^8 l~s2 ^ 

s^-2  82*2 

X  ^2®12  ®22  x12^2~sl^x22^82~2^x32^""sl+s2^  1 


^1^2  -1 

2  2  M2  l  xn(Sl)*21^-s2^x31^sl”82^ 

812  822 


-S1  -s2 

X  xi2^"'8l^x22^s2-^^x32^~sl+s2^®ll^®12^  ^821^®22^  ^ 


where  x^j(«),  Xjjf*),  and  x^jf-),  j-1,2,  are  defined  by  (3.7),  (3.8), 
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and  (3.9)  respectively.  Recombining  like  product  terms  and  using 

(3.12),  Equation  (3.25)  is  obtained. 

For  the  case  where  j*N1 j»P1 j-Qi j-0,  i-2,3,  j-1,2,  then  the 

bivariate  H-functions  reduce  to  univariate  H-functions,  and  (5.25) 

reduces  to  Carter's  Theorem  4-8,  (3:62),  for  the  distribution  of 

quotients  of  independent  H-function  variates.  Unlike  Theorem  5.4, 

this  property  does  not  hold  for  i»l,3.  This  is  due  to  the  fact  that 

s2  “8  2 

for  i“2,  the  H-function  is  defined  as  y  instead  of  y  . 


h  ' 


5. 3. 2. 3  The  Distribution  of  Powers 

For  the  same  reasons  given  in  section  5. 3. 1.3,  an  H-function 
of  the  form  ^H[z,v]  is  not  generally  obtainable  for  rational  power 
transformations  of  the  H-function  variates  X  and  Y.  However,  an 
H-function  of  the  form  jHfz^^.w*^ ]  is  obtainable  and  for  the  special 
case  where  p-q  an  H-function  of  the  form  jHfz.wJ  can  be  derived. 
Theorem  5.8:  If  X,Y  are  random  variables  with  bivariate  probability 
density  function  f^  y(x,y),  where 


8 


^V- 


VVj 


fX,Y(x»y) 


k  2Hf8ix»s2yl  »  x»y  >  0 


0  ,  otherwise 


then  the  bivariate  probability  density  function  of  the  random 
variables 


p,q  rational  and  p,q  >  0,  is  given  by 


fZ,W(z,w) 


m1,n1,m2,n2,m3,n3 

Pl,Ql*^2*^2,P3*^3 


§lz 

g2w 


1/P 

i/q 


(ei+(q-p)Ei»Ei) 

(a^+A^—pAj (Aj ) ; (c^—C^+qC^  *C^ ) 
(f1+(q-p)Fi,F1) 

(bi+Bi -pBj , Bt ) ; (d1-Di+qD1 , ) 


z,w>0 


otherwise 
where  K  -  kgjP^g^^/pq 

Proof  of  Theorem  5.8:  From  Theorem  2.3,  letting  tj-Sj,  t2«s2, 

-1, 


(5.26) 


f7(z)  -  M,  *[  M 
-1 


Zv  '  “2  1  ‘  ptrp+l,qt2-q-fll  X, Y 


t  fY  Y(x»y)  }  1 


'  M2  I  k2HiBlx.e2y)  )  1 

and  from  equation  (3.21) 
fZ,W<z»w> 


-  M2_1 

kgl  Ug2  Vxi(u)x2^_v^x3^u_v^ 

u-ptj^-p+1 

- 

v»qt2-q+l  -* 

fZ,W<z’W> 


(2wi  y 


IS  x1(u)x2(-v)x3(u-v) 


-t 

(g2qw) 


2 


dtjdt 


2 


u-ptj^p+1 

v«qt2-q+l 


where  K 


vvl 


and  IS  is  the  double  contour  integral  for  the 


an 


inverse  Mellin  transform  as  defined  by  (2.6).  Let  s^ptj,  s2«qt2. 


ds1ds2-pqdt1dt2,  then 


fZ,w(z,w)  x.(u)x2(-v)x-(u-v) 

*  pq(2*i)2  1  2  3 


-s,/p 


p.,  i.  ■s’/'' 


(gx  z>  (g2  w)  dSjds. 


u*SfP+l 

v»s2-q+l 


»  -  — // 
pq(2ui)Z 


xl(u)x2(~v)x3(u-v) 


(g,z  ,H) 


1/p 


“8 


1(g2»Uq)  ^dSjds. 


u-s1-p+l 

v-s2-q+l 

u-v-sl-s2+q-p 


from  which  (5.26)  follows. 

For  the  case  where  p«q,  then  from  (3.15),  the  H-function  in 
Equation  (5.26)  becomes 


m1,n1,m2,n2,m3,n3 

*^l»  ^2  ‘^2’  ^3*^3 


(e^ ,pEj ) 

(ai+A1-PAi  ,pAt ) ;  (^-Cj+pt^  .pCj  ) 


<bi+Bi -pBj .pBj ) ; (dj-Dj+pDj ,pD1 ) 


z,w>0 


where  K  -  Mg^)?'1  (5.27) 

If  (2.23)  is  applied  to  (5.27),  the  distribution  of  Z'-(ZW)P 
may  be  found.  It  can  be  shown  that  such  an  application  will  result  in 
a  solution  identical  to  that  given  by  Theorem  5.2,  case  III. 


I 


Example  5.7 

fx,Y<x*y>  - 


Consider  the  Kellogg-Barnes  III  distribution  given  by 
Xc  e"<«  +  x,y>0 


I 

I 


0 


r<c)  2 


H 


0,0,0, 1,1,0 
0,0, 1,0,0, 1 


ax 


(c,l) 

-  ;  (i,D 


(5.28) 


From  (5.27),  the  bivariate  probability  density  function  of  Z-X1^, 

w,yl/2  ls  gjven  as 


fZ,W(z»w) 


gl/2  0,0, 0,1, 1,0 
^C)  2  0,0, 1,0, 0,1 


/a  z 


7Z. 


w 


(c , 1/2) 

-  ;  (1/2, 1/2) 


(5.29) 


From  (3.15),  k*2,  equation  (5.29)  may  be  written  as 


fz,w(z,w) 


0,0, 0,1, 1,0 
0,0, 1,0, 0,1 


(c.l) 

-  5  d/2,1) 


(5.30) 


Equation  (5.30)  may  also  be  obtained  by  application  of  Theorem  5.8, 
p*q*l/2.  From  (3.17),  m«n*l/2,  (5.30)  may  be  written  as 


fZ,W<z’w> 


0,0, 0,1, 1,0 

1 

0,0, 1,0, 0,1 


az 


-  ;  (i,D 


(5.31) 


43 

r(c) 


2W  2H 


iw2 

a 


Comparing  (5.31)  with  (5.28),  the  bivariate  probability  density 
function  for  Z  and  W  may  be  given  by 


fZ,W<z*w> 


*£  zw  [  ac  uc  e_(au  +  Suv> 

r(C) 


1 


u«z 

v-w 


2 

2 


.  Z>BaC  wz2c+1  e"'(Qlz2  +  S(zw)2) 

F(c )  WZ 

The  solution  obtained  above  may  be  obtained  through  the  use  of 
Theorem  1.4.  Let  Z-X1/2  and  W-Y1^2,  then  X«Z2,  Y-W2,  and  the  Jacobian 
is 


3x/3z 

3x/3w 

2z  0 

J  - 

* 

3y/3z 

3y/3w 

0  2w 

From  Theorem  1.4, 

fZ,W<z,w)  "  fX,Y^z2,w2^  lJl 


4Bac  2c+l  -(az2  +  8(zw)2) 

r<Z) 


5.3.3  The  Distribution  of  a  Mix  of  Product  and  Quotient 


In  the  previous  two  sections  the  bivariate  distributions  of 
random  variables  which  were  products  or  quotients  of  pairwise 
independent  random  variables  from  two  or  more  bivariate  distributions. 
While  no  theorems  are  presented,  this  section  discusses  the  bivariate 
density  functions  for  random  variables  that  result  as  a  mix  of 
products  and  quotients  of  H-function  variates. 

As  an  example,  suppose  the  bivariate  density  function  of  the 
dependent  random  variables  Z»X  and  W-l/Y  is  desired.  From  (3.13)  it 
is  clear  that  if  X  and  Y  are  ^[x.y]  variates,  then  the  resulting 
bivariate  density  function  is  of  the  form  2H[x,y).  Conversely,  if  X 
and  Y  are  2H[x,y]  variates,  then  the  resulting  density  function  is  of 
the  form  jH[x,y].  Theorems  for  such  transformations  can  be  readily 
deduced,  but  are  not  presented  here  since  they  are  not  readily 
generalized  to  two  or  more  bivariate  H-function  distributions  as  is 
shown  in  the  following  paragraphs. 

Now  suppose  the  bivariate  density  function  of  the  dependent 
random  variables  Z»X^X2  and  W»Yj/Y2»  where  XpYj  ;  X2»Y2  are  pairwise 
independent,  is  desired.  From  Theorem  2.5,  the  bivariate  density 
function  of  Z  and  W  is  given  by 


fZ,W(z*w>  "  M2_1  f  Ms1’82{  fl(*l,yi>  }  Msi.2-s2{  f2(x2’y2)  }  ] 


If  Xj,Yj  »  are  iH(x,yJ  variates,  then  from  (3.19)  the 
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distribution  of  Z  and  W  is  given  by 
-1  -si  ~s2 

^Z,W^z,w)  "  M2  ^  ^1®11  ®21  xll(sl^x21's2^x31^sl+s2^ 

s^-2  s2~2 

X  ^2®12  ®22  x12^sl^x22^-s2^x32^+sl-s2^  1 

Comparing  the  x3i(si+s2)  term  and  the  x^^+s^-Sj)  term,  it  is  clear 
that  it  is  not  possible  to  make  a  change  of  variables  substitution  and 
meet  the  restriction  in  the  ^[x.y]  definition,  (3.6),  that  terms 
containing  both  s^  and  S2  must  have  s^  and  S2  be  of  the  same  sign.  It 
also  fails  to  meet  the  restriction  in  the  2H[x,y]  definition,  (3.12), 
that  s ^  and  s2  must  be  of  opposite  sign  for  terms  containing  both  s^ 
and  s2.  The  same  problem  occurs  if  Xj^Y^  ;  x2,Y2  are  2H[x,yl 
variates.  If,  however,  X^Y^  are  jH[x,y]  variates  and  X2,Y2  are 
2H[x,y)  variates,  then  the  density  function  of  the  random  variables  Z 
and  W  is  representable  as  a  bivariate  H-function  and  will  be  of  the 
form  jH[z,w).  Conversely,  if  X^,Y^  are  2H[x,y]  variates  and  X2,Y2  are 
^H[x,y]  variates,  then  the  density  function  of  Z  and  W  is  also 
representable  as  a  bivariate  H-function  and  will  be  of  the  form 

2H[x »y ] • 

To  eliminate  the  problems  described  above,  a  more  general 
definition  of  the  bivariate  H-function  would  need  to  be  introduced. 
Such  a  definition  would  be  similar  to  definition  (3.6)  but  would  have 
an  added  x^(s^-s2)  term.  Civen  such  a  definition,  it  can  be  shown 
that  jH[x,y]  and  2H[x,y]  would  be  special  cases  of  this  more  general 


definition.  Therefore,  Theorems  5.3  and  5.6  would  be  special  cases  of 
a  more  general  theorem  for  products  of  H-function  variates. 

Similarly,  Theorems  5.4  and  5.7  would  just  be  special  cases  of  a  more 
general  theorem  for  quotients  of  H-function  variates.  The  primary 
difficulty  with  such  a  definition  is  the  problem  of  identifying  the 
appropriate  contours  in  the  inversion  integrals  to  separate  poles  in 
the  left  and  right  half  planes  for  the  s^  and  Sj  variables.  While  it 
would  seem  that  such  a  definition  is  viable,  additional  study  on  this 
problem  would  be  required  befor  the  definition  could  be  implemented. 


CHAPTER  6 


Evaluation  of  the  Bivariate  H-function 
6.1  General  Remarks 

In  Chapter  5,  it  is  shown  that  the  probability  density 
function  of  the  product  or  ratio  of  two  dependent  H-function  variates 
is  a  univariate  H-function.  Further,  it  is  shown  that  the  univariate 
H-function  results  from  simple  algebraic  manipulations  of  the 
parameters  of  the  bivariate  H-function  probability  density  function. 
Then  to  find  the  density  function  of  the  product  or  ratio  of  two 
depndent  H-function  variates,  one  need  merely  find  the  univariate 
H-function  representation  by  application  of  Theorems  5*1  or  5.2.  The 
computer  program  given  by  Cook  (5:154-176)  may  then  be  applied  to  give 
a  numerical  evaluation  and  plot  of  the  univariate  H-function  density. 
For  these  types  of  problems,  evaluation  of  the  bivariate  H-function 
density  is  neither  required  nor  desired. 

For  problems  where  the  bivariate  density  function  of  products 
or  quotients  of  pairwise  independent  variates  from  two  or  more 
bivariate  density  functions  is  desired,  evaluation  of  the  bivariate 
H-function  is  required.  The  procedure  for  such  an  evaluation  is  to 
invert  the  Mellin  transform  of  the  bivariate  U-function  by  applying 
the  residue  theorem  in  an  iterative  fashion.  The  procedure,  as 
outlined  in  the  following  sections,  is  based  on  the  assumption  that 
any  bivariate  probability  density  function  of  interest  is  continuous 
in  each  variable  throughout  its  defined  region.  While  such  an 


assumption  may  have  some  theoretical  implications,  it  provides  little 
or  no  practical  limitation  to  the  evaluation  of  bivariate  H-function 
densities . 

By  applying  the  residue  theorem  iteratively  to  find  the 
analytic  form  of  the  bivariate  probability  density  function,  the 
procedures  outlined  by  Eldred  (7:112-136)  and  Cook  (5:115-118)  apply 
to  the  bivariate  inversion  techniques  as  well.  However,  the 
convergence  conditions  given  by  Cook  (5:61-83)  can  not  immediately  be 
generalized  to  the  bivariate  case.  While  in  most  instances,  one  can 
invert  the  Mellin  transform  of  the  bivariate  H-function  density  by 
summing  residues  first  in  either  the  left  or  right  half  s^  plane 
followed  by  a  summation  of  residues  in  either  the  left  or  right  half 
plane,  a  set  of  bivariate  convergence  conditions  must  be 
established  before  an  analysis  of  the  analytic  form  of  the  bivariate 
H-function  density  can  be  completed. 

6.2  Complex  Analysis  in  Multiple  Dimensions:  (91:25-40) 

Definition:  A  function  f(z^,...,zk)  in  a  domain  of  its 
variables  is  analytic  if  in  some  neighborhood  of  every  point 
(z^ ' , . . . ,zk' )  of  the  domain  it  is  the  sum  of  an  absolutely  convergent 
power  series 


repeated  application  of  the  ordinary  Cauchy  formula  leads  to  the 
formula 


*<*i 


•  > 


(2ni)‘ 


dpx 

^1 


dp2 

*2“*2 


.../ 


f(pi 


.Pu) 


pk"zk 


dpk  (6.2) 


As  in  the  case  of  one  complex  variable,  and  by  Theorem  6.1, 
since  f(z)  is  analytic  in  D,  equation  (6.2)  can  be  differentiated: 


f(z) 


“1  nk 

3z.  ^..az, 

1  k 


n^ ! . • ! 
(2*i)k 


f (pjdpj.. .dpk 


C  nl+1 

Lk  (PrZl) 


(pk‘zk> 


n^+1  (6.3) 


6.3  Inversion  of  the  Bivariate  H-f unction  Integral 

The  bivariate  H-function  is  presented  in  Chapter  3  in  the  form 
of  a  Mellin  transform  inversion  integral  and  is  given  as 


jHfx.yJ 


__1 _ 

(2ui)2 


J  ;  xl<s1)x2(s2)x3(s1+s2)x 

C1  c2 


-s2  , 

y  dsids2 


(6. A) 


where  x^(s^),  x2(s2),  and  x3(si+82^  are  defined  by  (3»7),  (3.8),  and 
(3.9)  respectively.  C^  is  a  contour  in  the  complex  s^  plane  running 
from  h-i®  to  h+i®  and  is  a  contour  in  the  complex  s2  plane  running 

from  w-i®  to  w+i®.  Both  C^  and  C2  are  indented  if  necessary  to  avoid 
the  poles  of  the  integrand. 

From  Theorem  6*2,  equation  (6.4)  may  be  written  as  an  iterated 
integral  and  may  be  given  as 
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KM 


1  “8 1  “s  2 
iH[x,y]  *  5  f  x,  (s,  )x  x» ( s n )x-j ( s i +s 5 )y  ds-j  ds»  (6.5) 

(2tt1)2  C.  C„ 


Consider  first  the  Inner  integral.  From  the  H-function 
definition,  poles  of  l(di+Dis2),  i«l...M2,  I(ei+Ei(s1+s2)) ,  i«l...M3 
lie  to  the  left  of  C2  and  poles  of  r(l-ci-Cis2),  i»l...N2, 

1  (l~f ^-Fi(s^+s2)) ,  i“l...N3  lie  to  the  right  of  C2*  Poles  of  a  gamma 
function  occur  at  nonpositive  integer  values  of  its  argument.  Hence, 
the  poles  for  the  left  half  s2  plane  may  be  given  by: 


,  j  -  0,1,2,. 


ei+Ei8l+j 


,  j  *  0,1,2,... 


m 


for  the  factors 


m2  m3 

d  i(d.+D.s~)  and  ii  r  (e,+E.  (s.+s,)) 
i-1  i-1 


respectively.  Similarly,  poles  for  the  right  half  s2  plane  may  be 


given  by 


for  the  factors 


1_ci+J 

®2ij  "  ~Ct  ’  i  "  °»1*2* 


l-fi-Fisl+J 


,  j  ■  0,1,2,... 


U  r(l-c  -C.s,)  and  11  T(l-f . -F.  (s,+s,))  •• 

i-1  1  i  i-1  i 

y 

respectively. 

In  evaluating  the  residues  of  the  poles  given  above,  s,  must 
be  carried  as  a  constant  to  be  treated  as  a  variable  in  the  next 
integral.  Because  of  this,  an  analytic  form  after  inversion  in  the  s^ 
plane  must  be  presented  to  the  plane  before  an  analysis  of  the 
location  of  the  poles  in  the  s^  plane  may  be  conducted.  However,  an 
analytic  form  can  not  be  obtained  for  the  general  case  until  a  valid 
convergence  proof,  similar  to  that  given  by  Cook.,  (3),  for  the 
univariate  H-function,  is  obtained.  However,  analysis  of  the  location 
of  the  poles  in  the  plane  is  possible  given  a  specific  set  of 

parameters  and  will  be  shown  in  examples  at  the  end  of  this  section. 

While  a  complete  analysis  of  the  location  of  the  poles  in  the 
s^  plane  is  not  possible  without  a  convergence  proof,  the  poles  from 
the  x1(s1)  term  may  be  identified  and  are  given  as 

j  ”  0,1,2,... 

j  "  0,1,2,... 


bi+j 

8lij  "  “ 

and 

l-ai+j 
8lij  '  ~ aT 


for  the  factors 


respectively*  The  top  set  of  poles  are  valid  for  the  left  half  plane 
and  the  bottom  set  are  valid  for  the  right  half  plane. 


Similar  results  may  be  derived  for  2H[x,y).  2Hlx»yJ  is 
presented  in  Chapter  3  in  the  form  of  a  Mellin  transform  inversion 


integral  and  is  given  as 


1  ~si  -s2 

,H{x,y]  -  - .  /  f  x1(81)x2(-s2)x3(s1~s2)x  y  ds.ds2  (6.6) 

1  (2*i)2  Cj  C2 


where  x^(sj),  x2(-s2),  and  x3(si“s2)  are  defined  by  (3.7),  (3.4),  and 
(3.20)  respectively.  From  Theorem  6.2,  equation  (6.6)  may  be  written 


as  an  iterated  integral  and  may  be  given  as 


2H(x,y] 


1  ,  _8l  -»2 

- ,  ■*  xi(s,)x  x2(-s2)x3(s1-s2)y  ds2]  ds,  (6.7) 

(2«i)Z  Cx  C2 


Consider  first  the  inner  integral.  From  the  2^[x>y] 
definition,  poles  of  l'(di-Dis2),  i»l...M2,  I (ei+Ei(s1-s2)),  i-l...M3 
lie  to  the  right  of  C2  and  poles  of  IXl-c^+Cj^)*  i*!...^* 
r(l*fi”Fi(s1-s2)).  i"l • • .Nj  lie  to  the  left  of  C2»  Since  poles  of  a 
gamma  function  occur  at  nonpositive  integer  values  of  its  argument, 
the  poles  for  the  left  half  s2  plane  may  be  given  by: 


. 

t 

I 

» 

i 

* 

1 

j 

* 
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respectively. 

Comparing  the  poles  In  the  S2  plane  for  2H[x,y]  to  those  in 
the  82  plane  for  ^H[xty]t  it  is  clear  that  for  cases  where  ^H[x,y]  and 
2H[x,y]  have  identical  parameters  the  poles  for  both  are  exactly 
identical  except  that  they  are  in  opposite  half  planes.  Poles  in  the 
8^  plane  are  exactly  identical  for  both  jH[x,y]  and  2H[x>y]*  If  82  is 


replaced  by  -s^  in  (6.7),  Chen  the  Inversion  integral  may  be  written 
as 

2H[x,y]  -  i  x.(Sl)x  1  V  x2(s2)x3(s1+s2)y  2ds2]  ds,  (6.8) 

2  (2ni)2  Cx  1  1  C2  2  J  1  2  1 

From  (6.8),  the  poles  in  both  the  s^  and  s2  planes  are 
identical  for  ^[x.y]  and  2H (x » y ]  *  Therefore,  the  residues  of  ^H[x,yJ 
will  be  identical  to  those  of  2H[x,y]  except  that  ^H(x,y]  will  have 
variables  (x,y)  and  2H(x,y]  will  have  variables  (x,l/y).  This  could 
also  have  been  derived  directly  from  (3.13).  From  this  analysis,  it 
is  clear  that  if  a  set  of  convergence  conditions  can  be  derived  for 
^H[x,y],  then  it  will  also  be  valid  for  2H [x ,y ] • 

All  of  the  bivariate  H-function  probability  density  functions 
given  in  section  4.3  were  inverted  to  reproduce  their  analytic  forms 
using  the  techniques  outlined  above.  The  exception  to  this  is  the 
Kellogg-Barnes  I  distribution.  It  can  not  be  inverted  for  the  general 
case  because  specific  values  of  the  8  parameter  determine  different 
pole  sequences.  However,  once  6  is  set,  the  residue  sequence  may  be 
established  and  the  distribution  may  be  inverted.  The  Kellogg-Barnes 
distribution  was  inverted  for  values  of  6  *  1  and  8  *  2. 

The  following  examples  demonstrate  the  inversion  techniques 
outlined  above.  The  bivariate  gamma  and  beta  distributions  are 
demonstrated  in  Examples  6.1  and  6.2  because  they  are  the  most 
difficult  to  invert  and  because  they  demonstrate  the  validity  of  the 
inversion  process  when  f^  Y(x,y)  *8  positive  for  only  a  portion  of  the 
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posicive  quadrant.  Example  6.3  demonstrates  the  inversion  techniques 
for  2H[x,y]  variates. 

Example  6.1:  Consider  McKay's  bivariate  gamma  distribution  given  by 
p+q  i  , 

fx, *<*•*>  ■  f<f5TcT)  *  lJ"x>'1  •  y>x>0 


2  1,0, 0,0, 1,0 

rf-.  .  H 

1  1  1,1, 0,0, 0,1 


ax 

ay 


(p+q-2,1) 
(p+q-1,1)  ;  - 


(P-1,1)  ; 


l  r(p-i+s1)J (p+q-2+s1+s2)  -sl  -s2 

— —  J*  -  (ax)  (ay)  ds  ds. 

(2wi)z  IXp+q-l+Sj) 


(6.9) 


From  Theorem  6.2,  the  Mellin  transform  inversion  integral  in  (6.9)  may 
be  treated  as  an  iterated  integral  and  may  be  written  as 


fx,Y(x,y)  - 

2  .  r (p-l+s,)(ax)  1  -s2 

iTp)  2Ti;  i(p4q-l+9l)  l2Fi;  rCP+9-2+81+82)(ay)  ds2]  dSj 

(6.10) 

Consider  now  the  inner  integral  in  (6.10).  The  integral  may  be 


»  w\  A  •\'IV  f  1 


4* 


-.'■A 


evaluated  by  summing  the  residues  In  the  left  half  &2  plane.  The 
poles  of  the  Integrand  are  given  by 


s2j  “  _(P+(l_2+8i+j)  »  J  *  0,1,2,... 

Summing  over  j  the  residues  of  the  poles  given  by  and  using  the 
property 


r(x) 


1  (x+n) _ 

(x+n-l)(x+n-2). . .x 


yields 


2  residues  - 


-s 


2j 


”  j  (p+q-l>j+s1-»^2j^8l'>,s2.i4'p'fq"'2'>'J^y/a) 

j-0  (si+s2j+P+tl-2^sl+s2j+P+<l”^^*  ’  *  (®l+s2  j+P+^+j) 


s2j“"( j+P+q-2+Sj) 


«,  s^+q+j-2 

-  1  (ay)  (-D3 


j-0 


3T 


si+p+q-2 

■  (ay) 


*  (-l)j  (ay)j 

j-0 


Recognizing  that  the  summation  is  just  the  power  series  expansion  for 
e~ay,  the  sum  of  the  residues  is  now  given  by 


l  residues 


(ay)P^-2  Cay)*1  e_ay 


1 


Substituting  the  solution  for  the  inner  integral  back  into  (6.10) 
gives 


fx,Y<x’y>  *  hj)  yP+q~2  e~ay 


.  r(p-i+s,)  -s, 

W  r'drtq-l-Mj)  dSl>  <6' 


The  inversion  integral  inside  the  [  ] * s  may  be  evaluated  by  summing 
the  residues  in  the  left  half  s^  plane.  The  poles  of  the  integrand 
are  given  by 


Slj  -  -(p-l+j)  ,  j  -  0,1,2,... 


Summing  over  j  the  residues  of  the  poles  given  by  yields 


I  residues 


-  T  (p+l+Sj  j+jXp-l+Sj  j+j)(x/y) 

(p-l+s1j)(p+s1j)...(p-l+s1;j+j)r(p+q-l+s1j) 

si i-l-p-j 


Using  the  property  T(q-j)  -  T(q )/(q-l)(q-2) . . . (q-j)  and  recognizing 


1 _  .  (q-D! 

(q-l)(q-2)...(q-j)  (q-l-j)! 
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E  residues 


cP'V'P  "  (-l)j(q-l)!(x/y)^ 

“W  jlo  (q-T^TTT  j! 


But,  from  (94:1),  the  series  above  is  just  the  binomial  expansion  for 
(l-x/y)q  *  , where  for  convergence,  (x/y)^  <  1,  or  restated,  y  >  x. 


„P-lyl-P 

E  residues  ■  ~ —  (l~x/y)q 

1  (q; 


-  y2~p~<1  (yO^1 


Substituting  this  identity  back  into  (6.11)  yields 

f  (x  v).  xp_1  (y-x)q_1  e"ay 

X,Yu,y;  r(p>r<q)  iy  ; 


,  y>x>0 


which  is  identical  to  the  analytic  form  given  in  (6.9). 


Example  6.2 :  Consider  the  bivariate  Beta  distribution  given  by 


r  %  r(Pl+p2+P3)  Pl-i  p2-i  p2-i 

f x, Y^x,y >  “  r(Pl)r(p2)r(p3)  x  y  (  _x’y) 


f(p1+p2+p3)  1,0, 1,0, 0,0  x  — — —  •  — • 

F(P1)F(P2)  1  o.l, 0,1, 1,0  y  (Pl+p2+p3-2,l) 

(Px-i.i)  ;  (p2-M). 


r<p1+p2+P3>  1  ~si  r  1  f  r(s2+P2-1)y 

r(P1)r(p2)  2ii;  r(si+pr1)x  r(s1+s2+p1+p2+p3-2)ds2l 


where  x,y>0  and  x+y£l.  (6.1 

From  Theorem  6.2,  (6.12)  may  be  inverted  by  the  residue  theorem 
applied  to  the  s2  plane,  treating  as  a  constant,  followed  by 
inversion  in  the  s^  plane  by  residues.  Consider  the  inner  inversion 
integral.  The  poles  of  the  integrand  are  given  by 

s2j  m  -(P2~l+j)  >  j  *  0,1,2,... 

from  which  the  sum  of  the  residues  may  be  given  by 

Z  residues  ■ 


*  _ r(s2j+p2H-j)(s2j->-p2-H-j)y  J _ 

j»0  j"^2— ^ ^s2  j+*>2 ^  *  *  *  ( ®2  j"*"P2— j )^ ( sl+s2 j+pl+p2+p3~2 ) 


P2’1+j  * 

-  y  (-1)J 

j-0  r(si+Pl+P3'1-J) 


s2j-1_P2-J 


p2-1  -  (~y)J 

7  j-0  J!  r(8i+Pi+P3_1“j) 


Substituting  this  series  back  into  (6.12)  yields 


■'  ••  s" « 


*  1 

r(Pl+P2+P3^  P?”1  t  "  (-y)Jr(sl+Pl_1)x 

fx,Y(x»y)  “  r(Pl)r(p2)  y  l2r1/j^0  j!r(s1+p1+p3-i-j)  dsi)  (6*13) 

The  inversion  integral  inside  the  [  ] * s  may  be  inverted  by  summing  the 
residues  in  the  left  half  s^  plane.  From  the  integrand  of  the 
inversion  integral,  the  poles  of  the  integrand  by 

su  *  ~(p^~l+i )  ,  i  *  0,1,2,... 

from  which  the  sum  of  the  residues  may  be  given  by 


£  residues 


(-yJ-^s^+pjKs^+pj-l+iJx 


-s 


li 


i-0  j-0  j!(sii+Pi-1)(sii+P1)...(8ii+Pi-l+1)r(sii+Pi+p3-l-J) 


«l-p,-i 


'li^-Pl 


£  £ 
i-0  j-0 


pl“1+i  . 

x  (-1) 

i!j!  r(p3-i-j) 


(-y)j 


pl“l  •  00  (p3-l,i+j)(-x)i(-y)J 

r (P3 )  j«o  i! 

Xpi_1  -  -  (i-p3,i+j)xi  yj 

i-0  j-0  i! 


where 


r(i-p3+i+j) 

d-p3,i+j)  r(i-p3) 


From  (B.4),  the  double  series  above  may  be  represented  as  an 
Appell's  hypergeometric  function  of  two  variables.  The  sum  of  the 
residues  may  then  be  given  as 


I  residues 


Ft7T>  ;b.b*5*,y> 


where  the  parameters  b  and  b’  are  arbitrary  real  positive  constants 
and  for  convergence  x+y<l.  From  Appendix  B,  special  cases,  it  is  seen 


;t>,b' ;x,y)  *  (i-x-y) 


Substituting  these  values  back  into  (6.13)  yields 


rcpj+pj+pj)  Pl-i  p2-i  p3-i 

*  r'(Pl)r(p2)rtp3)  x  y  (w'y)  •  xyy<1 


which  agrees  with  the  analytic  form  for  the  bivariate  beta 
distribution  given  in  (6.12). 


Example  6.3:  The  Kellogg-Barnes  III  distribution  is  given  by 

fX  "  f x°  e‘(ax+exF)  ,  x,y>0 

X*Y  F(c)  a,6>0  ,  c>2 


(c,l) 

0,0, 0,1, 1,0  ax  I - ;  (1,1) 


F(c)  2  0,0, 1,0,0, 1  £ 


F(c)  2iri^"  ^2ui^  ^^c+si+s2^ax^  ^sll  ds2  (6.14) 


Inverting  first  with  respect  to  Sj,  the  integrand  has  poles  at 


s!j  *  -(s2+c+j)  ,  j  *  0,1,2,... 

Summing  over  j  the  residues  in  the  left  half  s^  plane  given  by  s^ 
yields 


Z  residues  * 


*  r(c+s1j+s2+l+j)(c+s1j+s2+j)(ax) 

j-0  (c+si  j+S2^c+sl  j+s2+1^'  *  *^c+sl  j+s2+^ 


S1 j*-(s2+c+j) 


Sn+C+j  So+C+j  . 

o  X  ^  (-l)J 


-  ac  xc  (ax)82  Z 

j-o  y 


ac  xc  (ax)  2  e"ax 


Substituting  this  back  into  (6.14)  gives 


This  agrees  with  the  analytic  form  given  in  (6.14). 


CHAPTER  7 


Conclusion  and  Recommendations 

The  main  purpose  of  this  dissertation  has  been  to  demonstrate 
a  practical  technique  for  determining  the  probability  density  function 
and  the  cumulative  distribution  function  of  products,  quotients,  or 
powers  of  two  dependent  H-function  variates.  This  has  been 
accomplished  in  section  5.2.  While  trying  to  accomplish  this  purpose, 
other  contributions  have  resulted. 

Fox  (10)  and  Subrahamian  (19)  show  how  to  find  the  probability 
density  function  for  a  simple  product  or  ratio  ot  two  dependent 
variables  or  of  pairwise  independent  variables  from  two  bivariate 
distributions  using  double  Mellin  transforms.  This  work  has  been 
extended  in  section  2.4  to  account  for  arbitrary  rational  powers  of 
the  variables.  This  section  also  includes  extensions  of  Fox's  work  to 
n  sets  of  pairwise  independent  variables.  To  facilitate  manipulations 
of  Integral  transforms,  extensions  to  the  univariate  Mellin  transform 
properties  have  been  established  and  are  presented  in  section  2.2. 

A  second  type  of  H-function  which  is  strongly  related  to  the 
first,  is  defined  in  section  3.2.  Associated  properties  for  both 
H-function  types  are  given  in  section  3.3.  While  many  applications 
can  incorporate  the  two  definitions  into  one  by  simply  allowing  one  of 
the  variables  to  be  inverted  in  the  H-function  definition,  it  is 
simpler  in  other  applications,  such  as  probability  and  statistics,  to 
keep  the  two  types  of  H-functions  separated. 
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A  remarkably  rewarding  area  of  study  has  been  in  the  area  of 
the  cumulative  distribution  function  of  an  H-function  probability 
density  function.  First,  by  section  4.4,  the  improved  Laplace 
transform  for  the  univariate  H-function  given  by  Cook  (5)  has  been 
shown  to  be  unnecessary  if  the  H-function  is  a  probability  density 
function.  Specifically,  the  Laplace  transform  of  the  univariate 
H-function  distribution  given  by  Carter  (3)  is  sufficient  and 
complete.  This  not  only  leads  to  a  simplification  of  the  form  of  the 
cumulative  distribution  function,  but  it  also  provides  great  insight 
as  to  the  range  of  values  the  parameters  of  the  H-function 
distribution  can  undertake.  Second,  the  cumulative  distribution 
function  of  a  bivariate  H-function  probability  density  function  has 
been  shown  to  be  a  bivariate  H-function.  And,  third,  the  study  of  the 
cumulative  distribution  function  has  led  to  a  formula  for  finding  the 
constant  for  the  H-function  distribution,  given  in  section  4.5. 

As  stated  in  the  first  paragraph,  the  main  thrust  of  this 
dissertation  has  been  to  find  the  density  function  of  a  product  or 
ratio  of  two  dependent  H-function  variates.  A  natural  extension  to 
this  type  of  problem  is  to  find  the  bivariate  density  function  of  two 
H-function  variates  which  are  products  or  ratios  of  pairwise 
independent  variates  from  two  or  more  bivariate  H-function 
distributions.  This  has  been  accomplished  in  section  5.3.  However, 
unlike  the  univariate  H-function  combinations,  section  5.3  shows  that 
the  bivariate  density  function  of  some  combinations  of  H-function 


variates  will  not  result  in  a  bivariate  H-function  distribution. 

The  following  recommendations  are  made  for  directions  of 
future  work  on  bivariate  H-function  distributions: 

1.  A  set  of  convergence  conditions  similar  to  those  given  by 
Cook  (5)  for  the  univariate  case  are  needed  for  the  bivariate 
H-function.  While  several  results  have  been  obtained  on  the  bivariate 
H-function  distribution,  further  significant  advancements  would  be 
difficult  without  such  a  set  of  conditions. 

2.  Use  of  the  double  Laplace  transform  to  find  the 
probability  density  function  of  the  sum  of  two  dependent  random 
variables  is  another  area  of  possible  research.  Sneddon  (106)  devotes 
a  section  to  the  double  Laplace  transform  and  provides  some  insight  to 
the  solution  of  this  type  problem. 

3.  An  extension  of  Prasad’s  theorems  (57)  to  the  bivariate 
case  would  allow  a  researcher  to  develop  a  formula  for  the  Laplace 
transform  of  the  bivariate  H-function.  Such  a  formula  could  prove 
useful  in  the  study  of  sums  of  dependent  H-function  variates. 

4.  In  Chapter  6,  the  analytic  form  of  the  bivariate 
H-function  is  analyzed  by  performing  the  contour  integration 
iteratively.  While  this  works  well  for  special  cases,  it  is  difficult 
to  utilize  this  procedure  to  write  a  general  program  to  Invert  the 
bivariate  H-function.  It  is  believed  that  if  the  location  of  the  two 


dimensional  poles  can  be  systematically  analysed,  a  more  efficient 
inversion  technique  can  be  accomplished  by  using  the  properties  of  the 


multiple  Laurent's  series  given  in  (6.1). 

5.  Given  a  set  of  convergence  conditions,  a  study  of  the 
various  methods  of  numerical  inversion  of  the  H-function  should  be 
considered  a  must.  This  work  could  include  a  study  of  the  univariate 
computational  efficiency  as  well  as  a  generalized  code  to  numerically 
invert  the  bivariate  H-function. 

6.  Another  possible  area  is  the  use  of  H-functions  to  study 
probability  density  functions  defined  over  the  entire  real  line.  This 
area  of  research  may  be  applied  to  univariate  as  well  as  the  bivariate 
distributions.  The  positive-negative  component  methods  developed  by 
Epstein  (9),  Springer  and  Thompson  (18),  and  Fox  (10)  should 
accommodate  such  usage,  particularly  for  distributions  symmetric  about 
zero. 

7.  The  application  of  bivariate  H-functions  to  the  fitting  of 
contours  to  correlated  data  is  another  unexplored  realm.  Being  the 
most  general  of  the  special  functions  of  two  variables,  the  bivariate 
H-function  appears  to  be  as  suitable  for  contour-fitting  as  it  has 
been  for  analyzing  probability  density  functions  of  products  and 
ratios  of  dependent  random  variables. 
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APPENDIX  A 


Classical  Bivariate  Probability 


Distributions 


fX,Y(x,J0 


r(pl  +  p2  +  p3>  pl_1  p2_1  ,, 
-  x  y  (1  -  x  -  y 

r(p1)r(p2)r(p3) 


x,y  >  0,  x  +  y  £  l 


Gamma 


I)  McKay 


fv  v(x,y)  - 


r<p)r<q) 


xp_1  (y  -  x)'!*1  e"3? 


y  >  x  >  0 


II)  Cherian 


Y(x,y)  -  - ^ 

1  T* /  _ 


-(x-y)  nin(x,y)  p3-l  p.-l 

-  /  z  J  (x-z)  1  (y-z 


r<p1)r(p2)r(p3)  a 


x,y  >  o 


III)  Wicksell  -  Kibble 


fX,Y(x’y)  “  p-i_. 

cp  Lr(p) 


1  e-(x  +  y)/(I  -  c)  E  (cxy)^ 

r(p)  k-0  k! T(p  +  k)(l 


x,y  >  0 


*y 


■  *■  '■ »/ 


1  u  9y  • 


>  r  r. 
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Exponential 

I)  Gumbel 

fX  Y^x,y^  "  ^(1  +  axHl  +  ay)  -  a}  exp(-x  -  y  -  axy) 

x,y  >  o 

II)  Marshall  -  Olkin 

Fx>Y(*»y)  ■  exp{-XjX  -  X2y  +  ^3®ax(x,y)} 
x,y  >  0 


Pareto 


fXYu,y).  Bi£+12t«!» 

A>  /  I  ... 


P+1 


(bx  +  ay  -  ab)1*"2 


x  >  a,  y  >  b 


Student's  t 


*  — /  ■  — -o  - 

2*/l-pZ  (f-2) 


f  )  i  +  _ L 


(f-2)(i 


,  1  -(f+2)/4 
__  (xz  -  2pxy  +  y£)  > 

-PZ)  J 


-•  <  xfy  <  - 


F-distribution 


„  2  vi 

fXfY(x,y)  -  r(v/2)v0"v/z  n  JL_ 


v.  /2 


(vi/2)-l  (v2/2)»1 

x  y 


1-0  r(v1/2)  [(1  +  Vtx  +  v2y)/v0] 


v/2 


X,y  >0,  V  -  v0  +  vt  + 


Va'.\V,V 


■fc. -at.  JsR  <l  1 
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APPENDIX  B:  Hyper geometric  Functions 


Appel  s  Functions:  <21:281-302;7:13-36) 

The  hypergeometrlc  series  of  one  variable  is  given  as 


•  r 

*  t-  L . -  . \  _  p  (a,n)(b,n)x 
2Fl(a,b,c,x)  *n^Q  (cTnYn'r- 


(B.l) 


The  symbol  (u,k),  where  u  denotes  any  number  real  or  complex, 
and  k  is  any  real  integer,  is  defined  by 


(u,k)  -  -  u(u+l)...(u+k-l),  k  >  0 

T(u)  - 


(u.-k)  -  llil  ,  k  <  0 

(l-u,k) 


Elements  a,  b,  and  c  are  the  parameters  of  the  series  and  x  is 
the  variable  of  the  series •  The  series  is  not  defined  if  c  is  a 
non-positive  integer,  unless  a  or  b  is  also  a  negative  Integer  such 
that  -c  <  -a  or  -c  <  -b. 

Appell  derived  an  expression  for  a  hypergeometrlc  function  of 
two  variables  by  considering  the  simple  product  of  two  Gauss 


functions 
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2^^ (s »b jc ; x) 2F ^ (a ’ ,b 1 5 c  jy) 


m„n 


l  £  (a,m)(a',n)(b,m)(b',n)x  y 

m«0  n-0  (c,m)(c' ,n)m!n! 


(B.2) 


and  replacing,  in  turn,  each  pair  of  products  (a,m)(a',n),  for 
example,  by  the  composite  product  (a,m+n).  By  considering  the 
possible  composite  product  combinations,  Appell  derived  the  following 
four  hypergeometric  functions  of  two  variables. 


m«n 


F.Ca.b.b'jcjx.y)  -  l  £  (a,m4n)(b^Bj(b' ,n)x  y 

1  m-0  n-0  (c ,m+n)m!n! 


(B.3) 


where  for  convergence,  x  <  1,  y  <  1 


F2<..k,b'1c,c*ix.y)  -  J  S  (B.4) 

m-0  n-0  (c,*)(c  ,n)m!n! 


where  for  convergence,  x  +  y  <  1 


F3(a,a' ,b,b*;c;x,y) 


s  j  (a ,m)(a 1  .nXb.mHb'  ,n)xmyn 
m-0  n-0  (c,m+n)m!n! 


(B-5) 


where  for  convergence, 


x|  <  1,  ly |  <1 


F4(a,b;c,c';x,y) 


j.  j.  (a,m+n)(b,m+n)xmyn 
m-0  n-0  (c ,m)(c ' ,n)m!n! 


where  for  convergence. 


(B.6) 


Special  Cases:  (12:160) 


F1(a,b,b' ;a;x,y)  -  (1  -  x)-b(l  -  y)-b 

F2(a,b,b;b,b' ;x,y)  -  (1  -  x  -  y)~a 

F2(a,b,b';a,b,;x,y)  -  (1  -  x  -  y)  b(l  -  y)b  3 

F2(a,b,b' ;b,a;x,y)  *  (1  -  x)b  a(l  "  x  -  y)  b 

Kampe'  de  Ferle't's  Function:  (7:29-33) 

Appells  functions  may  be  represented  as  special  cases  of  a 
more  generalized  double  hyper geometric  function  of  higher  order.  Thl: 
function  was  first  defined  by  Kampe  de  Feriet  in  1921  and  is  given  by 


— 

- 

A 

B 

bl»bl'»  •  •*bB»bB' 

x,y 

C 

•  •  »Cq 

D 

*^1  »^1 ' » *  * 

B 


m..n 


Z  Z 


il  (a^.m+n)  H  (b^,m)(bi' ,n)x  y 
1-1  i-i 


(B.7) 


m-0  n-0 


II  (c<tm+n)  II  (d,  ,m)(dj  '  ,n)m!n! 
j-1  J-i 


where  A  +  B<C+D  +  lotA  +  B<C+D  +  l  and 


x  +  y  <  mln(l,2C~D+l) 


3/2 


^AD-A145  598  ALGEBRAIC  FUNCTIONS  OF  H-FUNCTIONS  UITH  SPECIFIC 
DEPENOENCV  STRUCTURE(U)  AIR  FORCE  INST  OF  TECH 
MRIGHT-PATTERSON  AFB  OH  S  D  KELLOGG  HAV  84 
UNCLASSIFIED  AFIT/CI/NR-84-39D  F/G  12/1  NL 


APPENDIX  C:  Kellogg-Barnes  Distributions 


Kellogg-Barnes  Type  I_  distribution 


f  (x  .  V*1  (X 2  +  y2\S  -a(x2  +  Y2) 

rX,YU*y)  irr(B+l)  '  X  ; 


x,y>0 

u,8>0 


The  probability  density  function,  (p.d.f.),  of  X  for  8*1  is  given  as 


fx(x) 


/  (x2+y2)e~ay2  dy 

0 


<xz 


By  symmetry,  the  same  is  true  for  the  p.d.f  of  y. 
moments  of  the  distribution 


The  moments  of  the  distribution  may  be  derived  from  equations 


y  y 

(1.5)  -  (1.11).  By  symmetry,  it  is  clear  that  p  »p  and  o  *o  . 

x  y  x  y 


84-1  •  •  ?  ? 

*a  f  I  x(x^  +  y*)®  e_a^x  +  y  ^  dxdy 


x  »r(8+l)  0  0 


2  2  2 

Letting  r  «  x  +y  ,  x  ■  rcosQ,  and  y  ■  rsind  yields 


A_8+l  -  _  2  n/2 

p_  ■  *  e  /  cos6  d0dr 


x  xr(s+i) 


Performing  the  inner  integral  and  letting  z*r*  and  dza2rdr  gives 


x  nr(8+l) 


<•  c  K1C, 

0 


2IVB+3/2) 

xTorTa+T) 


The  variance  may  be  given  by 


S  1  *2fY  v(x,y)  dxdy  -  u  2 
0  0  ^  x 


4a 


0+1 


nr(B+l) 


«•  ir/2  2 

/  /  r28+3  e~«  cos2q  d3<ir  -  u  2 

0  0  X 


B+l  _  4  rr(3+3/2)1 2 
2a  w2al  T(0+1)  J 


From  (1.10),  the  covariance  of  X,Y  is  given  by 


cov(x.y)  -  /  /  xyf  v(x,y)  dxdy  - 

0  0’  x  y 

■  A  «-ar2  V-  —  *•*  -  v, 

_  B+i  _  _4_  rr(8+3/2)-i2 

na  w2a|.  r(0+l)  J 

Since  the  correlation  for  X,Y  is  the  covariance  divided  by  the 
standard  deviations  of  X  and  Y,  the  correlation  p(x,y)  may  be  given 


«-’>  -  KISS) 


where 

a  .  4  rr(g-t-3/2)-|2 
n2L  r(S+l)  J 

Kellogg-Barnes  Type  II  distribution 

£.  „(*,,)  -  8a2  .-<«  +  W"  ,  x,y>0 

**  a,B>0 

The  p.d.f.  of  X  Is  given  by 


£x(.)  -  bo,2  /*  +  By/,)  By 

-  a2  xe-ax 

Similarly,  the  p.d.f.  of  Y  is  given  by 

fY(y)  -  8a2  /"  e-<«  +  6*/x>  dy 
0 

Recognizing  that  e”ox  is  the  kernel  of  the  Laplace  transform,  the 
p.d.f.  may  be  given  by 

fy(y)  *  1-a  1  8°2«"8y/X  > 


and  from  (95:146  #25) 


fY(y)  -  2(0o)3/2  /y  2(a8y)1/2  ] 

where  Ky[z]  Is  the  modified  Bessel  function  as  defined  by  Erdelyi 
(95:371;8:5).  The  final  form  for  f^(y)  can  also  be  derived  by  taking 
the  Mellin  transform  of  f^  y(x,y)  anc*  then  taking  the  limit  as  s-»l, 
(15:27  #3.16). 
moments  of  the  distribution 

li.  -  /V  «Sa2  +  dxdy 
0  0 

-  2/o 


M  -  /V  yfio2  •-<«  4*/*>  dxdy 
7  0  0 

-  2/ofi 


o^2  ■  fQfQ  x2  So^  e”^ax  +  dxdy  -  ux 

m 

-  a2  /  x3  e-0nc  dx  -  4/a2 
0 

■  2/o^ 


Kellogg-Barnes  Type  III  distribution 


f  /  \  c  -(ax  +  3xy) 

fX  v(x»y)  “  r7^  *  * 

A*x  r(c)  a,3> 


The  p.d.f  of  X  is  given  by 


'*<«> '  fe  ‘c  «•“  ‘•6X1'  * 


a  vc-l  -ax 

r(c)  x 


Similarly,  the  p.d.f.  of  Y  is  given  by 


**w  ■  ?£>  '0 «'  '•<*f8y>x  *■ 


gqC  T(c+1) 

r(c)(o+3y)cH 


-  c3ac  (a  +  3y)"(c+1) 


moments  of  the  distribution 


„-•£//  xc+l  e-<ax  +  8xy)  dxdy 

x  oo 


APPENDIX  D:  Contour  Plots  of  the  Bivariate  H-f unction  Distribution 

The  following  contour  plots  were  accomplished  using  the  Vector 
General  3404  graphic  display  with  the  VAX  11/780  computer,  VAX/VMS 
operating  system  3.5,  Advanced  Graphics  Lab.,  University  of  Texas. 


McKay's  bivariate  gamma 


U  1  O  L  ^  lUUblUll 


f  (x  y)  *  aP  q  -  xP-1  (y-x)<*~1  e”a^ 
X,Y(x'y>  F(p)r(q)  ; 


y>x>0 

a,p,q>0 


a)  a®2.0  ,  p»q*0.5 

Figure  D.l  Contour  Plot  of  McKay's  Bivariate  Gamma  Distribution 


The  bivariate  beta  distribution 


fX  Y(x,y)  -  r(pl)r(p2)r(p3) 


r(p1+p,+p3)  pl-i  p2-i  p3-i 

- - -  x  y  (1-x-y) 


x,y>0  ,  x+y£l  ,  p1,p2,p3>0. 


r(P1+P2+P3)  1,0,1, 0,0,0 

rCPl)r<P2)  lHo,l,0,l.l,0 


( Pl+P2+P3~2 , 1 ) 

(prl,l)  ;  <p2-l,l) 


Plots  of  the  bivariate  Beta  distribution  are  shown  in  Figure  D.2. 


*>  Pl-P2-P3-0.5 

Figure  D.2  Contour  Plot  of  the  Bivariate  Beta  Distribution 


I  ff 

If  _ 

0  - 


Vv>\ 


b)  o«1.0  ,  8-1.0 


c)  a-l.O  ,  6-2.0 

Figure  D.3  Contour  Plot  of  the  Keilogg-Barnes  I  Distribution 


c)  a«2 .0  ,  0*1.0 
Figure  D.4  Contour  Plot  of  the  Kellogg -Barnes  II  Distribution 


Figure  D.S  Contour  Plot  of  the  Kellogg-Barnes  Ill  Distribution 
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